View Single Post
Old 05-12-2015, 04:24 AM   #14
Senior Member
Location: Vienna

Join Date: Oct 2011
Posts: 123

You can use samtools view to extract the mapped/unmapped reads by filtering the 'unmapped' flag:

samtools view -F 4 -bh lib4seq.sorted.bam > lib4seq.sorted.mapped.bam
samtools view -f 4 -bh lib4seq.sorted.bam > lib4seq.sorted.unmapped.bam
Samtools view will help a lot; just have a look at some tutorials, for instance Dave's wiki
Michael.Ante is offline   Reply With Quote