got wrong reads counts from bedtools multicov(microRNA)
0
0
Entering edit mode
9.8 years ago
xdkong • 0

I used bedtools multicov to count the number of reads in the regions of a bed file. But I think the results were wrong, like:

1   17369     17436     0
1   17409     17431     0
1   17369     17391     0
1   30366     30503     0
1   30438     30458     0
1   567705    567793    4
1   567762    567783    2
1   1102484   1102578   8
1   1102504   1102525   4
1   1102540   1102561   4

I get another number by samtools view -c, e.g., I ran the command: samtools view -c mirna.sort.bam 1:1102504-1102525 The result was 444, ran samtools view -c mirna.sort.bam 1:1102540-1102561 and got 2310, I don't know why bedtools got so small numbers, who know the reason?

Thank you

RNA-Seq next-gen-sequencing alignment • 2.0k views
ADD COMMENT
0
Entering edit mode

Who has used the bedtools multicov to count reads?

I found that add the parameter -D can get a number near the number got from samtools -c.

ADD REPLY

Login before adding your answer.

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