Question: gffread "Error (GFaSeqGet): end coordinate (76134) cannot be larger than sequence length 76132"
0
gravatar for qwzhang0601
3.4 years ago by
qwzhang060170
United States
qwzhang060170 wrote:

Hello: I am trying to get coding sequences based on gff file and the genome fast file. I am using the function gffread. It works well with "gffread my.gff -g genome.assembly.fasta -x cds.fa". But after I add the parameter "-J" it reported an error "Error (GFaSeqGet): end coordinate (76134) cannot be larger than sequence length 76132". But I manually checked the annotation of the contig with the length 76132, and found there is no annotation with the end coordinate "76134". The maximum coordinate is 76132.

About the "-J" parameter

-J discard any mRNAs that either lack initial START codon or the terminal STOP codon, or have an in-frame stop codon (only print mRNAs with a fulll, valid CDS)

Anybody have met the same problem?

Thanks

gffread cufflink gff • 2.6k views
ADD COMMENTlink modified 3.3 years ago by elaak910 • written 3.4 years ago by qwzhang060170

Hi. I have the same issue. Did you find any solution?

ADD REPLYlink written 3.3 years ago by elaak910

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized.

ADD REPLYlink written 3.3 years ago by genomax89k

I had the same problem... For me the "easiest" way is to do it is to run gffread without the -J parameter. Then with a custom script check that each sequence with "ATG" and finish with one stop-codon. Additionally, you have to look for stop codons inside the sequence.

ADD REPLYlink written 15 days ago by JC0
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: 2045 users visited in the last hour