This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Where can I download HMM profile for different species?

I would like to ask is there any database allow me to download HMM for different species such as: 1. Pan troglodytes 2. Pongo abelii 3. Rhinopithecus rexellana

hmm profile species genome sequence
• 4,357 views • link
written by allyson1115ar • allyson1115ar

Which HMM do you need?

• 1 views • link
written by dago dago

Yes I understand, but do you need the HMM for the nucleotides seqs of all genes, for the proteins? Also, I am not sure you have hmm specific for these species, maybe you can try to use Pfam domains and start with that for gene annotation.

• 1 views • link
written by dago dago

I can't find any hmm link to download the above 3 species :/

• 1 views • link
written by allyson1115ar • allyson1115ar

Which kind of HMM profile are you looking for ? => I mean is it for coding genes ? For what purpose do you need it ? More precisely, which tool do you plan to use with ?

• 1 views • link
written by Juke34 Juke34

i would like to include the hmm in snap ab initio gene predictor since the package doesn't provide the other species

• 1 views • link
written by allyson1115ar • allyson1115ar

Each ab-initio tool has its own hmm profile format. And they don't necessarily use the same kind of information. So, I doubt you can use the profile coming from another tool as HMMER as suggested. I guess you are trying to annotate the gene coding. So the simplest is to create the profile yourself. Download the gff and the genome of your species (on Ensembl as example). Put the fasta sequence within the gff file as it should be (Introduced by the ##FASTA line at the end of the file). Then prepare the needed files to train snap by using the maker2zff.pl tool that is coming with the MAKER genome annotation tool package. You will get two output files:genome.ann and genome.dna Then to create the profile you will have to launch those commands:

fathom -categorize 1000 genome.ann genome.dna
    fathom -export 1000 -plus uni.ann uni.dna
    forge export.ann export.dna
    hmm-assembler.pl SPECIES . > SPECIES.hmm
• 1 views • link
written by Juke34 Juke34

It sounds like a reasonable approach!

• 1 views • link
written by dago dago

0 answers

No answers yet.

Log in to answer this question.

More posts like this