On Wed, Sep 9, 2020 at 11:23 AM Alexis Groppi ***@***.***> wrote: I have run Genomescope2 online with the default parameters My starting material is : 16 fastq files from 16 SMRTcells (PacBio RSII) I have run the following command lines : jellyfish count -C -m 21 -s 12000000000 -t 20 ./*.fastq -o myreads.jf jellyfish histo -t 20 myreads.jf > myreads.histo The genome size is around 220 Mb and very homozygous But in GenomeScope2 the results are : Genome Haploid Length 16,296,648 bp 16,748,017 bp Homozygous (aa) 95.0821% 95.7105% Heterozygous (ab) 4.28954% 4.91792% Full results are here : http://qb.cshl.edu/genomescope/genomescope2.0/analysis.php?code=RAe5oDd95jhgUdK9zZ8f <http://url> How can it be possible ? Thanks for your help — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#43>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AABP347GPMZXMQKBPHEGJCTSE6MWXANCNFSM4RCRHMJA> .

On Thu, Sep 10, 2020 at 5:02 AM Alexis Groppi ***@***.***> wrote: Thanks for your answer; I looked in the maze of the Falcon pipeline I used, and I found the corrected reads in fasta format my_falcon_directory/1-preads_ovl/db2falcon/preads4falcon.fasta It worked perfectly Heterozygous (ab) 0.211677% 0.250628% Genome Haploid Length 235,933,558 bp 236,368,387 bp — You are receiving this because you commented. Reply to this email directly, view it on GitHub <#43 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AABP34Y5WJDHW77HOR7THDLSFCI4DANCNFSM4RCRHMJA> .