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.
•
written
by
rob234king
1 answer
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.
Log in to answer this question.
-
MAKER has many FAILED
written by san96 •
-
Herald: The Biostar Herald for Monday, February 12, 2024
written by Biostar
-
Tutorial: Prediction of Ribosomal RNA Genes Using RNAmmer Software
written by Arslan Tariq •
-
Efetch from eutilities gives '400 Bad Request'; how to fix this?
written by Shraddha •
-
How to obtain a novel virus' complete genome from de novo assembly
written by brianna.flynn •
-
Blast non-contiguous sequences on same protein
written by michaelsacco •
-
Forum: Erik Garrison of FreeBayes uses this
written by Istvan Albert
-
Tool: Pyteomics - An Open-Source Python Framework For Mass Spectrometry (Ms) Based Proteomics
written by lev.levitsky •
-
Tutorial: Mapping Reads With Bwa And Bowtie
written by MSU course 2011 •
-
Score Protein Variants Based On Frequency Of Aa In Multiple Sequence Alignment
written by Tim
