When making read extraction with regions in BAM file, also mates are considered. What is the most efficient way also to collect them?
•
written
by
user191006 •
Just query-sort the files (samtools sort -n). You don't need to randomly access anything then.
This is not what I want, sorting is not the problem for me. I focused on mates, my aim is to collect all of the mates of reads, thanks for response.
By sorting on read names mates will be following each other in the sam file.
I am going over regions in BED file. I just focus on reads in the region at a time, I try to collect their mates even if they are not in the region.
If the mates aren't terribly close then you'll end up doing a lot of linear searches even if you get both mates into the same bgzf block. I suspect it'll be faster to query sort and search for overlaps with the BED file.
I would suggest performing your alignment step using BWA within the Speed-Seq suite. This will automatically extract discordant mate pairs (that map to different chromosomes, for example) in a separate bam file for you, which will enable you to look at them without having to go through a huge bam file of "normal" reads. You can load these into IGV or another browser and look at your discordant reads without having all of the "normal" reads bogging you down.
Speed-Seq can be downloaded here: https://github.com/hall-lab/speedseq.
Here is a simple example of how to run Speed-Seq alignment that will result in the files I just described: https://github.com/samuelwb/tumor-evolution/blob/master/Alignment/SpeedSeqAlign.sh.
Log in to answer this question.
Hi,
May I ask have you solve the problem? I also encounter the same problem. I want to extract the mate reads of specific reads even if the mate is unmapped. Actually, I have posted a question at How to extract unmapped mate for specific reads using pysam or other module in python?.
Thank you!
Best regards,
Xi Zeng