View Single Post
Old 12-09-2016, 06:27 AM   #8
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

Oh! Yes, technically is not part of the BBMap package. Rather, it gets produced by BBMap if you add the flag "", because it is a custom script specific to the name of the sam/bam file BBMap produced. simply calls Samtools with the correct (or, in the above case, apparently incorrect) commands to produce a sorted, indexed bam file.

Generally, you do this:

Code: in=reads.fq.gz out=mapped.bam; sh
That produces sorted_mapped.bam and sorted_mapped.bam.bai.

Last edited by Brian Bushnell; 12-09-2016 at 09:06 AM.
Brian Bushnell is offline   Reply With Quote