Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,789	0	C	A	40.4%	106.9 / 18.6	47	V24F (GTT→TTT)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base C (3/25);  new base A (19/0);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 2.21e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAT  >  NC_000913/1535640‑1535937                                                                                                                                                        |                                                                                                                                                      aCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAaaca                                                                                                                                                      <  1:127600/151‑1 (MQ=255)      aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAaacaa                                                                                                                                                     >  1:766103/1‑147 (MQ=255)          cGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCa                                                                                                                                             <  2:746757/151‑1 (MQ=255)              aCTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGaa                                                                                                                                          >  2:682168/1‑150 (MQ=255)                ttCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGGAACAAAAcccgaccaat                                                                                                                                         >  2:791968/1‑141 (MQ=255)                      cGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGa                                                                                                                                 <  1:595678/151‑1 (MQ=255)                          ggCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAAc                                                                                                                              <  2:552446/150‑1 (MQ=255)                               aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt                                                                                                                        <  1:595556/151‑1 (MQ=255)                               aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtg                                                                                                                         <  1:209172/150‑1 (MQ=255)                                   tAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCCGGAATTG‑CGCGGAAAGGGTGTgggg                                                                                                                    >  2:509290/1‑151 (MQ=255)                                    aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAgg                                                                                                                    <  2:588262/150‑1 (MQ=255)                                      cGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCAGAAATGCCGCGGAAAGGGTGTgggg                                                                                                                    >  2:494206/1‑148 (MQ=255)                                       ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAACACCCGCCGAAATGCCGCGGAAAGGGTGTgggg                                                                                                                    >  2:251780/1‑147 (MQ=255)                                         aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAg                                                                                                              <  1:122763/151‑1 (MQ=255)                                          tGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCAGAATTGCCGCGAAAAGGCTGTGGGGGaaaa                                                                                                               >  2:331973/1‑149 (MQ=255)                                           gcgcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGc                                                                                                             <  1:690547/150‑1 (MQ=255)                                           gcgcACATCGCTGTCGTCGCTGAGCAGCAGAAAGGCAGTTCCGGAATGG‑CTAACGTCCGGGGCGATAAGCGGCATTGGCGTTTCCTTAGGATGAAGGGC‑GGGGGAGAAAAAAACGGCGGAATTATCCCGaaaatgctgtggggg                                                                                                                   >  1:756065/1‑132 (MQ=255)                                                 aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAACCGGCAGAAATGCCGCGGAACGGGTGTGAGGTAAAAGCATCg                                                                                                         >  2:962825/1‑148 (MQ=255)                                                   cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGt                                                                                                     <  1:629795/150‑1 (MQ=255)                                                    gCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTgg                                                                                                   <  2:90773/151‑1 (MQ=255)                                                        tcgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCAGAAATGCCGCGAAAAGGTTGTGAGGGAAAAGAATCGACGTGGGac                                                                                                >  2:572059/1‑150 (MQ=255)                                                            cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTg                                                                                           <  1:85308/151‑1 (MQ=255)                                                             gCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGc                                                                                          <  2:319078/151‑1 (MQ=255)                                                                    agcagAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCa                                                                                   <  2:281881/151‑1 (MQ=255)                                                                          aaaGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCGGAAAAGTCGCGGAACGGGTGTGGGGGGAAAACAAAGACGCGGTGCACCTGCggggtctttta                                                                               >  2:199536/1‑140 (MQ=255)                                                                             gCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAACCCGCAGAATTTCCGCGAAACGGTTGTGGGGGAAAAAAATCGACGTGGGACACCCGCGGgtttttttaaattt                                                                          >  2:298406/1‑140 (MQ=255)                                                                                  ttCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTTCCGCGGAAAGGGTGTGGGGGAAAAACATCGACGGGGGACACCCGCGGGTTTATTAAAGTTTTcc                                                                       >  2:4574/1‑149 (MQ=255)                                                                                     cGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAAATGCCGCGAAAAGGGTGTGGGGGAAAAACAACGACGGGGGACACCCGCGGGTTTATTAAAGtcctcc                                                                       >  2:73799/1‑146 (MQ=255)                                                                                        aTGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGGCGGAA‑TGTCGCGGAACGGGTGTGGGGGAAAAAAAACGACCGGGGACACCCGCGGGTTTATTTAAATTTTCCCGTGAc                                                                 >  1:968746/1‑149 (MQ=255)                                                                                               aaC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAATTGCCGCGGAAAGGGTGGGGGGGAAAAACAACGACGTGGGACACCCGCGGGTTCATTAACGTTCTCCTGTGAATGgag                                                            >  1:135471/1‑147 (MQ=255)                                                                                               aaC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACt                                                        <  2:522362/151‑1 (MQ=255)                                                                                                    cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAt                                                      <  2:484155/149‑1 (MQ=255)                                                                                                    cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAt                                                      <  2:341655/149‑1 (MQ=255)                                                                                                    cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑Ca                                                   <  2:830826/151‑1 (MQ=255)                                                                                                     gcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CAt                                                  <  1:251775/151‑1 (MQ=255)                                                                                                       gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAAATGCCGCGCAAAGGTTGTGGGGGAAAAAAATCCACGCGGGACACCCCCGGGTTCATTTAAGTTCTCCCGTGGCTGGAGGACTa                                                       >  2:22947/1‑145 (MQ=255)                                                                                                        ggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGc                                               <  2:904036/151‑1 (MQ=255)                                                                                                            aTAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGGAATGTCGCGGAAAGGTTGTGGGGGAAAAAAATCGACGTGGGACACCCGCGGGTTTATTAAAGATCTCCCGTGACTGGAGAAATAT‑CATaaaccg                                            >  1:883715/1‑145 (MQ=255)                                                                                                              aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCa                                          <  2:553515/150‑1 (MQ=255)                                                                                                                   gCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTg                                      <  2:459032/149‑1 (MQ=255)                                                                                                                     aTTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGGAATGCCGCGGAAAGGGTGTGGGGGGAAAAAAACGACGGGGGGCAACCGCGGGTTCATTAAAGTTTTCCTTTGGCTGGAGGACTAT‑AATAACCTGCacaggggc                                   >  1:295489/1‑143 (MQ=255)                                                                                                                       ttGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCgg                                 <  1:979668/150‑1 (MQ=255)                                                                                                                        tGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCCGAATTGCCGCGCAAAGGGTTTGGGGGAAAAGAATCGACGCGGGGCAACTCCGGGTTTATTTAAGTTTTTCTGTGGCCGGGGAAAAAT‑CATAACCCGCCAGGGGGCgg                                 >  1:711245/1‑149 (MQ=255)                                                                                                                        tGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAACAACGGCAGAAATGGCGCGGAAAGGGTGTGAGGGAAAAAAATCCACGGGGGACAACCGCGGGTTCTTTAAAGTTTTCCCTTGGCCGGAGAAAAATCCATAACC‑GCAAGTGGGCggg                                >  2:384044/1‑149 (MQ=255)                                                                                                                           gcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGc                            <  2:228124/151‑1 (MQ=255)                                                                                                                                    aGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                   <  2:909548/151‑1 (MQ=255)                                                                                                                                                  ggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg      >  1:648118/1‑151 (MQ=255)                                                                                                                                                    aaaaCAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGa     >  2:378228/1‑150 (MQ=255)                                                                                                                                                     aaaCAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaaa   <  1:437527/151‑1 (MQ=255)                                                                                                                                                      aacaacACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt  >  2:505092/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      ACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAAC‑TCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAT‑CATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAT  >  NC_000913/1535640‑1535937

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