Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG44.7% 62.8 / 11.4 47N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (7/19);  new base G (21/0);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 9.93e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677892‑2678184
                                                                                                                                                   |                                                                                                                                                 
aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGAGGgttt                                                                                                                                                >  1:978087/1‑151 (MQ=255)
   tgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGAGGgttt                                                                                                                                                >  1:434462/1‑148 (MQ=255)
          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)
                                                                                                                                                   |                                                                                                                                                 
AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677892‑2678184

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: ‑