Fishing nucleotide sequences from genome files given aligned protein FASTA files using biopython
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:
• 1,210 views
•
link
updated
by
GenoMax
•
written
by
mirali1025 •
•
written
by
mirali1025 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
NCBI Relevant files not being returned as query results. Trying to build exhaustive dataset.
written by chordata •
-
Protein filtering for annotation
written by Ric
-
Fasta with multiple sequences to alignment object which can by used to build a phylogenetic tree
written by mac03pat •
-
Intron extraction with Biopython.
written by Evergreen •
-
How to extract intronic regions from a Genbank file?
written by Evergreen •
-
How to use for loops to compare each DNA sequence in one file to each sequence in another file?
written by rmartson •
-
Splitting and Extracting Features in fasta format from Genbank Files using Biopython
written by swknudsen •
-
Download complete bacterial genomes and associated plasmid sequences from NCBI
written by wanderingstefan •
-
Phylogeny Study Based On Whole Genome Sequences
written by GAO Yang
-
Refseq Genbank To Fasta Format Failing With Contig Fields
written by Käpt'N Silico •
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.