This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Adding CDS features to Gff3 file

Hello,

I would like to add the CDS feature to my Gff3 file, that has exons in it. I saw that genometools can do this for you

gt cds -startcodon -finalstopcodon -seqfile foo.fasta -o foo.gff3 [GFF3-FILE]

I am looking for the CDS to be the longest ORF with an ATG start codon per mRNA. Will the above command give me that? Secondly I did not see a flag for specifying the genetic code. Is this by default "standard vertebrate"?

Thanks for the help
Abhijit

genome-annotation
• 6,331 views • link
updated by Ram Ram • written by abhijit.synl abhijit.synl

Though not directly but can be acheived using a combination of Genometools and PASA

• 0 views • link
written by abhijit.synl abhijit.synl

This is resolved. please close thread

• 0 views • link
written by abhijit.synl abhijit.synl

hi，i had the same problem recently, can you tell me the solution detail

• 0 views • link
written by pigeon0411 • pigeon0411

thanks for your prompt reply, these help me a lot. but these still had a problem, after i processed the data accroding your method. the problem is : "gt cds: error file is not sorted." i think that may cause by absent some scaffold in my gff3 file, gt cds need all scaffold according to the order. gt cds will thought the file is not sorted, if some scaffold miss, even gff3 file had been sorted before. what's your opinion? thanks again!

• 0 views • link
written by pigeon0411 • pigeon0411

Please attach a sample of your Gff3 file. Maybe the first 250 lines.

• 0 views • link
written by abhijit.synl abhijit.synl

i'm sorry i don't known how attach a file in this site, can you give me your email?

• 0 views • link
written by pigeon0411 • pigeon0411

abhijit.synl & pigeon0411 :

Please do not add answers unless you're answering the top level questions. You should have been using Add Comment and Add Reply appropriately.

I'm moving the extraneous "answers" to comments, but please be more careful in the future.

• 0 views • link
written by Ram Ram

3 answers

You will need to install the following and make sure they are in your PATH variable

1) Gffread utility 2) Genometools 3) GAAS Toolkit

a) Tidy your GFF3

gffread -C -F -T --no-pseudo -o FOO1".gtf" FOO.gff3
gt gtf_to_gff3 -tidy -force -o FOO".gff3" FOO1".gtf" (gt = GenomeTools executable)
gt gff3 -sort -tidy -retainids -addintrons -setsource "Genebank" -force -o FOO"_modi.gff3" FOO".gff3"

b) Keep Longest Isoform

gff3_sp_keep_longest_isoform.pl --gff FOO"_modi.gff3" --output FOO"_modi_longest.gff3" gff3_sp_keep_longest_isoform.pl = From GAAS Toolkit)

Hope this helps. This was a simpler solution than using PASA as mentioned in my previous post Abhijit

• 0 views • link
written by abhijit.synl abhijit.synl

hi，i find GAG-Genome Annotation Generator can add startcodon&stopcodon to gff file, thanks for everyone!

• 0 views • link
written by pigeon0411 • pigeon0411

An easy way to do so would be to use TransDecoder

• 0 views • link
written by Juke34 Juke34

Log in to answer this question.

More posts like this