Entering edit mode
9.7 years ago
kaston
▴
40
I'm using the bedtools genomecov command to find coverage stats for a human genome alignment:
bedtools genomecov -ibam foo.bam
The bedtools documentation says that you need to specify a .genome file so the tool knows how big the chromosomes are using the "-g filename.genome" option. I haven't done this, but the command still works.
How does it know what the chromosome sizes are? Is it pointing to a .genome file somewhere in the system that I'm not aware of or is this simply not necessary contrary to the documentation?
Thanks
Please add your reply to a comment to the answer.