Extract single transcript sequences from genome assembly
Hi,
I want to extract one transcript per gene from the genome assembly and its GFF3 file registered in Refseq. And I want to extract a transcript that has a longest or a certain frame. Anyone know of an efficient method?
I am looking for this method to do ortholog comparisons.
Best regards,
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Tool: ORFanage: by-reference protein annotation and comparison for transcriptome assembly
written by Ales •
-
Restricting ncbi-genome-download (bacteria, viral, fungi, archaea, protozoa) to a specif host
written by saruman •
-
Tools to check the length of isoforms in reference transcript
written by shinyjj •
-
How to remove bases out of codon frame
written by kaz0816 •
-
Reference of Whole exome sequecing analysis
written by Shahzad •
-
How to create custom gtf annotation file?
written by John
-
RSEM error preparing index file with mouse mm10 genome reference files and Refseq
written by torihikitorti •
-
How to make a custom gtf or gff3 file with a list containing transcript loci
written by donger1103 •
-
script to Extract longest isoform from stringTie assembly
written by Bioinfonext
-
How to get the protein sequence of a transcript variant (Ensembl Perl API) ?
written by nathan.alary •
Python and biopython can help you on that
https://biopython.org/wiki/GFF_Parsing
https://pypi.org/project/gff3/
https://pythonhosted.org/gffutils/