Is there a way to remove reads whose mate does not map as expected in a Bam file directly? I know we can run command line "samtools view -h -f 2 in.bam out.sam" to exclude the reads in the output sam file.
However, I wonder if I could exclude those reads directly from the in.bam file.
•
written
by
devinliao0918 •
1 answer
In other words, you want to know if it's possible to edit a BAM file in-place with samtools. The answer is no (well, samtools will let you by just specifying the same file for input and output). While it's theoretically possible, it'd be a good way to shoot yourself in the foot if not done correctly (such as if the writer file pointer caught up to and overwrote where the reader file pointer was pointed) so it's not something anyone will recommend you frequently do.
•
written
by
Devon Ryan
Log in to answer this question.
-
Samtools merge bam issue with header tags
written by sam •
-
Separating mapped and unmapped reads from BAM file
written by marongiu.luigi
-
Extract reads where both pairs maps properly and none of the pairs are split from a bam file
written by rimjhim.roy.ch
-
How do you remove reads from .bam file with FLAG = 0?
written by josolvason •
-
RNA-Seq Analysis: To sort or not to Sort?
written by kanika.151
-
Mapping shorts reads using Bowtie
written by bioinfo
-
mapped/unmapped mate with samtools
written by Flavia •
-
Samtools: remove reads that are masked as duplicates
written by devinliao0918 •
-
How to use Samtools to remove reads whose mate does not map as expected
written by devinliao0918 •
-
Extracting reads with low mapping score from Bam file
written by devinliao0918 •
Hi, Take some time to go through the samtools manual. It will definitely help you. Samtools works for both sam and bam format and can generate output files in both sam and bam format. Also, you can't just change the original file with samtools. You will have to create a new sam or bam file using the original sam/bam file. The command mentioned above can be used to generate a new "bam" output if
-bparameter is used. Check http://samtools.sourceforge.net/samtools.shtml.Thanks. I'm not sure whether "-b" works but will definitely try that.
I tried to read samtools manual carefully but failed to understand all of the details. I would like to see more examples.
It will work. But there will be a new bam file instead of modifying your original bam file.