Get discordant read pairs from bam
Hi,
I would like to know if there is a way to isolate discordant read pairs from a bam file. For example, when I use samtools stats in my bams, I can see the number of outward oriented read pairs and I would like to work just with them, is there a way to do that with samtools?
Thanks
• 1,887 views
•
link
updated
by
benformatics
•
written
by
paulaotero.sanchez •
•
written
by
paulaotero.sanchez •
1 answer
I found using samflags can get a bit confusing depending on the specific aligner-added flags and downstream BAM processing (i.e. marking duplicates).
https://github.com/GregoryFaust/samblaster has worked well for me in the past
• 0 views
•
link
written
by
benformatics
Log in to answer this question.
More posts like this
-
Sorting the Input Bam file by Read Id.
written by jigarnt •
-
Add or Simulate Noise in Short-Read Paired End Data
written by Suraj •
-
How to find average read length and number of reads in a bam file using samtools in linux?
written by MobiusT •
-
Confused about read pair orientation in bam
written by malthesr •
-
samtools stats result different from samtools flagstat
written by snishtala03 •
-
Extract reads based on pair orientation
written by Emilie •
-
samtools - Finding Average Fragment size for pair?
written by oars
-
Inward oriented pairs and outward oriented pairs from a BAM file - Bamcheck vs Samtools
written by Prakki Rama
-
programs on BAM file
written by zhoub •
-
Is There An Elegant Way To Extract Only The Properly-Paired Reads In A Sam/Bam File?
written by KCC
Identifying discordantly mapped reads