Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,199	0	T	G	45.0%	50.3 / 10.8	40	T35P (ACC→CCC)	thiM	hydoxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (6/16);  new base G (18/0);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.48e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGAC  >  NC_000913/2185064‑2185336                                                                                                                                        |                                                                                                                                          gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTATCACCCCATCAtt                                                                                                                            >  1:230485/1‑151 (MQ=255)    cgcgTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTgg                                                                                                                          >  2:518673/1‑150 (MQ=255)      cgTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCa                                                                                                                       <  2:31209/151‑1 (MQ=255)        tGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACAACATTGGTCATg                                                                                                                     >  1:248408/1‑151 (MQ=255)        tGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGGACCACATCATTGGTCAt                                                                                                                      >  2:157961/1‑150 (MQ=255)         ggCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCAACAAATCATTGGTCATg                                                                                                                     >  2:373413/1‑150 (MQ=255)            tGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGGACCACATCATTGGGCATGCAg                                                                                                                  >  2:272615/1‑150 (MQ=255)                   gtgtGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAag                                                                                                          <  1:381761/151‑1 (MQ=255)                     gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCTCATCATTGGTCATGCAGTGCACAaga                                                                                                         <  1:585818/150‑1 (MQ=255)                     gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAagag                                                                                                        <  2:121530/151‑1 (MQ=255)                     gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGCATGCAGTGCCCAagcg                                                                                                        >  1:70213/1‑149 (MQ=255)                       gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgg                                                                                                       <  2:722124/150‑1 (MQ=255)                       gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgg                                                                                                       <  1:537466/150‑1 (MQ=255)                          aaCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGaa                                                                                                   <  1:433683/151‑1 (MQ=255)                            cGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAggggga                                                                                                   >  2:495456/1‑147 (MQ=255)                               tAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGt                                                                                                >  2:22317/1‑149 (MQ=255)                                aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGAATGTTgg                                                                                             >  2:597032/1‑151 (MQ=255)                                aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGtt                                                                                               >  1:317245/1‑149 (MQ=255)                                  tCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGGACCACCTCCTTGGTCATGCAGGGCACAAGGGGGGAAAGGTgg                                                                                             >  1:591808/1‑149 (MQ=255)                                   caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGTCATGCAGTGCCCAAGAGGGGGATGTTGGTgg                                                                                          >  2:17723/1‑150 (MQ=255)                                    aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                                                                                          <  1:541849/150‑1 (MQ=255)                                      caaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCCCATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTgggggaa                                                                                        >  1:319634/1‑145 (MQ=255)                                                gCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGGGGAAAAgggggaa                                                                               >  1:617580/1‑144 (MQ=255)                                                   aTAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGt                                                                          <  1:88822/151‑1 (MQ=255)                                                     aGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTg                                                                         <  2:133080/150‑1 (MQ=255)                                                            aaCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGAtt                                                                 <  2:567048/151‑1 (MQ=255)                                                                   tgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTg                                                           >  2:480326/1‑150 (MQ=255)                                                                          ttCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCCTGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAAccc                                                      >  1:277119/1‑148 (MQ=255)                                                                           tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcag                                                    <  2:121121/149‑1 (MQ=255)                                                                           tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGCATGCAGTGCACAAGGGGGGGATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAAacccg                                                    >  2:264356/1‑147 (MQ=255)                                                                             ggTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGAACCCAGCAgg                                                >  2:560883/1‑151 (MQ=255)                                                                                ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                              <  2:687662/150‑1 (MQ=255)                                                                                     tAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCCTGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACtt                                         >  1:632788/1‑150 (MQ=255)                                                                                           tCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAATGTTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGGGAAAAAAGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACCTTCAt                                     >  1:518056/1‑148 (MQ=255)                                                                                            cGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGtt                                 <  1:402891/151‑1 (MQ=255)                                                                                              cTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGGGAAACAGGTGTAACGCGTGCGCCGAATGCGCTGAACCCAGCAGGTCGACTTGCATAGttt                                >  1:165431/1‑150 (MQ=255)                                                                                                        aCCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCAt                     <  2:342128/151‑1 (MQ=255)                                                                                                          cGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGCACCACATCATTGGTCAGGCAGGGCACAAGAGGGGGATGTTGGTGAAAAAGGGGGAACGCGTGGGCAGAATGCGCTGAACCCAGCAGGTCGGATTGCATAGTTTGCTCCTGCCATa                    >  1:108816/1‑150 (MQ=255)                                                                                                                   gCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTgaagaa          <  2:546221/151‑1 (MQ=255)                                                                                                                           ttGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGAc  >  2:716682/1‑151 (MQ=255)                                                                                                                                        |                                                                                                                                          GAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGAC  >  NC_000913/2185064‑2185336

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