Dear all,
After aligning protein sequences of a gene of interest (using mafft) and generating a reference tree (using IQTREE), I wanted to align a set of environmental sequences into the obtained reference tree using PaPaRa. I used the following command line:
papara -a -t ${ref_tree} -s ${ref_alignment} -q ${query_seq} -j 12 -n report
But I obtained the following error message :
terminate called after throwing an instance of 'std::out_of_range'
what(): vector::_M_range_check
/var/spool/slurm/d/job866081/slurm_script: line 39: 36795 Aborted
Both reference tree and alignments files are in the write format and try to run this command line on my personal laptop and on the cluster of my institute but got the same error message.
Any idea of the problem here?
Thanks !!