Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,255	0	A	G	25.4%	81.6 / 10.2	56	Q152R (CAG→CGG)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (28/13);  new base G (14/0);  total (43/13)
Fisher's exact test for biased strand distribution p-value = 2.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCA  >  NC_000913/3113‑3393                                                                                                                                                |                                                                                                                                           gcggcgCTGATGGCGATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCa                                                                                                                                            >  2:399291/1‑143 (MQ=255)       ccgATGGCGATGAATGCACACTGCGGCAGGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTtgat                                                                                                                                      <  1:1131142/141‑1 (MQ=255)         gATGGCGATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTtgatga                                                                                                                                    <  1:802373/143‑1 (MQ=255)          aTGGCGATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTgatgat                                                                                                                                   >  1:879985/1‑143 (MQ=255)              cGATGAATGAGCCGGGCGGCAAGCCGCTGCAGTGCCA‑TCGTTTGCTGGCTTTTCTGCGCGAGCTGGACGGCCTTATCTCCGGCCGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCgaa                                                                                                                               <  1:409481/143‑1 (MQ=255)               gATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGAACgaag                                                                                                                              >  2:820461/1‑143 (MQ=255)               gATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATCCAATTGATGATCgaag                                                                                                                              >  2:546694/1‑143 (MQ=255)               gATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGGGGTATGCAGTTGATGATCgaag                                                                                                                              >  1:1049955/1‑143 (MQ=255)                  gaatgaaCACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGaaa                                                                                                                           >  2:608010/1‑143 (MQ=255)                       accacTGCGGCAAGCCGCT‑TAATTATACTCGTTTGCTGTCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGAc                                                                                                                      <  1:476980/141‑1 (MQ=255)                          acTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGGTGGAAGGCCGTATCGCGGGCAGCATGCAGTACGACAACGTGGCACCGTGGTTTCGAGGGGGTGTGACGTTGAGGACGGAAGAAGACGAcgtc                                                                                                                   >  1:983699/1‑143 (MQ=255)                          acTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTGCGGTTGATGATCGAAGAAAACGAcctc                                                                                                                   >  2:1127573/1‑143 (MQ=255)                             gCGGCAAGCCGCT‑GAATGACACTCGTGTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGTACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACatcatc                                                                                                                <  1:138164/143‑1 (MQ=255)                                     ccGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTCATGATCGAAGAAAACGCCATAATcagccagga                                                                                                        >  2:788744/1‑141 (MQ=255)                                        cT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTAGGCGGTTGATGAACGAAAAACCCGACAACATCAGGCAGCAAGt                                                                                                     >  2:4674/1‑143 (MQ=255)                                         t‑tAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTg                                                                                                    >  1:492322/1‑143 (MQ=255)                                             aTGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCa                                                                                                 >  1:1052012/1‑143 (MQ=255)                                                acacTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTAGGCAGTTCATTAAAGAAAAAAAGGACGTCATGCACCAGCAAgagacacgg                                                                                              >  2:116285/1‑135 (MQ=255)                                                   cTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGGTGATCGAAGCAAACGACATCATCAGCCAGCCAGTGCCCgggcct                                                                                           >  1:407662/1‑140 (MQ=255)                                                      gTTTGCTGGCTTTGATGGGCGGGTTGGAAGGCCGTATCTCCGGCAGCATTCATTCCTACATCGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTtgat                                                                                        <  1:654356/143‑1 (MQ=255)                                                            tGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTgg                                                                                  >  1:112133/1‑143 (MQ=255)                                                              gCTCTGATGGGCGAGCTGGCAGGTCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAACCGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCt                                                                                <  1:597170/143‑1 (MQ=255)                                                                  ggATGGGCGAGCTGCAAGGCCCGATCTGCGGCATCATTCATTACGACAACGTGGCACCGTGTTTTCTCGTTGGTATGCAGTTGATCATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTgg                                                                            <  1:895611/142‑1 (MQ=255)                                                                  gactgcgcgcccTGGCAGGCCGTGGCTCGGGCAGATTTCATTACTACAACCAGCGACCGTCTGTTCACGGTGGCATGCAGTTGATGTTCGAAGAACGAGACAGCATGAGCCAGCAAGTGCAAGGGTTTGATGAGTGGCTGTgg                                                                            <  2:850012/132‑1 (MQ=255)                                                                        gCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGAAAGTGCGTTTCAGGGTTTGATGAGTGGCTGTGGGTGCTg                                                                      >  1:210800/1‑143 (MQ=255)                                                                          gAGCTGGAAGGCCGTATCTCCGGCAGCATTGATTACGGCGACGTGGCACCGTGTGTGCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGc                                                                    <  1:1024586/143‑1 (MQ=255)                                                                          gAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACGTCATCAGCCAGCAAGTGACAGGGTCTGAAGAGAGGCTGTGGGAGATGGc                                                                    >  2:733920/1‑143 (MQ=255)                                                                             cTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATCCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTa                                                                 >  2:158197/1‑143 (MQ=255)                                                                              tGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTAGCCGGTTGATGAACGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGGTGAGTGGCAGCGGGTGCgggactat                                                                >  2:313544/1‑138 (MQ=255)                                                                                gAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTGCAGTTGATGATAGAAGAAAACGACATCATCAGCCAGCCAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATct                                                              >  2:267852/1‑142 (MQ=255)                                                                                gAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCGGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGAAGAGTGGCTGTGGGTGCTGGCGTAGcc                                                              >  2:254656/1‑143 (MQ=255)                                                                                gAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGCTGATCGCAGAAAAAGTGGTGGTCAGCCAGCAAGTGCCAGGGTTTGATGAGGGGCGGTGGGTGCTGGCGTATcc                                                              >  2:816011/1‑143 (MQ=255)                                                                                gAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGGTTTCTCGGTGGTGTGCACTTGGGGATAGGAGAAGACGACGGCAGCAGCCAGCAAGTGGCAGGGTTTGATGAGTGGCTGTGGGCTCTGGCGTAttc                                                              >  1:432300/1‑141 (MQ=255)                                                                                 ggggCCGGAGCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGAGCGAAGAAGACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCg                                                             <  1:911984/141‑1 (MQ=255)                                                                                    gCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCgggg                                                          >  1:733742/1‑143 (MQ=255)                                                                                      cGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCggggct                                                        >  1:509193/1‑141 (MQ=255)                                                                                          tctcCGGCATCATTCCTGACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGCTGATCGAAGAATACGACATCATTAGCTAGCAAGTGCTAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTaaa                                                    <  1:899517/143‑1 (MQ=255)                                                                                            tcCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTGCGGTTCATTATCAAAGAAACCGACATCATCGGCCAGCAAGTGCCAGGGTTTGATGAGCGGCTTTGGGGTCTGGCGCACGCGGGGacccaagc                                                  >  2:646294/1‑136 (MQ=255)                                                                                              cGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTGCGGTTGATTATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGCGGGTGCTGGCGTACCCGGGGATTcaacacc                                                >  2:590657/1‑139 (MQ=255)                                                                                                 cagcaTTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAGAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATAAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTGACGTGTCGa                                             >  1:677308/1‑143 (MQ=255)                                                                                                    cattcattACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTCCGGTTCATTAACGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTACAGTCTCGACgg                                          >  2:714157/1‑143 (MQ=255)                                                                                                    cattcattACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACgg                                          >  1:822327/1‑143 (MQ=255)                                                                                                     attcattACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTCCGTATCATTAACGAAGAAAACGACATCCTCAGCCAGCAAGCGCCAGGGCATGATGAATGGCGGTGCGAGCTCGCGTATCCGGGGAACAAAGtccagccata                                         >  2:120444/1‑135 (MQ=255)                                                                                                      ttcattACGACAACGTGGCACCGTGTTTTCTCGGAGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAACGTCTCGAAGGCa                                        >  1:154213/1‑143 (MQ=255)                                                                                                         attACGACAACGTGGCACCGTGTTTTCTCGGTGGTCTGCGGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTAGACGGccggg                                     >  2:392548/1‑139 (MQ=255)                                                                                                          ttACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCGGTTCATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAg                                    >  2:444395/1‑143 (MQ=255)                                                                                                                cAACGTGGCACCGTGTTTTCTCGGTGGTCTGCGGTTCATTATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGCGATTACAGTCTCGACGGCAGAAGCCAggg                              >  2:103640/1‑143 (MQ=255)                                                                                                                   cGTGGCACCGTGTTTTCTCGGTGGTCTGCGGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGGCTCGACGGCAGAAGCCAGGGCTa                           >  2:640285/1‑143 (MQ=255)                                                                                                                   cGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTa                           >  1:1149149/1‑143 (MQ=255)                                                                                                                   cGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTACCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTa                           >  2:38093/1‑143 (MQ=255)                                                                                                                          ccGTGGTTTCTCGGGGGTATGCAGTGGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGGGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTAcc                    <  1:788744/143‑1 (MQ=255)                                                                                                                             tgtTTTCTCGGTGGTCGCCGGATGATTAAAGAAGAAAACGACATCATCAGCCAGCAAGTGCGAGGGTATGAAGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATGAAAGTCTCGACGGCAGAAGCCGGGGCTAGTTGACCGgc                 >  2:594055/1‑143 (MQ=255)                                                                                                                               ttttCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCACACGCCAGGGCTATTTTACCGgcgc               >  1:297269/1‑143 (MQ=255)                                                                                                                                 ttCTCGGTGGTCGCCGGATCATTATCGAAAAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAg             >  2:161626/1‑143 (MQ=255)                                                                                                                                    tcGGTGGTCGGCAGATCATTAAAGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTAGATAAGTGGCTGAGGACGCGGTCGTATCGCGGGAACAAAATCTCGACGGCAGCAGCCAGGGCGATTTTACCGGCGAAGTAt          >  2:510338/1‑143 (MQ=255)                                                                                                                                     cGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTGTCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATc         >  1:588752/1‑143 (MQ=255)                                                                                                                                            aTGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCGGTATCGCCGCCa  >  2:437517/1‑143 (MQ=255)                                                                                                                                                |                                                                                                                                           GCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCT‑TAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCA  >  NC_000913/3113‑3393

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