I'm mapping a virus from a large bam that is mostly host sequence. Bowtie2 outputs ALL of the reads, not just those that aligned to the reference sequence.
Is there a way to get bowtie to *directly* output only mapped reads, or do I have to post-process them through Samtools?
Thanks!
Is there a way to get bowtie to *directly* output only mapped reads, or do I have to post-process them through Samtools?
Thanks!
Comment