This is a test version of Biostars. For the public version, visit https://www.biostars.org.
find position of 3' utr at entrez

Hi

I am using the BIOPYTHON library Bio to retrieve sequences from Entrez, I need to know the ending position of 5' UTR and the starting position of 3' UTR. How do I get them?

thank you

p

sequence entrez biopython
• 1,688 views • link
written by paumarc • paumarc

Are you fetching the records in the Genbank format? If yes then you can parse the features to identify the UTR tags.

• 0 views • link
written by Sej Modha Sej Modha

Can you tell me any example of how to get or tutorial them? it is driving me mad

thanks

• 0 views • link
written by paumarc • paumarc

Following code can be modified to check for UTR instead of CDS for example:

from Bio import SeqIO
records = SeqIO.parse("sequence.gb", "genbank")
for seq_record in records:
    for feature in seq_record.features:
        if 'CDS' in (feature.type):
            if 'methyltransferase' in str(feature.qualifiers["product"]):
                print(feature.qualifiers["locus_tag"],feature.qualifiers['product'])

See more on the BioPython webpages: http://biopython.org/DIST/docs/tutorial/Tutorial.html#htoc130

• 0 views • link
written by Sej Modha Sej Modha

0 answers

No answers yet.

Log in to answer this question.

More posts like this