Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,4940CT41.5% 35.9 / 17.7 42R231C (CGC→TGC) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (9/15);  new base T (17/0);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 3.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTC  >  NC_000913/3745362‑3745644
                                                                                                                                    |                                                                                                                                                      
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                           <  2:434227/150‑1 (MQ=255)
                                                                                                            ccAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCgcggcg                           <  1:78254/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:452261/151‑1 (MQ=255)
                                                                                                                  atcatTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAacca                    <  1:133838/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)
                                                                                                                                    |                                                                                                                                                      
CTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTC  >  NC_000913/3745362‑3745644

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