From the log file it seems snippy doesn't do assembly. It only does mapping:

fastq --> SAM --> BAM --> VCF --> consensus_seq/snps

But if multiple samples are given (with snippy-multi option) it outputs an aligned fasta format. My understanding so far is it doesn't go through any MSA alignment algorithm, but use the pairwise alignment information from SAM files.

How is the aligned fasta file generated? The log file doesn't contain information about how snippy calculates it. Is there any other tool that I can use to get similar aligned fasta file from SAM/VCF file without going through any actual MSA algorithm?

  • $\begingroup$ I tried to make the question bit more clear, but I am not sure if I got right what you were asking for (maybe that's why I had a need to try to improve it), please, don't hesitate to edit the question if I the changes I made were inappropriate. $\endgroup$ Mar 28, 2019 at 10:31


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.