Hi Oskar,
I got some problems about the depth calculation. I saw that BEDtools is used as the calculator for depth, and there is a perl script to calculate average depth. Previous calculation showed that for some of my datasets, in the assembly_contigs.depth.txt the merged average depth was 0.00000000. Then I manually checked the depth calculation using the BEDtools and perl script:
This is the output of BEDtools taking one sample's bam file. The assembler used was megahit, the value are generally low. So is this the reason that the output would be too low that by rounding to 8 decimal places the average depth is hence 0.00000000?
So Sorry for bothering again.
Andy
Hi Oskar,
I got some problems about the depth calculation. I saw that BEDtools is used as the calculator for depth, and there is a perl script to calculate average depth. Previous calculation showed that for some of my datasets, in the assembly_contigs.depth.txt the merged average depth was 0.00000000. Then I manually checked the depth calculation using the BEDtools and perl script:
This is the output of BEDtools taking one sample's bam file. The assembler used was megahit, the value are generally low. So is this the reason that the output would be too low that by rounding to 8 decimal places the average depth is hence 0.00000000?
So Sorry for bothering again.
Andy