MGSE
MGSE copied to clipboard
estimated genome size is not being displayed
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