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%	96.6 / NA	31	D666D (GAC→GAT)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

TCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAA  >  NC_000913/1954765‑1955236                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        tcacgCACGCATTCCACGCCTTTTTTGATCCCGGCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCTCCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCtcatcg                                                                                                                                                                                                                                      <  2:30155/245‑1 (MQ=255)       tcgcataataCGCCTTTTTTAATCCCGTCAGTCAGCACAGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCAACAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGAGAATAGTGATTGGTTCACGATCAGCGTACGCATATTTTTTACGCTGTTGCGCATAGTTAAGCTGTCTGTGTAAACGGTGGGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTAATCGGCGGGAcc                                                                                                                                                                                                                              >  2:48442/10‑251 (MQ=255)                 cGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATtca                                                                                                                                                                                                                    <  1:441702/251‑1 (MQ=255)                                cGTCAGTCATCACCGCGCCTGTCAAAATCTGACCACGTTGGGTCTAGCCACGCACCAGATCGCCATTCGCAATATCAAAGCGGGCTGCATCTTCGGTGTGAATAGTGATTGGTTCACGAGCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTGAAGCTGACTGTGTAAACGGGGTGCCGGATGTGCGGTCAGAAGCCGTAACTGCTTCTCATCGGCGTTACCTTTCAATTCATCACTAGccagcc‑ag                                                                                                                                                                                                    <  2:336875/251‑1 (MQ=255)                                    agtcaCCACCGCGCCTGTCAAAATCTGACCACTTTGGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCACTCATCAGGAGCCAGCC‑AGGTTg                                                                                                                                                                                                <  2:86128/251‑1 (MQ=255)                                                  cTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTCCCTTTCCATTCATCAGGagaccgcacgggtgggtgtggcgggga                                                                                                                                                                                    >  2:441745/1‑225 (MQ=255)                                                       aaaaTCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAGAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGCTCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGACGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATACGGAGCCAGCC‑ATGTTGGGGGTGCCGGG‑CAATccct                                                                                                                                                                            >  1:62244/1‑249 (MQ=255)                                                                   cGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCCTAAGGAGCCAGCC‑AGGGTGGGTGGGGCGGG‑CAAACCTTATAG‑A‑CCAAttt                                                                                                                                                              >  1:359747/1‑251 (MQ=255)                                                                        tttCCAGAAACGCCCCAAATCGCCATTCTCTATTCCAAAGCGTGCAGCATCTTCGGTGTTAATAGTGATTGGTTCACGATCGGCGATCGGATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTcc                                                                                                                                                          <  1:492074/249‑1 (MQ=255)                                                                         ttCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTTCCTTTTTCAGGAGCCAGCC‑AGGTGGGGTGTGACGGG‑GACTTCTTATTG‑C‑CCAATTTTTCCa                                                                                                                                                         >  2:266367/1‑250 (MQ=255)                                                                               acacGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCCGGAAGCAGCC‑AGGGTTGGTGTGCCGGG‑CCATCCTTTTAA‑C‑CAAAATTTTTCAACGtttt                                                                                                                                                  >  1:173283/1‑251 (MQ=255)                                                                                              ccATTCGCAATACCAAAGCTGGTAGCAGCTTCGGTGTGAATAGTGATTGGTTCACAGTCGGCGATCGCATAGTTTTTACGCAGTTCCGCATAGTTAATCTGACTGTGTAAACGGTGTGCCTGATGTGCGGTCAGAAGCTGTACCTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGTCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGAGTTCAAt                                                                                                                                   <  2:3501/251‑1 (MQ=255)                                                                                                  tCGCAATTCCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGATCACGATCGGCGTTCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCTGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCTATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTGTTCCAACGTTTTTGAATAGATTTCAAtttt                                                                                                                                <  1:552223/250‑1 (MQ=255)                                                                                                    gCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGACAGCC‑AGGGTGGGGGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTTCAACGTTTTTGAATAGATTTCAATTTTGc‑c                                                                                                                            >  1:219341/1‑251 (MQ=255)                                                                                                         aCCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTg                                                                                                                       <  1:6992/251‑1 (MQ=255)                                                                                                                    cagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCCACGTTTTTTTATAGATTTCAATTTTTC‑CGCTTTGCGTACCCAgc                                                                                                            >  1:424706/1‑251 (MQ=255)                                                                                                                                        gTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTTCTCAGGAGCCCGCC‑CGGTTGGGTGGGCCGGGGCAA‑ACTATTAG‑C‑CAAATTTTTTCACATTTTTTGGATAGATTTAAATTTTTCTCCTTTGGC‑TACCCAGCCCATTTTTccacagctcggg                                                                                        >  2:264864/1‑240 (MQ=255)                                                                                                                                                        cGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCCTTCATCCGGGGCCAGCC‑AGGGTGGGTGTGCCGGG‑CAATCCTTATAGAC‑CAAATTTT‑CCCACGTTTTTGAATAGATTTCAATTTTTC‑CGCTTGGCGTACCCAGCGCCTTTTTCACCGGCTCGGCGCGGGAATCAACATa                                                                         >  1:233830/1‑250 (MQ=255)                                                                                                                                                                  atatTTTTTACGCAGTGCCGCATAGTTATGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCGTGGCGTACCCAGCGCATTTTTAACCGGATCGGCGCGGAAATCACCATAACGAACGTACt                                                              <  1:356668/251‑1 (MQ=255)                                                                                                                                                                       ttttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTTAAG‑C‑CAAATTTTTCCAACGTTTTTGTATAGATTTCAATTTTTC‑CGCTTTGCGTACCCAGCGCATTTTTCTCCGGGTCGGCGCGGAAATCACCATAACGAACGTCCTgttct                                                         >  1:341333/1‑250 (MQ=255)                                                                                                                                                                        tttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTACATTGCTCCGAAGCCACCCCAGGT‑GGGTGTGGCGGG‑CAATACTTTAAG‑CACCAATTTTTC‑ACCGTTTTTAAGTAAATTTTAAATTTTA‑CTCTTGTGGGAACCAGCGCATTTTTTACCGGAGCGTCGCGGGAATCACAATacacgacggtattgtctt                                                        >  2:398115/1‑234 (MQ=255)                                                                                                                                                                         ttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGGGTGCCGGG‑CAATCCTTATAG‑C‑CCAATTTTTTCAACGTTTTTTTATAGATTTCAATTTTTC‑CGCTTGGCGTTCCCAGCGCATTTTTTACCTGATCGGCGCGGGAATCACCATAACGAACGTACTgttcgtt                                                       >  1:79092/1‑251 (MQ=255)                                                                                                                                                                               aGTTCCGCATATTTAAGCGGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCg                                                                                                                                                                                       <  2:304055/121‑1 (MQ=255)                                                                                                                                                                               aGTTCCGCATAGTTAAGCTGACTGTGTAACCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGAGCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTtctc                                                 <  1:329002/251‑1 (MQ=255)                                                                                                                                                                               aGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTTCATTCATCAGGAGCCAGCC‑AGGTTGGGGGTGCCg                                                                                                                                                                                       >  1:304055/1‑121 (MQ=255)                                                                                                                                                                                                      gtgtAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTGCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCTCCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGtt                          <  1:48442/251‑1 (MQ=255)                                                                                                                                                                                                                     cGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGa           <  1:434788/251‑1 (MQ=255)                                                                                                                                                                                                                          gtgCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCa      <  1:441745/251‑1 (MQ=255)                                                                                                                                                                                                                            gCGGTCAGAAGCTGTAACTGCCTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGACCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCAtt    <  2:549622/251‑1 (MQ=255)                                                                                                                                                                                                                             cGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCCGCC‑AGGGTGGGGGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTTC‑CGCTTGGCGTACCCAGCGCATTTTTTACCGGATCGGCGCGGGAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGGTCAGTTTATTTTTCTTCCAGAAGGGAtt    >  1:24325/1‑250 (MQ=255)                                                                                                                                                                                                                              ggTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCCTTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTTTAT‑C‑CAAATTTTTTCAACGTTTTTTAATAGATTTCAATTTTTC‑CCCTTTGCGTACCCCGCGCCTTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCCCTGTTGCCCATTTCGATCCGTTTTTTTTTCTTCCAGAAGGccttt   >  1:30930/1‑246 (MQ=255)                                                                                                                                                                                                                               gTCAGAAGCTGTAACTGCTTCTCATCGGCGGAACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGGGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTTC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAAAGAAAGTACTGTTCGTTCTTCTCGCAGTGGCGCATTTCGGTCAGTTTATTTTGCTGCCAGAATGCATTaa  >  1:405946/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        TCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCC‑AGGTTGGGTGTGCCGGG‑CAATCCTTATAG‑C‑CAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGC‑CGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAA  >  NC_000913/1954765‑1955236

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