Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,955,000	A→G	L670L (TTA→CTA)	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,000	0	A	G	100.0%	109.1 / NA	32	L670L (TTA→CTA)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/12);  total (20/12)

TACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCG  >  NC_000913/1954868‑1955131                                                                                                                                     |                                                                                                                                    tACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTtc                                                                                                                              >  1:245651/1‑140 (MQ=255)    cAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCtcg                                                                                                                           <  2:89413/140‑1 (MQ=255)        gCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTcgg                                                                                                                       <  2:758379/140‑1 (MQ=255)           tGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTcggcgg                                                                                                                    <  2:245651/140‑1 (MQ=255)                caTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCt                                                                                                               >  1:520219/1‑140 (MQ=255)                           tgAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATtcatc                                                                                                    >  2:408533/1‑140 (MQ=255)                             aaTAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAg                                                                                                  >  1:476193/1‑140 (MQ=255)                               tAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGa                                                                                                <  2:476193/140‑1 (MQ=255)                                     ttGGTTCACGATCGGCGCTCGCTTGTTTTTTACGCAGTTTCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTACGGTCAGACGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGagccagc                                                                                          >  2:146947/1‑140 (MQ=255)                                       ggTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAgccagcca                                                                                        >  1:806124/1‑140 (MQ=255)                                          tCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGt                                                                                     <  2:520219/140‑1 (MQ=255)                                                  ggCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGtgtg                                                                             >  1:357097/1‑140 (MQ=255)                                                          cATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCaa                                                                     <  2:357097/140‑1 (MQ=255)                                                                   tACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAg                                                            >  1:398817/1‑140 (MQ=255)                                                                               cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTcc                                                >  1:958901/1‑140 (MQ=255)                                                                                    ttAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAATGt                                           >  2:30257/1‑140 (MQ=255)                                                                                    ttAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGt                                           >  1:694686/1‑140 (MQ=255)                                                                                    ttAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGt                                           >  2:237697/1‑140 (MQ=255)                                                                                       aGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGtttt                                        >  2:10449/1‑140 (MQ=255)                                                                                              tgtgTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAg                                 >  2:776751/1‑140 (MQ=255)                                                                                                          gtgCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGGATAGATTTCAATTTTg                     >  2:896198/1‑140 (MQ=255)                                                                                                          gtgCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAtttt                      <  1:284259/139‑1 (MQ=255)                                                                                                          gtgCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAAtttt                      >  2:284259/1‑139 (MQ=255)                                                                                                          gtgCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTg                     >  1:778323/1‑140 (MQ=255)                                                                                                          gtgCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTg                     >  1:674963/1‑140 (MQ=255)                                                                                                          gtgCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTg                     >  1:888787/1‑140 (MQ=255)                                                                                                            gCCGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGcc                   >  1:531581/1‑140 (MQ=255)                                                                                                              cGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGc                 <  2:398817/140‑1 (MQ=255)                                                                                                              cGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGc                 <  2:778323/140‑1 (MQ=255)                                                                                                              cGGATGTGCGGTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGc                 <  2:958901/140‑1 (MQ=255)                                                                                                                        gTCAGAAGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTAcc       <  2:806124/140‑1 (MQ=255)                                                                                                                              aGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAgcg  <  1:1001414/139‑1 (MQ=255)                                                                                                                              aGCTGTAGCTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAgcg  >  2:1001414/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                    TACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCG  >  NC_000913/1954868‑1955131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑