Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Mafft alignment failed to save #24

Open
ALongmuir01 opened this issue Sep 23, 2020 · 1 comment
Open

Mafft alignment failed to save #24

ALongmuir01 opened this issue Sep 23, 2020 · 1 comment

Comments

@ALongmuir01
Copy link

ALongmuir01 commented Sep 23, 2020

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.

@Tabima
Copy link
Member

Tabima commented Sep 23, 2020

Hi,

It seems to be that the path for either MAFFT or HMMER are not correctly specified. Make sure that MAFFT is at the c:/mafft/ path, also try and remove the last backslash (i.e. c:/mafft). Same with HMMER.

In addition, if you have problems after you test the new paths and make sure the programs are running correctly (as mentioned in the README file of effectR), can you send me a more comprehensive summary of the output of the run, pleas?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants