Clustalo failing to complete multiple sequence alignments

Clustalo failing to complete multiple sequence alignments

0

Hi,

I have been using clustalo in the terminal to create multiple sequence alignments for the past few weeks, but every time I try lately it fails. I’ve tried on two computers, with different files. I’ve tried looking up the errors but haven’t found anything. Does anyone know what this means, or what might be going on?

Thanks in advance!

Here’s my command:

clustalo --in=fl_L1PA2.fa --out=L1PA2.aln --force --outfmt=clustal --wrap=175 --threads=6 --verbose

And here’s the output. It goes along fine through the first few parts, and then..

Using 6 threads
Read 914 sequences (type: DNA) from fl_L1PA2.fa
Using 96 seeds (chosen with constant stride from length sorted seqs) for mBed (from a total of 914 sequences)
Calculating pairwise ktuple-distances...
Ktuple-distance calculation progress done. CPU time: 1035.86u 2.24s 00:17:18.09 Elapsed: 00:03:29
mBed created 31 cluster/s (with a minimum of 1 and a soft maximum of 100 sequences each)
Distance calculation within sub-clusters done. CPU time: 387.80u 0.50s 00:06:28.30 Elapsed: 00:01:18
Guide-tree computation (mBed) done.
HHalignWrapper:hhalign_wrapper.c:1419: problem in alignment (profile sizes: 1 + 1) (chr1_142876163-142882314 + chr1_223568877-223574899), forcing Viterbi
        hh-error-code=4 (mac-ram=8000)
hhalign:hhalign.cpp:961: Problem Reading/Preparing profiles (len(q)=0/len(t)=0)
HHalignWrapper:hhalign_wrapper.c:1447: problem in alignment, Viterbi did not work
        hh-error-code=4 (mac-ram=64000)
hhalign:hhalign.cpp:961: Problem Reading/Preparing profiles (len(q)=0/len(t)=0)
FATAL: could not perform alignment -- bailing out


multiple


sequence


alignment


clustalo


msa

• 102 views

Read more here: Source link