Skip to content

Question about HMMER output and LTRs #337

@BioFalcon

Description

@BioFalcon

Hi,

I'm trying to run nhmmer (v3.4), trying to use a model I constructed that has a structure of an LTR. Bellow you can see how the self alignment looks like

Image

And the alignment itself looks like this (each color regresents a different base, gaps are white):
5' End
Image

3' End

Image

The alignment itself is around 2.1kb. I constructed the hmm model using the command

hmmbuild --symfrac 0 --dna --informat afa --fragthresh 1 -n TEST LTR.hmm LTR.aln.fa

However, the results always falls short at 200bp (the start of the alignment). I used the following command to do the search
nhmmer --cpu 8 --tblout LTR.out --notextw --noali --seed 1992 ./LTR.hmm HMMDB_Genome

When I trim the LTR on the 3'end manually, suddenly I start getting hits that extend onto the first 200bp. Any ideas how to solve this?

Thanks

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions