Skip to content

Commit e132f20

Browse files
committed
Update Bismark arguments
1 parent 09433e2 commit e132f20

1 file changed

Lines changed: 5 additions & 26 deletions

File tree

primeseq/src/org/labkey/primeseq/pipeline/BismarkWrapper.java

Lines changed: 5 additions & 26 deletions
Original file line numberDiff line numberDiff line change
@@ -136,9 +136,8 @@ public AlignmentOutput performAlignment(Readset rs, List<File> inputFastqs1, @Nu
136136

137137
args.add("--samtools_path");
138138
args.add(new SamtoolsRunner(getPipelineCtx().getLogger()).getSamtoolsPath().getParentFile().getPath());
139-
args.add("--path_to_bowtie");
140139

141-
//TODO: consider param for bowtie vs. bowtie2
140+
args.add("--path_to_aligner");
142141
args.add(getBowtie2Exe().getParentFile().getPath());
143142

144143
if (getClientCommandArgs() != null)
@@ -247,11 +246,11 @@ public IndexOutput createIndex(ReferenceGenome referenceGenome, File outputDir)
247246
}
248247
}
249248

250-
//first build for bowtie2
249+
// build for bowtie2
251250
List<String> args = new ArrayList<>();
252251
args.add(getWrapper().getBuildExe().getPath());
253252
args.add("--bowtie2");
254-
args.add("--path_to_bowtie");
253+
args.add("--path_to_aligner");
255254
args.add(getBowtie2Exe().getParentFile().getPath());
256255
args.add(indexOutputDir.getPath());
257256
getWrapper().execute(args);
@@ -263,21 +262,6 @@ public IndexOutput createIndex(ReferenceGenome referenceGenome, File outputDir)
263262
throw new PipelineJobException("Unable to find file, expected: " + bowtie2TestFile.getPath());
264263
}
265264

266-
//then build for bowtie
267-
List<String> args2 = new ArrayList<>();
268-
args2.add(getWrapper().getBuildExe().getPath());
269-
args2.add("--bowtie1");
270-
args2.add("--path_to_bowtie");
271-
args2.add(getBowtieExe().getParentFile().getPath());
272-
args2.add(indexOutputDir.getPath());
273-
getWrapper().execute(args2);
274-
275-
File bowtieTestFile = new File(genomeBuild, "CT_conversion/BS_CT.1.ebwt");
276-
if (!bowtieTestFile.exists())
277-
{
278-
throw new PipelineJobException("Unable to find file, expected: " + bowtieTestFile.getPath());
279-
}
280-
281265
File indexBaseDir = new File(localFasta.getParentFile(), getIndexCachedDirName(getPipelineCtx().getJob()));
282266
if (!indexBaseDir.exists())
283267
{
@@ -307,8 +291,8 @@ public static class Provider extends AbstractAlignmentStepProvider<AlignmentStep
307291
{
308292
public Provider()
309293
{
310-
super("Bismark", "Bismark is a tool to map bisulfite converted sequence reads and determine cytosine methylation states. It will use bowtie for the alignment itself.", Arrays.asList(
311-
ToolParameterDescriptor.createCommandLineParam(CommandLineParam.create("-L"), "seed_length", "Seed Length", "Sets the length of the seed substrings to align during multiseed alignment. Smaller values make alignment slower but more sensitive.", "ldk-numberfield", null, 30),
294+
super("Bismark", "Bismark is a tool to map bisulfite converted sequence reads and determine cytosine methylation states. It will use bowtie2 for the alignment itself.", Arrays.asList(
295+
ToolParameterDescriptor.createCommandLineParam(CommandLineParam.create("-L"), "seed_length", "Seed Length", "Sets the length of the seed substrings to align during multi-seed alignment. Smaller values make alignment slower but more sensitive.", "ldk-numberfield", null, 30),
312296
ToolParameterDescriptor.createCommandLineParam(CommandLineParam.create("-N"), "max_seed_mismatches", "Max Seed Mismatches", "Sets the number of mismatches to be allowed in a seed alignment during multiseed alignment. Can be set to 0 or 1. Setting this higher makes alignment slower (often much slower) but increases sensitivity. Default: 0.", "ldk-numberfield", new JSONObject(){{
313297
put("minValue", 0);
314298
put("maxValue", 1);
@@ -844,9 +828,4 @@ private static File getBowtie2Exe()
844828
{
845829
return SequencePipelineService.get().getExeForPackage("BOWTIE2PATH", "bowtie2");
846830
}
847-
848-
private static File getBowtieExe()
849-
{
850-
return SequencePipelineService.get().getExeForPackage("BOWTIEPATH", "bowtie");
851-
}
852831
}

0 commit comments

Comments
 (0)