Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0430CT43.8% 44.6 / 25.5 48E107K (GAA→AAA) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (8/19);  new base T (21/0);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCC  >  NC_000913/2677902‑2678190
                                                                                                                                             |                                                                                                                                                   
gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTg                                                                                                                                             <  1:148974/150‑1 (MQ=255)
       gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCt                                                                                                                                     <  2:206049/151‑1 (MQ=255)
        cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGGAGAGGGTTTCCGCCACTGTCAGcccct                                                                                                                                    >  2:685036/1‑149 (MQ=255)
           aGCGTTGAAACAAGCTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGCCCCCGGGGGAGGGTTTCCGCCCCTATCAACCCCTc                                                                                                                                   >  1:417773/1‑149 (MQ=255)
            gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTcca                                                                                                                                 <  1:813417/150‑1 (MQ=255)
              gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAAGCAGAGGGTTTTCGCCACTGTCCGCCCTTccccc                                                                                                                               >  2:880682/1‑150 (MQ=255)
                tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGCGGCGGAGGGTTTTCCCCCCCGTCAGCCCTTccccc                                                                                                                               >  1:25714/1‑148 (MQ=255)
                 gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGa                                                                                                                            <  1:598008/150‑1 (MQ=255)
                  aaaCAAGTTGTTCAGGACTAACGTCAGCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACt                                                                                                                          <  2:1085120/151‑1 (MQ=255)
                  aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGTTTTTCGCCACTGTCCGCCCCTccccc                                                                                                                               >  1:410981/1‑146 (MQ=255)
                        gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTc                                                                                                                    <  1:550082/151‑1 (MQ=255)
                                  aCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGCAGAGGGTTTTCGCCCCTTTCAGCCCCTCCCCCAGACTTAACACCTGATaaa                                                                                                            >  2:1006456/1‑149 (MQ=255)
                                    tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGTCCCCGGGAGAGGGTTTTCGCCCCTGTCCGCCCCTCCACCCGACTTTAATCCTGATaaa                                                                                                            >  1:428273/1‑147 (MQ=255)
                                      aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGGCCGAGGGGGGGGGTTTTCGCCACTGTCCGCCCCTCCCCCCGACTTTACTccctatttaacgcg                                                                                                       >  1:1112571/1‑138 (MQ=255)
                                          cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTTAGCCCCTCCCCCCGACTTAACTCCTGATAAACCGCGCGCAc                                                                                                  >  1:999299/1‑151 (MQ=255)
                                           aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAcc                                                                                                 >  2:876380/1‑151 (MQ=255)
                                               cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGGTTTTCTCCACTGTCAGCCCCTCCACCCGACTTAACTCCTCATAAACCCCGCGCAccccc                                                                                              >  2:1061993/1‑150 (MQ=255)
                                                      ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTc                                                                                         >  2:547413/1‑148 (MQ=255)
                                                       cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCTACCCAGGCAGAGGGTTTTCGCCCCCGTTAGCCCCTCCACCCGACTTAACTCCTGATAAACCCCGCGCAccccccgctccacg                                                                                     >  2:272385/1‑142 (MQ=255)
                                                          tAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                     <  1:638245/148‑1 (MQ=255)
                                                          tAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGc                                                                                   <  1:684870/150‑1 (MQ=255)
                                                              cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGCAGAGGTTTTTCGCCACTGTCAGCCCCTCCCCCCGACTTAACTCCTGATAAAACGCCCGCACCCCCCGTTCAGCGGCACgg                                                                               >  1:443350/1‑150 (MQ=255)
                                                               aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGc                                                                              <  1:340779/150‑1 (MQ=255)
                                                                 cGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCTCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCTCCCCCAGTTCGGCGGCACGGCGAg                                                                           <  1:283430/151‑1 (MQ=255)
                                                                  gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGACTTAACTCCCGATAAACCCCGCGCACCCCCCGTTCAGCGGCACGGCGAg                                                                           >  2:462868/1‑150 (MQ=255)
                                                                    tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCCCCAGACTTTACTCCTGATAAAACCCGCGCACCCCCCCTTCAGCGGCACCGCGAg                                                                           >  2:757862/1‑148 (MQ=255)
                                                                        cATCCTGCGCCATTTGCAGGTAATCAATCATGCTGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGt                                                                    <  2:979040/151‑1 (MQ=255)
                                                                         aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGGAGAGGGTTTTCGCCACTGTCAGCCCTTCCCCCCGACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGc                                                                      >  2:802759/1‑148 (MQ=255)
                                                                           ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCCCCCCTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCCCGCACCCCCAGTTCAGCGGCACGGCCAGGCAGCg                                                                     >  2:52098/1‑147 (MQ=255)
                                                                            cTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGGCCCCGGCGGAGGGTTTTCGCCACCGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAAACGCGCGCCCCCCCCGTTCAGGGGCGCCGCGGGGCAGCGGAg                                                                  >  1:959497/1‑149 (MQ=255)
                                                                                    tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGGCGAGGGTTTTCCCCACTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCGCTCAGCGGCACCGCGAGGCAGCGGAGCCTCCTcccc                                                        >  1:1109509/1‑151 (MQ=255)
                                                                                    tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCAc                                                         <  1:641615/150‑1 (MQ=255)
                                                                                      tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGAGCGGGGCAGAGGGTTTTTGCCACCGTCAACCCTTCCACCCCACCTTACTCCTGATAAACCCCGCGCACCCCCCCTTCCGCGGGCCGGCGAGTCAGCGGAGGCTCGTCAcc                                                        >  1:1095787/1‑149 (MQ=255)
                                                                                       gCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTc                                                     <  1:55088/151‑1 (MQ=255)
                                                                                                caatcaTGCCGTTGCGGCGGGGGCCCTGACCGAGGCAGAGGTTTTTCGCCCCTGTCAGCCCTTCCCCCCCACTTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCCCCGCGAGGCAGCGGAGCTTCGTCACCCTCCAGTc                                                >  1:273243/1‑147 (MQ=255)
                                                                                                caatcaTGCCGTTGCGGCGGGGCCCCTGATCGGGGCAGAGGGTTTCCCCCACCGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGGTCAGCGGCACGGCGGGGGAGCGGAGCTTCGTCACCTTcc                                                    >  1:140590/1‑143 (MQ=255)
                                                                                                caatcaTGCCGTTGCGGCGGGGCCACTGAACCAGGCAGAGGGTTTTCCCCACTGTCAGCCCTTCCACCAGACTTAAATCCTGATAAAACCCGCGCACCCCCAGTTCAGCGGCACGGCGAGACAGCGGAGCCTCGTCACCCTCCAATcgc                                              >  2:387289/1‑149 (MQ=255)
                                                                                                    caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                         <  2:434568/150‑1 (MQ=255)
                                                                                                    caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                         <  1:1019551/150‑1 (MQ=255)
                                                                                                            tGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                                 <  1:373420/150‑1 (MQ=255)
                                                                                                             gcggcgGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAg                               <  2:273308/151‑1 (MQ=255)
                                                                                                             gcggcgGGGCCACTGAACGAGGCAGAGGTTTTTCGCCACTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAAACGCGCGCCCCCCCAGTTCAGCGGCACGGCGAGTCAGCGGAGCTTCGTCCCCTTCCAGGCGCGTCTCCCCCATcc                                 >  1:970565/1‑149 (MQ=255)
                                                                                                                           gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTa                   <  1:356900/149‑1 (MQ=255)
                                                                                                                           gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc                 >  1:691166/1‑151 (MQ=255)
                                                                                                                             ccGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg               <  1:876178/151‑1 (MQ=255)
                                                                                                                                aGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc            >  2:1112035/1‑151 (MQ=255)
                                                                                                                                    agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc        >  1:661534/1‑151 (MQ=255)
                                                                                                                                          tttCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGcc  >  2:270980/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCC  >  NC_000913/2677902‑2678190

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