Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,955,010	G→A	D666D (GAC→GAT)	torZ ←	trimethylamine N‑oxide reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,955,010	0	G	A	100.0%	91.7 / NA	29	D666D (GAC→GAT)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/14);  total (15/14)

CATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAG  >  NC_000913/1954766‑1955226                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              cATGCACGCATACAAAGCCTTTATTGATCCCGTCCGTCACCACCGAGCCTGTCCAAATCTGCCCACGTTTGTTCCAGCCACGTACCCGATCGCCCTTCGCTATACCAAAGCGTGCAGCATCTTCGGAGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGAGTCAACGGTGTGCCGGGTGTGCGGTCAGAAGCTGTAACTGCTTCTCATcggcg                                                                                                                                                                                                                         >  2:262651/1‑251 (MQ=255)    gCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGtg                                                                                                                                                                                                                      >  2:351098/1‑250 (MQ=255)          ctACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGCTCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACCGCGTGCCGGATGTGCGATCAGAAGCTGTAACTGCTTCTCATCGGGTGTACC‑tctc                                                                                                                                                                                                               >  2:209272/2‑248 (MQ=255)                       tttGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGa                                                                                                                                                                                                   <  2:281801/250‑1 (MQ=255)                          gATCCCGTCAGTGACCACCGAGCGTGGCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGGTCGCCATTCGCAATACCAAAGCGTGAAGAATCTGCGGTGTAAATAGTTTTTGGCTCACGATCGGCGAACGCATATATTTTACGAAGATCAGCGTAGTTAAGCTGACTGTGTAAACAGTGTGCCGGATGGGCTGTCAGAAGCTGTAACTGCTACTAATCGGCTGTAAC‑TTTcagttacacagtaaccc                                                                                                                                                                                               >  2:316470/1‑235 (MQ=255)                                gtcagtcaCCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAgg                                                                                                                                                                                       >  1:526788/1‑251 (MQ=255)                                                             gACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAg                                                                                                                                                             <  1:444699/248‑1 (MQ=255)                                                             gACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAg                                                                                                                                                             <  1:532572/248‑1 (MQ=255)                                                                    tttGTTCCAGACACGCACCAGATCGCCATTCGAAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAGGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTTCCTTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑Attttt                                                                                                                                                  >  1:146807/1‑251 (MQ=255)                                                                        ttCCAGACACGCAACAGATCGCCATTCGCAATACCAAAGAGTGCAGCATCTTCGGTGTGAATAGAGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTAAC‑CTTCCCTTTTTCAGGAGACCA‑C‑CAGGGTGGGGGTGCCGGGCACCACTTATAACCAA‑ATTTTT‑CCa                                                                                                                                              >  2:579873/1‑250 (MQ=255)                                                                                                                    agcaTCGTCGGTGTGAATAGGGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCATTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGAGGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCATC‑CAGGTTGGGTGTGCCGGGCAATCCGTATAGCCAA‑ATTTTT‑CCAACGTTGTTGAATAGATTTCAATTTTGCCGCTTGGCTTACCCAgcg                                                                                                 <  2:239133/251‑1 (MQ=255)                                                                                                                           tCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTGTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAttt                                                                                            <  1:369745/249‑1 (MQ=255)                                                                                                                            cGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAttttt                                                                                          <  1:351098/250‑1 (MQ=255)                                                                                                                             ggTGTGAATAGTGATTGGTTCACGCTCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCCTTTCC‑TTCCTCGGGGG‑GCAGCACAG‑CTGGGGGGGACGGGCAAACCCTATAGACAA‑ATTTTTTCCAAATTTTTTGAATA‑AATTCAATTTTGCCGCTTTGCGTATCCAGCGCATTTTTc                                                                                         >  2:429113/1‑250 (MQ=255)                                                                                                                                 tgcaTAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCCTTCATCAGGCG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAAAGTTTTTGAATAGATTTCAATTTTTCCGCTTGGCGTAGCCCGCGGATTTTTTACCgg                                                                                    >  1:461540/4‑251 (MQ=255)                                                                                                                                 tgAATAGTGATTGGTTCACGATCGGAGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTACATTCATCCGGAG‑CCAGC‑CAGGGTGGGGGTGCCGGGCGATCCTTATAGCCAA‑AATTTT‑CCAACGTTTTTGAATAGATTTCATTTTTGCCGCTTTGCATACCCAGCGCCTTTTTCACCgg                                                                                    >  1:480670/1‑251 (MQ=255)                                                                                                                                              gTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTACCTGCTTCTAATCGGCGGTACA‑CTTCCATTCATCAGGAG‑GCAGC‑CACGCTGGGTGTGCCGGGCAATCCGTATA‑CCAACATTTTT‑CCATCGTTTTTGAATAGAATTAAATTTTTCCGCTTTGCTTACACAGCGCATTTTTCCCCGGATTGGCGCCGaa                                                                        >  1:320164/1‑250 (MQ=255)                                                                                                                                                       cGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATa                                                               <  1:130670/250‑1 (MQ=255)                                                                                                                                                         gCGATCGCATATTTTTTACGCAGTTCCGCCTAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTAAC‑TTTCCATTCATCAGGAG‑CCAAC‑CAGGTTGGGTGTGCCGGGCCATCCTTATAGCCAA‑ATTTTT‑CCCACGTTTTTGAAAAGATTTTAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATTGGCGCGGGAATCCACATaacg                                                            >  2:390347/1‑251 (MQ=255)                                                                                                                                                           gATCGCATATTTTTTACGCAGTTCCGCATAGGTAAGCTGACTGTGTAAACGGTGTGCCGGATGTTCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTGCC‑CTTACCTTCATCAGGAG‑CCAAC‑CAGGGTGGGTGTGCCGGGCAAGCCTTCTAACCAA‑ATTTTT‑CCAACGCTTTTTAGTTGGATTCAATTTTGCCTCCGTGCGGACCCCGCGCATTTTTCCCCGGGTCGGCGCGGAAATCACCATaaac                                                            >  2:548093/1‑247 (MQ=255)                                                                                                                                                                       tttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTgttcg                                               <  1:570295/250‑1 (MQ=255)                                                                                                                                                                       tttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTG‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCTCCATAACGAACGTACTgttcg                                               <  1:402754/250‑1 (MQ=255)                                                                                                                                                                                    gCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGTCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCtgtg                                 <  1:531104/251‑1 (MQ=255)                                                                                                                                                                                       tAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCTGTCAGAAGCTGTAACTGCTTCGCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATGTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCGTGGCGTACCCAGTGCATTTTTCACCGGGTCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGgc                               <  1:583092/250‑1 (MQ=255)                                                                                                                                                                                             aGCTGAATGGGTAAACGGTGTGCCTGATTGGGGGTCATAAGCTGTATCGTCTTCGCATCGGCGGGACC‑TGTCCACTCAGCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCATTTTTTCTGAATAGATTTCAATGTTGCCGCGTGGCGTGCCCAGCGCATTTTTCACCGGATCGGCGCTGAAATCTCCATAACGAACGTACTTTTCGTTCTTCGCGCTTTGGCGCATTTc                        <  2:575092/251‑1 (MQ=255)                                                                                                                                                                                                         aaaCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCCTTTCTCAGGAG‑CCAGC‑CAGGTTGGGGGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAAAGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTTACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCCTTTTGATCAGTTTAtt            >  1:315785/1‑251 (MQ=255)                                                                                                                                                                                                         aaaCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCACGAG‑ACAGC‑CAGGTTGGGTGTGCCCGGCAATTCTTATAGCCAA‑ATTTTT‑ACAACGTTTTTGCATCAATTTTACATTTGCCGCTTTGCCGACCCAGCCCATTTTTCCCCGGATCCGCGCGGAAATCACCATAACGAAAGGACTGTTCGTTCTTCTTGCTGTggaggcttttggacagtttatt            >  1:562140/1‑231 (MQ=255)                                                                                                                                                                                                          aaCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGATTAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGGTCTTCCCTCGGTGGCGCATTTCGATCAGTTTAtt            <  2:113578/250‑1 (MQ=255)                                                                                                                                                                                                                   ccGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAg  <  2:127438/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              CATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACC‑TTTCCATTCATCAGGAG‑CCAGC‑CAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAA‑ATTTTT‑CCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAG  >  NC_000913/1954766‑1955226

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑