How to count mapping rates for ChIP-seq?
1
0
Entering edit mode
2.4 years ago

Hi, so how should we count for mapping rates for ChIP-seq? I used "samtools flagstat", but it always showed 100% mapping. For example, I have a fastq file, there are 136,192,060 lines in the file, so there should be 136,192,060/4=34,048,015 reads in this file. After mapping the reads to reference genome, Iused 'samtools flagstat' to see the mapping rate, it showed "11,892,792 + 0 mapped (100.00% : N/A)". But the mapped reads(11,892,792) is far lower than the original reads(34,048,015), how could the mapping rate be 100%? Thanks! And also, after trim-galore, the QC result showed the total sequences of the file is 9,477,948. So actually where does this number come from? Thanks so much!

Chip-seq mapping reads • 1.5k views
ADD COMMENT
0
Entering edit mode

Hello. What did you use to map the reads? I use Bowtie2, the log file of which gives useful information.

ADD REPLY
0
Entering edit mode

yes, you are right. Thanks!

ADD REPLY
0
Entering edit mode
2.4 years ago

It all depends on how the file was filtered.

Various alignments may be filtered out - for example, unmapped reads are often removed or not even reported if you pass certain flags to the aligner.

When the unmapped reads are filtered out then the alignment rate becomes 100% :-)

ADD COMMENT
0
Entering edit mode

I see, thanks so much! So how do you usually count the mapping rate?

ADD REPLY
0
Entering edit mode

It is the number of primary alignments divided by the number of reads.

ADD REPLY
0
Entering edit mode

Thanks! I mean what kind of tools you use to calculate this?

ADD REPLY
0
Entering edit mode

for example with samtools you would need to filter out unmapped (4), supplementary (2048) and secondary alignments (256) like so

samtools view -c -F 4 -F 2048 -F 256   alignments.bam

will produce a number of primary alignments

ADD REPLY

Login before adding your answer.

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