This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Genome sequence from protein

I've done tblast of a certain protein and got a list of hits. Is there a way to programmatically extract genomic regions for my hits?

I'm using biopython mostly

sequence alignment genome
• 1,450 views • link
written by tim.ivanov.92 4

Are you interested in extracting a specific region of the subjects or entire fasta sequence? What sort of databases are you using? Could you post your current code?

• 1 views • link
written by Sej Modha 528

I am interested in a specific region, the database is ncbi nr

here is is the code

from Bio.Blast.Applications import NcbiblastxCommandline
tblastn='/uge_mnt/home/tim_ivanov/soft/blast/ncbi-blast-2.7.1+/bin/tblastn'
nt_db='/uge_mnt/home/tim_ivanov/DATA/fasta/nt'
blastx_cline = NcbiblastxCommandline(cmd=tblastn, query="mrp_NM_001316390.faa", db=nt_db, outfmt=5, out="opuntia.xml")
print blastx_cline
• 1 views • link
written by tim.ivanov.92 4

0 answers

No answers yet.

Log in to answer this question.

More posts like this