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

ACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTG  >  NC_000913/2185053‑2185342                                                                                                                                                   |                                                                                                                                                aCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGgtt                                                                                                                                               >  2:256488/1‑149 (MQ=255) aCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGc                                                                                                                                             <  2:323657/151‑1 (MQ=255)           tGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCAcatc                                                                                                                                     >  2:304350/1‑149 (MQ=255)             aGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCAtt                                                                                                                                  >  2:47654/1‑150 (MQ=255)                 cgTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCAATgg                                                                                                                                >  1:803649/1‑148 (MQ=255)                     gCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCa                                                                                                                         <  2:385858/151‑1 (MQ=255)                            cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGcaca                                                                                                                  <  2:571384/151‑1 (MQ=255)                            cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGcac                                                                                                                   >  2:614349/1‑150 (MQ=255)                                gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGcaca                                                                                                                  >  2:599208/1‑147 (MQ=255)                                     aaCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGaa                                                                                                         <  1:139977/151‑1 (MQ=255)                                      aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAt                                                                                                        <  1:554573/151‑1 (MQ=255)                                      aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAt                                                                                                        <  1:321986/151‑1 (MQ=255)                                          tAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGtt                                                                                                     <  2:742692/150‑1 (MQ=255)                                          tAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGTTTTCAACACACCATTGGTCATGCAGGGCACAAGAGGGGGATggt                                                                                                     >  2:657244/1‑148 (MQ=255)                                            aTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGt                                                                                                  <  2:623329/151‑1 (MQ=255)                                              caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCCGTGCACAAGAGGGGGATGTTGGTg                                                                                                 >  2:651049/1‑150 (MQ=255)                                                      cACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtg                                                                                        <  1:187168/151‑1 (MQ=255)                                                                        actgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTCAAAGTTTTCACCCCATCATTGGGCATTCAGTGCAAAAGAGGGGGATGTTGGGGGAAAAGGGGGAAAGCGGGCGCAGATTg                                                                      >  1:573061/1‑151 (MQ=255)                                                                        actgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTg                                                                      <  1:690539/151‑1 (MQ=255)                                                                          tgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGAtt                                                                       <  1:479620/148‑1 (MQ=255)                                                                               ggCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGGCATGCAGGGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAAAGCGGGCGCAGATTgcg                                                                    >  1:608056/1‑146 (MQ=255)                                                                                       cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGGGCGCAGATTGCGCTGAACCcagc                                                         >  2:870657/1‑149 (MQ=255)                                                                                       cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagcag                                                       <  2:64955/151‑1 (MQ=255)                                                                                        ggTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGTGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAgg                                                      >  2:507617/1‑151 (MQ=255)                                                                                          tttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTTCTCCACAACCTTTGTCATGCCGAGTACAAGAGGGGGATGTTGGGGTAAAAAGGGTAACGCCTGCGCCGATTGCGCTGAAACCAGCCGGt                                                     >  2:944295/1‑150 (MQ=255)                                                                                           ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                    <  2:601750/150‑1 (MQ=255)                                                                                             cGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATTCAGGGCACAAGAGGGGGATGTTGGTGGAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCg                                                   >  1:482016/1‑149 (MQ=255)                                                                                                 aaCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAAAGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACtt                                               >  2:87452/1‑149 (MQ=255)                                                                                                   ccATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAt                                           <  1:613960/151‑1 (MQ=255)                                                                                                    cATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATa                                          <  1:178813/151‑1 (MQ=255)                                                                                                    cATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCCGTGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCa                                            >  2:783367/1‑149 (MQ=255)                                                                                                    cATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGTGTAAAGCCTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCAt                                           >  1:906515/1‑150 (MQ=255)                                                                                                          tGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGTGGAAAAAGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTTCATAGttt                                      >  1:949589/1‑149 (MQ=255)                                                                                                              gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTcc                                 >  2:507058/1‑150 (MQ=255)                                                                                                              gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCACATCATTGGTCCTGCCGTGCACAAGAGGGGGATGGTGGGGGAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAAACCACCAGGTCGACTTTGATAGTTTTCTcc                                 >  1:400749/1‑150 (MQ=255)                                                                                                                tGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTg                               >  2:567791/1‑150 (MQ=255)                                                                                                                  cACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCCCAGATTGCGGTGAACCCAGCAGGTCGGCTTGCATAGTTTTCCCCTGcc                             >  2:294141/1‑150 (MQ=255)                                                                                                                     cGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGGCATGCCGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCCGCAGGGCGACTTGCATAGTTTGCTCCTGcc                             >  1:194227/1‑147 (MQ=255)                                                                                                                       aGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGGCTTGCATAGTTTGCTCCTGCCATaa                         >  1:339489/1‑149 (MQ=255)                                                                                                                        gcgcCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGt                      <  1:571989/151‑1 (MQ=255)                                                                                                                         cgcCAGCAAGGTATTGGCGGTAAAGGGTTGTACCCCATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGGGGAAAAGGGGGAACGCCTGCGCAGATTGCGCTGAAACCAGCAGGGCGACTTTCATAATTTGCTCCTGCCCTAACg                       >  1:858500/1‑149 (MQ=255)                                                                                                                                aaGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGGCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGCCAAAACGGgaagaa                >  1:99371/1‑149 (MQ=255)                                                                                                                                    tATTGGCGGTAAAAGGTTGCACCACATCCTTTGTCATGCAGTGCACAAGAGGGGGATGTTGGTGGAAAAAGAGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCaa            >  1:562823/1‑149 (MQ=255)                                                                                                                                     aTTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGa         <  2:457809/151‑1 (MQ=255)                                                                                                                                        ggCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGAcc       <  1:104860/150‑1 (MQ=255)                                                                                                                                         gCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtg     <  2:803657/151‑1 (MQ=255)                                                                                                                                            gTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtggtg  >  2:672358/1‑151 (MQ=255)                                                                                                                                            gTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtggtg  <  1:256481/151‑1 (MQ=255)                                                                                                                                            gTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtggtg  <  2:504973/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                ACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTG  >  NC_000913/2185053‑2185342

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