bedtools multicov not running
1
1
Entering edit mode
8.6 years ago
Bnf83 ▴ 150

Hi guys, I'm trying to count the number of reads per gene (RNA Seq analysis) to perform the gene expression analysis. I'm using: bedtools multicov -bams my_file.bam -bed Human_genes.bed > gene_counts.txt Strangely, the gene_counts.txt file is exactly equal to the Human_genes.bed file. Since no Error or Warning message is generated I don't know how to deal with this issue. Does anyone have experience of this kind of error?

Thanks in advance

RNA-Seq bedtools multicov gene expression • 3.2k views
ADD COMMENT
1
Entering edit mode

Because you are giving it a bed file, it will tell you how many regions from that bed file are in the bam file and it seems like you have all of them.

ADD REPLY
0
Entering edit mode

Yes, you' re right, but as from here: http://bedtools.readthedocs.org/en/latest/content/tools/multicov.html I expect the count of alignments but this does not appear in the output file.

ADD REPLY
0
Entering edit mode

what bedtools version you are using?

ADD REPLY
0
Entering edit mode

I'm using bedtools v2.16.1

ADD REPLY
0
Entering edit mode
7.9 years ago
ariel.balter ▴ 260

I'm having the same issue. Is it a problem with the format of the input files? If you figure it out, post the answer!

ADD COMMENT

Login before adding your answer.

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