2
$\begingroup$

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?

Improve this question
$\endgroup$
1
  • $\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$
    – Kamil S Jaron
    Mar 28, 2019 at 10:31

0

Reset to default

Your Answer

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