Identifying discordantly mapped reads
How can I use samtools to identify discordantly paired reads or is it even possible? I can use awk to simply output reads which are mapped on different chromosomes, and also reads where the distance between the pairs is significantly larger than expected but is there a way to do this via samtools?
2 answers
There are different options:
- You can use
samtoolsand give as argument the following flag-F 1294. In this way you'll discard all reads having which are: read mapped in proper pair, read unmapped,mate unmapped,not primary alignment,read is PCR or optical duplicate. - As you've mentioned, you can use
awkto extract mates mapped on different chr. You should convert your bam to sam, and thenawk '($3!=$7 && $7!="=")'.
I would vote for the first approach, because it includes the distance between pairs issue.
using samjdk: http://lindenb.github.io/jvarkit/SamJdk.html
$ java -jar dist/samjdk.jar -e 'return record.getReadPairedFlag() && !record.getReadUnmappedFlag() && !record.getMateUnmappedFlag() && !record.getReferenceName().equals(record.getMateReferenceName());' in.bam
• 1 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Why Is My Uniquely Mapped Read Count Higher in Samtools Compared to Bowtie2?
written by daffodil •
-
What can cause broken read-pairs (chimeric read-pairs) in a sequencing run?
written by William
-
Uniquely mapped paired-end reads on the same chromosome
written by barzilayrom1 •
-
Filtration of paired reads with both length higher than 50 from bam file
written by Mgs3 •
-
Bowtie2 Alignment Summary vs Samtools
written by ATRX
-
samtools - Finding Average Fragment size for pair?
written by oars
-
Identification of short insertions in genome
written by valerie
-
Filtering both concordantly and discordantly mapped reads
written by valerie
-
Calculate the number of paired-end reads their middle size is greater than # KB
written by lisadavic66 •
-
Concordant pairs bam file larger than accepted_hits bam file?
written by bmpbowen •
just filter out flags 3854 ( http://broadinstitute.github.io/picard/explain-flags.html ) ? (
)