Question: Cuffdiff-output-couldnot find fasta record for a scaffold
0
gravatar for raj_262920
3.8 years ago by
raj_26292020
Germany
raj_26292020 wrote:

Dear Galaxy team

When i ran cuffdiff for my merged transcripts file using a custom reference genome, It ran successfully and i could find significant differences. But i see under the cuffdiff files, it says' couldn't find fasta record for scaffold 1 and scaffold 10, will not be considered for bias correction'.

I checked my reference genome for the respective scaffolds and they are present and i also checked the gtf file and looks like it is represented correctly. I don't know where this messages is coming from and why ?!

I think because of this it doesnt give any FPKM values for genes in those scaffolds which inturn just says NOTEST for these scaffolds. Could you please help me to possibly clear this ?

Thanks

rna-seq galaxy • 1.2k views
ADD COMMENTlink modified 3.8 years ago by Jennifer Hillman Jackson25k • written 3.8 years ago by raj_26292020
0
gravatar for Jennifer Hillman Jackson
3.8 years ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

Double check that the chromosome identifiers are an exact match between all of the inputs. Even slight differences can be problematic. 

Also double check that data for for these chromosomes actually exists in Cuffmerge output. 

More custom reference genome troubleshooting help can be found on this wiki. Formatting can be very important. It is often necessary to do this formatting as the first step, before running any tools in the analysis (before mapping). Sometimes you can adjust a bit and it will be OK, other times not. The "Troubleshooting" section has instructions for making changes, and most can be done within Galaxy.
http://wiki.galaxyproject.org/Learn/CustomGenomes

Best, Jen, Galaxy team

ADD COMMENTlink modified 3.7 years ago • written 3.8 years ago by Jennifer Hillman Jackson25k
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 16.09
Traffic: 169 users visited in the last hour