Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,5090AT45.8% 43.7 / 26.7 48T236S (ACC→TCC) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base A (8/18);  new base T (22/0);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 3.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCC  >  NC_000913/3745362‑3745652
                                                                                                                                                   |                                                                                                                                               
cTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCCGTTTTccc                                                                                                                                              >  2:168916/1‑151 (MQ=255)
       ggtggtggCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTAcc                                                                                                                                         <  1:335815/149‑1 (MQ=255)
          ggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGc                                                                                                                                    <  2:382244/151‑1 (MQ=255)
            tggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGGGCTGTTTACCCCTACCGAAGc                                                                                                                                    >  1:291531/1‑149 (MQ=255)
               cAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGGCCTCTTTTCCCCCTCCCAAGCCgggg                                                                                                                               >  2:498340/1‑151 (MQ=255)
                 ggCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGGCCGTTTTCCCCTACCCAAACCgggg                                                                                                                               >  1:538038/1‑149 (MQ=255)
                             tAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCCCCGGCCCGTTTTCCCCCACCCGAGCCCGGGCCGGTTCTGCttt                                                                                                                  >  2:572391/1‑150 (MQ=255)
                             tAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGCTTTTCTCCCGCCTGTTTTCCCCCCCCGAAGCCcggggc                                                                                                                             >  2:340546/1‑137 (MQ=255)
                                     cccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTTCCCCTACCCAAGTCGGGGGCGTTGCTGCTTTTTTTgcgc                                                                                                          >  2:285085/1‑150 (MQ=255)
                                       cGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGCCTCCGGCATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGGTATCTCCGGCCCCGTTTCCCCGCCCGAAACCCGGGAGGGTTCTGCCTTCTATGCGgtgt                                                                                                       >  2:380924/1‑151 (MQ=255)
                                          cAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTg                                                                                                    <  2:491211/151‑1 (MQ=255)
                                           aGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTTTCCCTCCCCAAACCGGGGCCGTTGCTGCTTTCTTTTCGCTGttt                                                                                                     >  2:102007/1‑149 (MQ=255)
                                             cAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTGTTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCg                                                                                                 >  2:520296/1‑151 (MQ=255)
                                             cAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCTCTCCGGGCCGTTTTCCCCTACCCAAACCGGGGGCCGTGGTGGTTTCTTTGCGCGTTTTGt                                                                                                   >  2:448991/1‑149 (MQ=255)
                                              aaaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTTTCCCGGCGTGTTTTCCCCTACCCAAGACGGGGCCGGTGCTGCTTTCTATTCGCCGtttt                                                                                                    >  2:452116/1‑146 (MQ=255)
                                               aaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTTCTCCCGCCTGTTTTCCCCCACCCAAACCGGGGGCGGTTCTGCTTTCTTTTCGCTGTTTTTCGcc                                                                                               >  2:229674/1‑151 (MQ=255)
                                               aaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCCCTACCGAAGCCGGGGGCGGTTCTGCTTTCTTTGCGCTGTTTGTCGcc                                                                                               >  1:245337/1‑151 (MQ=255)
                                                   gCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAg                                                                                           <  2:601108/151‑1 (MQ=255)
                                                        gATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGCTTTTCTTCAGGCTGTTTTACCCTGCCCAAAGCCGGGCCGGTGCTAGTTTCTTTTCCCTTTTTTTCGCcaaaaatat                                                                                       >  2:275314/1‑143 (MQ=255)
                                                                    cTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCCGCCTGTTTACCCCTACCGAAGCCGGGGCCGGTGCTGCTTTCTATGCGCTGTTTGTTGCCACAGTTATTTACCGTGaaa                                                                            >  1:597943/1‑149 (MQ=255)
                                                                     tGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTTTATGCGCTGTTTGTTGCCACAGTTATTTACCGTGAAAt                                                                           >  1:361164/1‑149 (MQ=255)
                                                                      ggCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGAc                                                                        <  1:399522/151‑1 (MQ=255)
                                                                       gCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCTCTCCGGCCTGTTTTCCCCCACCCAAGACCGGGGGGTTGCTGGCTTCTTTGCCGTGTTTTTCTCCCCACTTTTTTACCGGGAAATTAcc                                                                       >  2:107220/1‑151 (MQ=255)
                                                                              tttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGcc                                                                 <  1:102026/150‑1 (MQ=255)
                                                                               ttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCTGTTTTCCCCCCCCCAAGCCGGGGCCGCTGCTGCTTTCTATTCGCTGTTTTTCGCCCCCGGTTTTTTCCGTGAAATGACCTTCGcccc                                                               >  2:122863/1‑149 (MQ=255)
                                                                               ttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCCGCCTGTTTTCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTTGCCACAGTTATTTTACGTGAAAAGAACTTCGCcac                                                               >  1:489627/1‑151 (MQ=255)
                                                                                tGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTTCCCCACCGAAGCCCGGGCCGTTGCTGCTTTTTATGCGCTGTTTTTTGCCACAGTTATTTACCGTGAAATGACCTTCGCcac                                                               >  1:621883/1‑150 (MQ=255)
                                                                                    tcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTTTCGCCACAGTTATTTTCCGTGAAATGACCTTCGCCACACctt                                                          >  1:310771/1‑151 (MQ=255)
                                                                                      cGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTgg                                                        <  1:198941/151‑1 (MQ=255)
                                                                                       ggTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCTCCGGCCTGTTTTCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTTTTTGCGCTGTTTGTCGCCCCAGTTATTTACCGGGAAATGACCTTCGCCACACttt                                                          >  1:186383/1‑148 (MQ=255)
                                                                                         tATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCCCCACCGAAGCCGGGGCCGGTGCTTCTTTCTTTTCGCCGTTTGTCGCCACAGTTATTTTCCGTGGAATGACCTTTGCCACACTTTTGCAt                                                     >  1:223864/1‑151 (MQ=255)
                                                                                         tATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCCGTTTACCCCTACCGAAACCGGGGCCGTTGCTGCTTTCTATTCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTTGCAt                                                     >  1:615432/1‑151 (MQ=255)
                                                                                         tATCTGGGCGCTGTTGCTTCCAGTCATCATTATTGGCGGGTTTCTCTGCGGCCTGTTTTCCACCCCCCAAGCCCGGGCCGTTGCTGCTTTTTATCTGCTGTTTTTCGGCCCAGTTTTTTACCCTGAAATGACCTTTGCCCCCCTTTgtgct                                                     >  2:360293/1‑147 (MQ=255)
                                                                                                 cgcTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTACACCGGCCCGGTTTCCCCTACCGAAACCGGGGGCGTTTCTGCTTTTTAAGCGCCGGTTTTTGCCAACCTTATTTTCCGGGGAATGACCTTTGCCACACTTTTGGATgata                                                 >  1:531395/1‑144 (MQ=255)
                                                                                                            ccAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCgcggcg                                   <  1:78254/150‑1 (MQ=255)
                                                                                                            ccAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCgcggcg                                   <  2:434227/150‑1 (MQ=255)
                                                                                                             cAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGaa                                 <  1:415262/151‑1 (MQ=255)
                                                                                                              aGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGaaa                                <  2:574676/151‑1 (MQ=255)
                                                                                                                tcatcaTTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGaaaa                               <  2:567963/150‑1 (MQ=255)
                                                                                                                  atcatTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAacca                            <  1:133838/151‑1 (MQ=255)
                                                                                                                  atcatTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAacca                            <  1:452261/151‑1 (MQ=255)
                                                                                                                              gggTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTa                <  2:537864/151‑1 (MQ=255)
                                                                                                                                    cGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCCGTGGTTATGTTTc          >  1:329559/1‑151 (MQ=255)
                                                                                                                                    cGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTACTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTc          >  2:237594/1‑151 (MQ=255)
                                                                                                                                     gCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCt         <  1:308143/151‑1 (MQ=255)
                                                                                                                                        ccGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCtggt      <  1:99745/151‑1 (MQ=255)
                                                                                                                                          ggCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCtggtgg    <  1:349408/151‑1 (MQ=255)
                                                                                                                                            ccTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGcc  >  1:99479/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCC  >  NC_000913/3745362‑3745652

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