This is a test version of Biostars. For the public version, visit https://www.biostars.org.
tblastn extend beyond homology region

Hi,

I am performing a very simple task. I got a protein sequence, with conserved domains and I am mapping this protein and domains to a genome. I found the analogous domain. However I would like to get 100 amino acids downstream and upstream the domain. Because its a tblastn I am a bit stuck on how to do that properly to get everything on the same reading frame.

Any suggestion?

Thanks !

blast
• 1,003 views • link
updated by maxime.policarpo maxime.policarpo • written by Amaranta_Remedios • Amaranta_Remedios

1 answer

Are you trying to find the orthologous sequence from your protein in a genome ? If so, you can just extract your tblastn region with samtools faidx (like 100000bp upstream + 100000bp downstream) and then use exonerate to find the CDS of your gene on this genome, which will give you the exon sequences and also the intron borders.

• 0 views • link
written by maxime.policarpo maxime.policarpo

Log in to answer this question.

More posts like this