Question: search for class of proteins in unannotated genome
0
gravatar for rob234king
3.6 years ago by
rob234king570
UK/Harpenden/Rothamsted Research
rob234king570 wrote:

I would like to search an unannotated genome for all amino acid transporters (pfam:PF01490). Whats the most effective way of doing this, blast sounds easy but the proteins are multi exon, I was thinking using hmm's but not sure how and if best option given searching genome. Can anyone recommend and maybe example linux commands etc.

blast • 787 views
ADD COMMENTlink modified 3.6 years ago by Asaf6.4k • written 3.6 years ago by rob234king570
0
gravatar for Asaf
3.6 years ago by
Asaf6.4k
Israel
Asaf6.4k wrote:

You can use tblastn and search the proteins in your database against the genome. If you want to scan peptides using HMM you can use interproscan.sh from http://www.ebi.ac.uk/interpro/download.html You can configure it to search only protein families from Pfam to narrow down the search. You'll have to find ORFs before you do this or predict genes (which is harder), ORFs can be found with EMBOSS getorf.

ADD COMMENTlink written 3.6 years ago by Asaf6.4k
Please log in to add an answer.

Help
Access

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