Skip to content
Snippets Groups Projects
Commit ceaa8fd0 authored by Wai Yi Leung's avatar Wai Yi Leung
Browse files

Handle RNA seq files, sam filter only takes the primary alignment

parent b4a13e74
No related branches found
No related tags found
No related merge requests found
......@@ -69,7 +69,7 @@ class Gears(val root: Configurable) extends QScript with SummaryQScript {
// sambamba view -f bam -F "unmapped or mate_is_unmapped" <alnFile> > <extracted.bam>
val samFilterUnmapped = new SambambaView(qscript)
samFilterUnmapped.input = bamfile
samFilterUnmapped.filter = Some("unmapped or mate_is_unmapped")
samFilterUnmapped.filter = Some("(unmapped or mate_is_unmapped) and not (secondary_alignment)")
samFilterUnmapped.output = new File(outputDir, s"$outputName.unmapped.bam")
samFilterUnmapped.isIntermediate = false
add(samFilterUnmapped)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment