Question: How to extract 3', 5' UTR sequences from genbank records using python, PERL and R code?
0
gravatar for mathavanbioinfo
11 months ago by
India
mathavanbioinfo50 wrote:

Hello All, I have 1000 sequences of genebank records, I want to extract the only 3'UTR, 5'UTR sequences from the sequences and to store in excel format. Share your ideas and suggestion [using PERL or Python or R codes]

utr • 488 views
ADD COMMENTlink modified 11 months ago by zubenel110 • written 11 months ago by mathavanbioinfo50

Hi, please post a sample gbk file and define the headers that you want to see in your output file (Ex: seqID, locusTag, sequence ... ).

ADD REPLYlink written 11 months ago by hugo.avila160

Please take a look at the biopython cookbook and tutorial.

ADD REPLYlink written 11 months ago by WouterDeCoster44k
0
gravatar for padwalmk
11 months ago by
padwalmk100
padwalmk100 wrote:

Hi, It's unclear wither you have the gff file with you or fasta.

You can look in to following post

Extract coordinates of upstream region up to closest coding region in R

ADD COMMENTlink modified 11 months ago • written 11 months ago by padwalmk100
0
gravatar for zubenel
11 months ago by
zubenel110
zubenel110 wrote:

If you have gff file you might try to use gff2fasta.pl with option -feature set as "five_prime_UTR" or "three_prime_UTR" or something like that. Also you may read how to get sequences of specific features with BioPerl.

ADD COMMENTlink modified 11 months ago • written 11 months ago by zubenel110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1119 users visited in the last hour