estimated genome size is not being displayed #2
Description
Hi,
I've been trying to use MGSE to estimate the genome size of my species. Although I've followed the instructions to the program, it hasn't calculated the genome size for me. This is the output of the report.txt file
average coverage in reference regions (median): 0
average coverage in reference regions (mean): 0
total coverage (combined length of all mapped reads): 21971396415.0
total sequence length: 255888427
estimated genome size based on mean [Mbp]: 0
estimated genome size based on median [Mbp]: 0
This is the command I've run:
python2.7 MGSE.py --bam /home/buitracn/Genomes/corals/Pocillopora.verrucosa/01_Pver_lib4_DISCOVAR_genome_assembly_final/Pver_lib4_DISCOVAR_Assembly_Filtering/09_MGSE_GenomeSizeEstimation/Pver_alignment2haploidGenome_sorted_local.bam --bam_is_sorted --out /home/buitracn/Genomes/corals/Pocillopora.verrucosa/01_Pver_lib4_DISCOVAR_genome_assembly_final/Pver_lib4_DISCOVAR_Assembly_Filtering/09_MGSE_GenomeSizeEstimation --busco /home/buitracn/Genomes/corals/Pocillopora.verrucosa/01_Pver_lib4_DISCOVAR_genome_assembly_final/Pver_lib4_DISCOVAR_Assembly_Filtering/09_MGSE_GenomeSizeEstimation/full_table_6_Pver_lib4_filtered_ScaffCSAR_Gapfillep_A_ref_D_augustus_fly.tsv --name Pver
If you could please help me figure out what I've done wrong I'll greatly appreciate
Sincerely,
Carol