Does this output indicate my bam is paired end?
1
0
Entering edit mode
3.0 years ago
curious ▴ 750

I have bams from the same sequence batch, trying to find out if they are paired read

as per this post

samtools view -c -f 1 in.bam

output 53178051 plugged in to this decoder, which seems to say paired read because it contains 1

Oddly, the very next sample produces 47757395

I have several hundred bam from the same batch, about half of then output the 1

Wouldn't this be unexpected if they are all from the same batch?

samtools • 577 views
ADD COMMENT
1
Entering edit mode
3.0 years ago

I don't know if you quite understand what that command line is doing.

Whatever the output of that command is, that is how many paired end reads are in your bam.

ADD COMMENT
0
Entering edit mode

-> my embarrassment (-‸ლ)

ADD REPLY

Login before adding your answer.

Traffic: 4034 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