bedtools multicov displays only the last BED feature
1
0
Entering edit mode
7.0 years ago
ThePresident ▴ 180

I want to count the number of reads that fall in some predefined regions of a genome. Normally, this can be easily done with bedtools multicov. I have sorted and indexed bam files of my alignments. I created a BED file with my regions of interest (excel, then saved as tab-delimited, and added .bed extension).

Here's what BED file look like:

Chr1    500000  505001  
Chr1    1000000 1001001 
Chr1    2000000 2000251 
Chr1    2500000 2500101

I run bedtools with this command:

bedtools multicov -bams ./*.bam -bed bedfile.bed

I have 5 bam files in the same directory, as well as 5 corresponding .bai index files. I get read counts for all 5 bam files but only for the last feature of the bed file, whatever the last line is. If I add more features it will still give me the counts for the last one only. Any idea what's wrong?

Thanks, TP

bedtools multicov • 2.7k views
ADD COMMENT
0
Entering edit mode

I made another BED file in a text editor and this time it worked. I have no idea what was wrong with my BED file...

ADD REPLY
0
Entering edit mode

Was it made on Windows?

ADD REPLY
0
Entering edit mode

It was made on mac OS

ADD REPLY
0
Entering edit mode

I want to count histone modification reads that fall in some predefined regions of a genome using bedtools multicov command i cant understand what BED file refer to How can I know my regions of interest?? on what basis i will create it ? can you help me please?

ADD REPLY
0
Entering edit mode
5.2 years ago

I want to count histone modification reads that fall in some predefined regions of a genome using bedtools multicov command i cant understand what BED file refer to How can I know my regions of interest?? on what basis i will create it ? can you help me please?

ADD COMMENT

Login before adding your answer.

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