Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG40.4% 68.3 / 10.6 47N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (7/21);  new base G (19/0);  total (26/21)
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

CAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCG  >  NC_000913/2677891‑2678177
                                                                                                                                                     |                                                                                                                                          
caGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGttt                                                                                                                                          <  1:311318/151‑1 (MQ=255)
caGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGgtt                                                                                                                                          >  1:439909/1‑151 (MQ=255)
        tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGcc                                                                                                                                    <  1:167958/149‑1 (MQ=255)
         ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCAc                                                                                                                                  <  2:29889/150‑1 (MQ=255)
                 ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACCGACCCAGGCAGAGGGTTTTCCCCACCGTCAGcc                                                                                                                          >  2:866293/1‑150 (MQ=255)
                    ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACCGACCGCGGCAGAGGGTTTTCGCCCCTGTCCGCCCCTc                                                                                                                      >  1:765914/1‑151 (MQ=255)
                         gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTccacc                                                                                                                  <  2:411547/150‑1 (MQ=255)
                             aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCCGACCGCGGCCGAGGTTTTCCCCCCCTCTCCGCCCCTCCCCCCGACt                                                                                                             >  2:736416/1‑151 (MQ=255)
                               aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTCAGCGGGGCAGAGGTTTTTCGCCACCCTCAGCCCTTccccc                                                                                                                  >  2:698798/1‑144 (MQ=255)
                               aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCCGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCCTCCACCCGAcct                                                                                                            >  1:276354/1‑148 (MQ=255)
                                cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACCGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCCCCAGACTTa                                                                                                           >  1:910300/1‑150 (MQ=255)
                                  aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACt                                                                                                        <  2:327044/151‑1 (MQ=255)
                                  aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACt                                                                                                        <  1:251159/151‑1 (MQ=255)
                                         cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATa                                                                                                 <  1:712219/151‑1 (MQ=255)
                                          aGGACTAACGTCAACGCCCAGCGCCTGTAAGCACCGTTGGGCATCCTGCCCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGCCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGAt                                                                                                  <  2:224860/149‑1 (MQ=255)
                                          aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGACCGAGGGAGAGGGTTTCCGCCCCTGTCAGCCCCTCCACCCGaacctaatcctggta                                                                                                 >  2:499144/1‑136 (MQ=255)
                                            gACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAc                                                                                              <  1:466647/151‑1 (MQ=255)
                                               tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGC‑CCTTAACCGGGGAGAGGTTTTCCCCCCCCCGTGGCCCTTTCACCCGCCTTAAATCCTTATAAAcccc                                                                                           >  1:51924/1‑149 (MQ=255)
                                                aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATGTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTAAGCCCTTCCACCAGACTTAACTCCTGATAAACcg                                                                                            <  2:639588/149‑1 (MQ=255)
                                                   gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGCCTCACC‑AGGCCGAGGTTTTTCCCCCCCTTCCGGCCCTCCACCCCACTTAACACCTCCTaaaaacccc                                                                                         >  2:763216/1‑143 (MQ=255)
                                                    tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTCACCCAGGCAGGGGGTTTTCGCCACCGTCCGCCCCTCCACCAGACTTAACTCCCGATAAAAcgcgcg                                                                                        >  2:330420/1‑149 (MQ=255)
                                                            cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGGAGGGGGTTTCCGCCACTGTCCGCCCTTCCACCCGACTTAACTCCTGATAAAACCCGCGCAcccccc                                                                                >  2:793886/1‑148 (MQ=255)
                                                                    gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGCCCCCGGCAGAGGGTTTTCGCCCCCGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAACCCCGCGCACCCCCCGTTCAGCgg                                                                       >  2:763840/1‑150 (MQ=255)
                                                                     tAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGACCGAGGCCGAGGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCCCGCGCACCCCCCGGTCAGCgg                                                                       >  1:243880/1‑149 (MQ=255)
                                                                      aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                    <  2:320467/151‑1 (MQ=255)
                                                                       aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACg                                                                   <  2:792247/151‑1 (MQ=255)
                                                                        gCAACGTTGGGCATCCTGGGCTATTTGCAGGTAATCAACCAGGCCGTTGCGGCGGGGCC‑ACTGACCGAGGAAGAGTTTTTCCGCCACTGTCAGCCCTTCCAGCAGCCTTAACTCCAGATAAACCGCGCGCACCCCCAGTTAAGCGGCACgg                                                                  <  2:388498/151‑1 (MQ=255)
                                                                              ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGACCGGGGGAGAGGGTTTTCCCCCCCGTCAGCCCCTCCCCCAGACCTAAATCCTGGTAAAACCCCCGCACCCCCCGGTCaac                                                                         >  2:699054/1‑136 (MQ=255)
                                                                              ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑AACGCCCCGGGGAGAGGGTTTCCGCCCCCGCGATCCCCTCCCCCCGAACTAAATCCTGATAAACCCCGCGCCGCCCCCCGTCAGCAGCACGGCGAg                                                              >  1:447428/1‑149 (MQ=255)
                                                                                  gCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACCGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAAATCCCTATAAAACGCGCGCCCCCCCCGTTCCGCGGCACGGCGAGGCAGc                                                         >  1:19044/1‑150 (MQ=255)
                                                                                   cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGACCGAGGGAGAGGGTTTCCGCCCCTGTCAGCCCTTCCCCCCGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                           >  1:805769/1‑147 (MQ=255)
                                                                                      ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                      <  2:346279/149‑1 (MQ=255)
                                                                                         gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGAGCGAGGGAGAAGGTTTCCGCCCGAGTTAATCCGTCTTCCCGACTTAACTACTGGTAAAACCCGCCCACCCCCCGTGCAGACGCGCAGCGGGTCCCCGTAGCtt                                                  >  1:752641/1‑150 (MQ=255)
                                                                                                 tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑CCTGACCGAGGCAGGGGGTTTTCACCACTGTCAGGCCTTCCACCACACTTAAATCCTCATAAAACCCGCGCACCCCCCCCTCAGCCGGACGGGGAGCCAGCGGAGCTTCGTc                                              >  2:583334/1‑146 (MQ=255)
                                                                                                 tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCCGCCCCTCCACCAAACTTTACTCCTGTTAAAACGCGCGCACCCCCCGCTCAGCGGCACGGCGAGTGACAGTAGCTTCCCCAcc                                           >  1:976574/1‑149 (MQ=255)
                                                                                                  gCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACCCACCGAGGCCGAGGGTTTTCGCCACTTTCAGCCCCTCCCCCAGACATAACTCCTGATAAAACACGCGCACCCCCCCCTCCGCGGCACGCCGGGTCAGCGTAGCGTCCTTACCCTc                                        >  2:74136/1‑151 (MQ=255)
                                                                                                  ccAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTCCACCAGGAAGCGGTTTTTCGCCACTGTTGTCCCCTCCACCACCCTTAAATCCTGATAAACCGCACCCACCCCCCGCTCACCGGCAAGGCGCGTCCGCGGATCTCCCTCAcc                                           >  2:270016/2‑148 (MQ=255)
                                                                                                    aGGTAATCAATCATGCCGTTGCGGCGGGGGC‑ACTCAACGAGGCAGAGGGTTTCCGCCACGGTCCGCCCCTCCCCCCCACTTCAAACCTGATAAACCCCGCGCACCCCCCGTTCAGCGGCCCGGCGCGTCTACGGAGCCTCTTCGCCCTcc                                       >  1:366680/1‑150 (MQ=255)
                                                                                                         atcaatcaTGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                 <  1:472095/151‑1 (MQ=255)
                                                                                                               caTGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCg                           <  1:654345/151‑1 (MQ=255)
                                                                                                                aTGCCGTTGCGGCGGGGCC‑ACCGCCCTAGGCAGAGGGTTTCCACCCCTGCCAGCCCTTCCACCAGACTTAACTCCTGATCAACCCCGCCCACCCCCCAACCAGCGGCACGGCGGGGCAGCGGAGATTATCCACCTTCCCGACGCGTCCCGc                          >  2:555778/1‑151 (MQ=255)
                                                                                                                            cggGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCCCCTTCCAGTCGCGTCTCGTCAATCCAGatag              <  2:549780/151‑2 (MQ=255)
                                                                                                                            cggGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat              <  2:909962/151‑1 (MQ=255)
                                                                                                                                  c‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc        <  1:248853/151‑1 (MQ=255)
                                                                                                                                       gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc    <  2:977642/151‑1 (MQ=255)
                                                                                                                                        aCCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATcc   >  1:725758/1‑151 (MQ=255)
                                                                                                                                        aCCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATcc   <  1:696494/151‑1 (MQ=255)
                                                                                                                                        aCCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATcc   >  1:43026/1‑151 (MQ=255)
                                                                                                                                         ccGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg  <  2:439683/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
CAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCC‑ACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCG  >  NC_000913/2677891‑2678177

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