From 78ebe881b930e613b579673bd4fc071eeafdf2cd Mon Sep 17 00:00:00 2001 From: Niema Moshiri Date: Wed, 7 Feb 2024 14:51:32 -0800 Subject: [PATCH] Update README.md --- README.md | 3 +++ 1 file changed, 3 insertions(+) diff --git a/README.md b/README.md index c26ce81..ac2a448 100644 --- a/README.md +++ b/README.md @@ -61,6 +61,9 @@ Users will likely want to keep the output of the read mapper (e.g. as a compress minimap2 -t 1 -a -x sr reference.fas reads.fastq.gz | tee >(viral_consensus -i - -r reference.fas -o consensus.fas) | samtools view -b -@ 1 > reads.bam ``` +# Output Files +ViralConsensus produces 3 output files: the consensus sequence (FASTA), base counts at every position (TSV), and insertion counts before every position (JSON). You can visualize the position base counts TSV file in [SamBamViz](http://niema.net/SamBamViz) to look at coverage across the genome as well as to look at base counts at each position. + # Citing ViralConsensus If you use ViralConsensus in your work, please cite: