To retrieve the transcript from gff file with the longest isoform
1
0
Entering edit mode
Aneesha • 0
@9e2ddcf8
Last seen 12 weeks ago
India

How to retrieve the transcript from gff file with the longest isoform? I tried using R , Biophython and other software also in ubuntu, but got so many installation problems.

If you anyone could help in solving this, it will be helpful.

The R script as follows,

 proteome <- biomartr::getProteome(db = "refseq", organism = "Arabidopsis thaliana")
annotation <- biomartr::getGFF(db = "refseq", organism = "Arabidopsis thaliana")
# retrieve longest isoforms and store in new file
retrieve_longest_isoforms(proteome_file = proteome, 
                          annotation_file = annotation, 
                          new_file = "Athaliana_pep_longest.fa")
# import new file into R session                          
Athaliana_pep_longest <- Biostrings::readAAStringSet("Athaliana_pep_longest.fa")

Thanks in advance. Aneesha

FunctionalAnnotation StatisticalMethod • 172 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 17 hours ago
United States

Neither biomartr, nor orthologr are Bioconductor packages, so this isn't the place to ask questions about them.

However, the question you ask and the code you present are different things. You appear to want the 'transcript with the longest isoform,' which could mean lots of different things, really, but the code you present is a way of getting the amino acid sequence for the longest isoform. Those are different things! It is probably pretty trivial to get what you want using Bioconductor, but you need to first define exactly what you are trying to get, and perhaps what you plan to do with it, after which we may be able to provide pointers.

0
Entering edit mode

Okay, Thank You

ADD REPLY

Login before adding your answer.

Traffic: 432 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6