gffread
gffread copied to clipboard
An error was found in the generated protein sequence
Hi @alephreish @gpertea
I extracted CDS sequences and protein sequences based on the CDS contained in the GTF.
gffread Y2.gtf -g Genome.fa -x cds.fa -y pep.fa
But I found the cds sequence in this ID is normal, but the protein sequence contains some .
(For the record: I do not maintain gffread).
The "CDS" sequence you pasted (I had to OCR it - pls do not paste sequence screenshots!!) has multiple stop-codons. Stops are indicated by dots by default in gffread (check out the -S argument).
@alephreish
Okay, thanks.