xavierdidelot / ClonalFrameML

ClonalFrameML: Efficient Inference of Recombination in Whole Bacterial Genomes
GNU General Public License v3.0
109 stars 27 forks source link

option -fasta_file_list true? #152

Closed goshng closed 6 months ago

goshng commented 6 months ago

Hi,

I want to use -fasta_file_list opton because I have three alignment files in FASTA format: l1.fa, l2.fa, and l3.fa, and want to use them as blocks in ClonalFrameML. File l.txt has a line: l1.fa l2.fa l3.fa.

ClonalFrameML tree.nwk l.txt cf_l -fasta_file_list true

ClonalFrameML v1.13
fasta_file_list = true
Finished reading in control file.

ERROR: could not find listed file l.txt

I've attached a compressed file with the files that I tested. I tried the source code version and conda package version, and either produced the same error.

$ ls
l1.fa  l2.fa  l3.fa  l.txt  run  run.log  tree.nwk

Did I misunderstand the usage of the option -fasta_file_list?

Thank you, SangChul

t.tar.gz

xavierdidelot commented 6 months ago

The fasta_file_list option is useful if you have each single sequence in a separate fasta file. But what you have here is separate alignments of genes in each fasta file. To handle this, ClonalFrameML needs you to first put all the genes in a XMFA file and then use the xmfa_file option: cat l1.fa > l.xmfa;echo = >> l.xmfa;cat l2.fa>>l.xmfa;echo = >> l.xmfa;cat l3.fa>>l.xmfa;echo = >> l.xmfa Then you should run ClonalFrameML tree.nwk l.xmfa cf_l -xmfa_file true