schatzlab / genomescope

Fast genome analysis from unassembled short reads
Apache License 2.0
258 stars 56 forks source link

smaller genome size and lower heterozygosity #112

Open Xoey0314 opened 1 year ago

Xoey0314 commented 1 year ago

Hi,

I used Illumina sequencing data to assess the genome size. First, I used jellyfish to get the histo data. The I utilized the Genomescope for evaluation and visualization. But the results seem to have some problems.

The size of the genome in the same species is about 900Mb to 1Gb, and the heterozygosity is about 2% to 3%. However, my genome is about 700Mb, and heterozygosity is only 0.23% .

I don't know where the problem is.

My input commands are as follows: zcat 1.fq.gz ; zcat 2.fq.gz |jellyfish count -m 21 -s 1000000000 -t 12 -C -o sample.jf /dev/fd/0

jellyfish histo -t 12 sample.jf > sample.histo