How to get CDS.fasta from gff file + sequenceFile.fasta
I have two files:
- An annotation file in gff format
- A fasta formatted sequnce file for a whole genome.
I want to have a resulting file : myGenomeCDS.fasta
How can I go about doing this via python or Tools ?
Thanks
• 9,034 views
•
link
2 answers
If your gtf file has the annotation of CDS's you can extract CDS sequences using gffread in this way:
gffread -g genome.fa -x CDS.fa annotation.gtf
• 0 views
•
link
There are so many ways to do this. If you don't want to write your own code, use the utilities from preexisiting packages such as BedTools (fastafrombed), Glimmer (extract) and even Galaxy as a tool (extract features).
• 0 views
•
link
Log in to answer this question.