Question: bedtools multicov displays only the last BED feature
0
gravatar for ThePresident
22 months ago by
ThePresident110
ThePresident110 wrote:

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

multicov bedtools • 890 views
ADD COMMENTlink modified 27 days ago by rania.hamdy10 • written 22 months ago by ThePresident110

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 REPLYlink written 22 months ago by ThePresident110

Was it made on Windows?

ADD REPLYlink written 22 months ago by WouterDeCoster39k

It was made on mac OS

ADD REPLYlink written 22 months ago by ThePresident110

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 REPLYlink written 27 days ago by rania.hamdy10
0
gravatar for rania.hamdy1
27 days ago by
rania.hamdy10 wrote:

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 COMMENTlink modified 27 days ago • written 27 days ago by rania.hamdy10
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1036 users visited in the last hour