Question: Fishing nucleotide sequences from genome files given aligned protein FASTA files using biopython
gravatar for mirali1025
12 months ago by
mirali10250 wrote:

Hello everyone:

I have written a script using biopython to try and fish out nucleotide sequences of aligned proteins (within FASTA files). My task is to use local genbank files of particular strains and try to extract the nucleotide sequences based on the aligned FASTA files. I'm having some trouble finishing the script, however, but I think it is almost complete:

Biopython script

fasta assembly gene • 346 views
ADD COMMENTlink modified 12 months ago by genomax70k • written 12 months ago by mirali10250

One has to request permission to view your script and, sorry, I do not want use my personal email account, nor create a phony account just to try to help someone at a forum. Why don't you use a github gist?

If you make your alignment output gff or bed, you can then use bedtools to obtain the underlying sequence.

ADD REPLYlink written 12 months ago by h.mon26k
Please log in to add an answer.


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