Hi,

I have been trying to replicate the test analyses in the discosnp manual, but I am running into trouble with the extending the central sequence with contigs and unitigs.

Specifically, using v2.0.0 I get identical output (i.e. a .fa file with central sequences only, pasted below) when I run these three commands:

./run_discoSnp++.sh -r "data_sample/reads_sequence1.fasta.gz data_sample/reads_sequence2.fasta.gz" -p test
./run_discoSnp++.sh -r "data_sample/reads_sequence1.fasta.gz data_sample/reads_sequence2.fasta.gz" -p test -T
./run_discoSnp++.sh -r "data_sample/reads_sequence1.fasta.gz data_sample/reads_sequence2.fasta.gz" -p test -t

I have checked that it's nothing to do with passing the arguments through the script: it makes no difference to the output if I change the default in the script instead of passing the argument via the command line, i.e. changing

extend=""to extend="-T" or extend="-t" and running the script without passing additional options still produces the same .fa file as pasted below. (Note, I delete the file after each run, just to make sure..., and the script finishes successfully each time).

Any tips? Am I doing something wrong?

Thanks,

Rob

Output from all commands in this file:

test_k_31_c_4_D_0_P_1_b_0_coherent.fa

>SNP_higher_path_3|high|nb_pol_1|C1_0|C2_134|rank_1.00000
CAAGTGCACTTCCACAGAGCGCGGTAGAGAGTCATCCACCCGGCAGCTCTGTAATAGGGAC
>SNP_lower_path_3|high|nb_pol_1|C1_124|C2_0|rank_1.00000
CAAGTGCACTTCCACAGAGCGCGGTAGAGACTCATCCACCCGGCAGCTCTGTAATAGGGAC
>SNP_higher_path_2|high|nb_pol_1|C1_0|C2_86|rank_1.00000
ACTAATAGGGCCGGGCTACATGTTAACTACTAGGCTATAACCTATTGATGACCCGGTCCAT
>SNP_lower_path_2|high|nb_pol_1|C1_74|C2_0|rank_1.00000
ACTAATAGGGCCGGGCTACATGTTAACTACAAGGCTATAACCTATTGATGACCCGGTCCAT
>SNP_higher_path_1|high|nb_pol_1|C1_110|C2_0|rank_1.00000
GCAGCGCAACAACGCAACAGCTCGAGGTGTTCTTCGCAGAGAAACCGCACGTCCAGTTCTA
>SNP_lower_path_1|high|nb_pol_1|C1_0|C2_114|rank_1.00000
GCAGCGCAACAACGCAACAGCTCGAGGTGTACTTCGCAGAGAAACCGCACGTCCAGTTCTA