Get discordant read pairs from bam
1
0
Entering edit mode
9 months ago

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

samtools NGS • 666 views
ADD COMMENT
1
Entering edit mode
ADD REPLY
1
Entering edit mode
9 months ago

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

ADD COMMENT

Login before adding your answer.

Traffic: 3340 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6