Question: genomescope profile failed to converge
0
gravatar for StudentBio
5 months ago by
StudentBio0
StudentBio0 wrote:

hello the community,

to estimate the characteristics of my genome I used GenomeScope, as a result I got an error like GenomeScope profile failed to converge and depending on what i found the problem is because low coverage or low quality reads ....so please if anyone can help me to solve this this my fastqc report after trimming enter image description here

fastqc report after trimming with sickle

jellyfish k-mers genomescope • 334 views
ADD COMMENTlink modified 5 months ago by lieven.sterck3.3k • written 5 months ago by StudentBio0

I assume you did use the jellyfish tool to count the kmer? If oyu used any other software (ntCard ?) you need to modify the input count table before uploading it to genomescope

It might also depend on the kmer size you used, so try to play around with that one a bit and see if it gets better?

ADD REPLYlink modified 5 months ago • written 5 months ago by lieven.sterck3.3k

I removed the -C option from jellyfish and it worked thank you very much jellyfish count -m 17 -s 1000000000 -t 10 /home/vshare/bowtie2/trimmed.fastq -o reads2.jf

ADD REPLYlink written 5 months ago by StudentBio0
1
gravatar for lieven.sterck
5 months ago by
lieven.sterck3.3k
VIB, Ghent, Belgium
lieven.sterck3.3k wrote:

I had this myself as well on several occasions. it's often due to the fact that you have too low coverage of data and that you thus not see the typical peaks appearing (yet at that coverage) . You should anyhow get a graph, no? on that you will see that you're lacking the peak the software needs to do it's calculations and thus not converges.

EDIT: is also should not have anything to do with FastQC processing, it should work as good (or bad) on trimmed or untrimmed data

ADD COMMENTlink modified 5 months ago • written 5 months ago by lieven.sterck3.3k

Same for me, lack of coverage so no peak, I've never seen anything else myself.

There's an unlikely chance that it's a highly fragmented genome, see Suppl. Figure 3 in 'The Plastid Genome in Cladophorales Green Algae Is Encoded by Hairpin Chromosomes'. But again, unlikely.

ADD REPLYlink written 5 months ago by Philipp Bayer5.8k

small world ... paper of good colleagues of me :)

ADD REPLYlink written 5 months ago by lieven.sterck3.3k

Tell them 'congratulations' from me, it's a really cool paper!

ADD REPLYlink written 5 months ago by Philipp Bayer5.8k

thanks of you all i can resolve it removing the -C option from jellyfish and it worked jellyfish count -m 17 -s 1000000000 -t 10 /home/vshare/bowtie2/trimmed.fastq -o reads2.jf

ADD REPLYlink written 5 months ago by StudentBio0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1225 users visited in the last hour