View Single Post
Old 10-28-2010, 04:40 AM   #2
drio
Senior Member
 
Location: 4117'49"N / 24'42"E

Join Date: Oct 2008
Posts: 323
Default

Study carefully the SAM specification, particularly the flag field. Once you know
what flag to use, go ahead and use the filter capabilities in samtools view (-f/-F INT)
to extract the reads you are interested on.

For example (show me reads where the mate is ummapped):
Code:
$  samtools view -f 8 foo.bam
__________________
-drd
drio is offline   Reply With Quote