# How snippy makes MSA-like aligned fasta output from multiple samples?

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?

• 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. – Kamil S Jaron Mar 28 '19 at 10:31