Skip to content

Mafft alignment failed to save #24

@ALongmuir01

Description

@ALongmuir01

I'm trying to run effectR to search for RXLR effectors. When I run the hmm.search function it runs the mafft alignment but then fails to run hmmer because there's not alignment file found. I have tried this again with the test infestans files included and it's still not working.
This is my work flow:

> fasta.file <- system.file("extdata", "test_infestans.fasta", package = "effectR")
> ORF <- seqinr::read.fasta(fasta.file)
> REGEX <- regex.search(ORF, motif="RxLR")
> candidate.rxlr <- hmm.search(original.seq = fasta.file, regex.seq = REGEX , alignment.file=NULL, save.alignment=T, mafft.path = "c:/mafft/", hmm.path = "c:/Users/alongmuir/Documents/hmmer3.0_windows/")

The alignment is then printed to the console, no candidate.rxlr file is made and then I get the following message:

MAFFT alignment finished!
Starting HMM
---
Creating HMM profile

Error in hmm.search(original.seq = fasta.file, regex.seq = REGEX, alignment.file = NULL,  : 
  No MAFFT alignment found

I'm not sure if I'm doing something wrong here. Any help would be appreciated.

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