Off topic:Quickest Way To Pull A List Of Reads Out Of A Bam File?
1
0
Entering edit mode
9.8 years ago

What is the quickest way find the final aligned locations of a list of reads in BAM file? I can call

samtools view my.bam | grep -e uniqueID1 -e uniqueID2 > output

to grep for uniqueID1 and uniqueID2 in my.bam, but crawling through the SAM is quite slow. I'd like to call:

handyTool extractID my.bam -l listOfIdentifiers > output

and get the SAM line for each identifier found in the BAM. Is there a better way?

samtools next-gen • 4.5k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 1820 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6