Entering edit mode
11.6 years ago
Nicolas Rosewick
11k
Hi,
Is there a rapid way to extract alignment from a sam file using read ids (about~100 read ids in average). If the read ids are in a file (one per line), I could do :
cat in.sam | grep -f idFile.txt > out.sam
but with a big sam file (~40Gb) it takes a lot of time.... so is there maybe a method to extract these alignments faster ?
Thanks,
N.
well, not really duplicate. It was BAM, not SAM.
duplicate of
Extracting subsets of reads from a BAM file