View Single Post
Old 08-24-2011, 01:32 PM   #3
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

bwa's sampe adds an XT tag to the sam file that indicates if the read is uniquely mapped, or repetitively mapped. So you could pull all of those reads out with awk, get the cluster coordinates, and then pull out all the reads from the .bam that had those coordinates.
swbarnes2 is offline   Reply With Quote