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

GCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCC  >  NC_000913/2677902‑2678190
                                                                                                                                              |                                                                                                                                                   
gcgcGGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACt                                                                                                                                              <  1:411760/149‑1 (MQ=255)
       gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCC‑TGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGccc                                                                                                                                      >  1:665616/1‑150 (MQ=255)
       gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCC‑TGACCCAGGCGGGGGTTTTTCGCCACTGGCAGccc                                                                                                                                      >  2:519027/1‑150 (MQ=255)
             cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGC‑CCGCCGAGCCAGAGGGTTTTCGCCCCCTTCAGCCCCTccccc                                                                                                                               >  2:53368/1‑151 (MQ=255)
               ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCAC‑CGACCGAGGCGGAGGGTTTTCGCCACTGTCAGCCCCTccacc                                                                                                                               >  1:23135/1‑149 (MQ=255)
                      aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGAACGCGGCAGAGGGTTTCCCCCCCTGTCAGCCCTTTCACCCGaact                                                                                                                         >  1:197609/1‑145 (MQ=255)
                              cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATa                                                                                                              <  1:915138/151‑1 (MQ=255)
                               aGGACTAACGCCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGCCCGAGGCCGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTg                                                                                                                 <  1:99230/147‑1 (MQ=255)
                                   cTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCACCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAcc                                                                                                          <  1:583395/150‑1 (MQ=255)
                                         tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGGCCCAGGGCGGGGGTTTTCGCCACCGTCAGCCCCCCCCCCCCACTTAACTCCCGATaaaaccc                                                                                                        >  2:493726/1‑142 (MQ=255)
                                               cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCC‑CGAACGAGGCAGAGGGTTTTCGCCACTTTCAGCCCCTCCCCCACACTTAACTCCTGATAAACCGCGCGCAccccc                                                                                              >  1:558304/1‑150 (MQ=255)
                                               cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑CGACCGAGGCAGAGGGTTTTCGCCCCCGTCAGCCCTTCCCCCCGACTTTACTCCCGATAAAACGCGCGCAccccc                                                                                              >  1:833816/1‑150 (MQ=255)
                                                  aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑CGGCCGAGGCAGAGGGTTTTCCCCCCCGTCAGCCCTTCCCCCCGACTTAACTCCCGATAAAACGCGCcccccccccctt                                                                                          >  2:329254/1‑147 (MQ=255)
                                                   gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGGCCGAGGCAGAGGGTTTTCGCCCCCCTCAGCCCTTCCACCAGACTTAAACCCTGATAAAAAGCGCGCAccccccgctc                                                                                         >  2:803653/1‑146 (MQ=255)
                                                   gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                         <  1:669284/151‑1 (MQ=255)
                                                    cgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCa                                                                                        <  2:528695/151‑1 (MQ=255)
                                                      ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCCAGGGAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAAATCCTGATAAAACGCGCGCACCCCCCGtt                                                                                          >  1:608233/1‑147 (MQ=255)
                                                         gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                     <  2:216361/149‑1 (MQ=255)
                                                            aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCC‑TGACCGAGGCAGAGGGTTTTCGCCCCTCTCAGCCCCTCCACCCGACCTAACTCCCGATAAACACCGCGCACCCCCCCTTCAGCGGccc                                                                                 >  2:599296/1‑148 (MQ=255)
                                                              cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGc                                                                              <  1:675670/151‑1 (MQ=255)
                                                                aCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCC‑CGCCCCAGGCAGAGGGTTTTCGCCCCCGTCCGCCCCTCCCCCAGACTTTACTCCTGATAAACCGCGCGCccccccccctcagc                                                                                      >  2:517372/1‑133 (MQ=255)
                                                                    tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTc                                                                         <  1:854983/150‑1 (MQ=255)
                                                                    tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTc                                                                         <  2:801328/150‑1 (MQ=255)
                                                                       gCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCC‑TCACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGAGTTTAATCCTGATTAAACGCGCGCACCCCCCGGTCAGCGGCACGGGGAGTc                                                                         >  1:498357/1‑147 (MQ=255)
                                                                        cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGTGCCCCTTTCCCAAGGAAGGGGTTTT‑CGCCACCGTCAGCCCTTTCCCCACACTTAAATCCCGGTAAAACGCCCGCACCCCCAGTTCAGCGGCACGGCggg                                                                           >  1:706249/1‑142 (MQ=255)
                                                                           ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAg                                                                  <  2:334456/150‑1 (MQ=255)
                                                                                 ccATTTGCAGGCAATCAATCATGCCGTTGCGGCGGGGCCCC‑TGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCCTCCACCAGACTTAAATCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGCGTCCGCGTAGCtt                                                               >  1:667221/1‑147 (MQ=255)
                                                                                    tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCC‑CGACCGAGGCAGAGGGTTTTCGCCACCGTCAGCCCTTCCACCCGACTTAAATCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGGACGGCGAGGCAGCGTAGCTTCGTcccc                                                        >  2:63720/1‑151 (MQ=255)
                                                                                          ggTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTTCGCCCCTGTCAGCCCTTCCCCCAGACTTAAATCCCGATAAAACGCGCGCACCCCCCGTTCAGCGGCCCGGCGGGTCAGCGTAGCCTCGTCACCCTcc                                                    >  2:359490/1‑149 (MQ=255)
                                                                                                  atcaTGCCGTTGCGGCGGGGCCCC‑CGACCGAGGCAGAGGGTTTTCCCCCCCGTCAGCCCCTCCCCCCGACTTTACCCCTGCTAAAACGCGCGCCCCCCCCCCTCAGCGGCACGGCGGGGTAGCGGAGGCTCGTCACCTTcc                                                    >  2:692589/1‑141 (MQ=255)
                                                                                                              cggcggGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                                 >  1:31538/1‑148 (MQ=255)
                                                                                                                gcggGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                             >  1:507574/1‑150 (MQ=255)
                                                                                                                  ggggCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                           >  2:867095/1‑150 (MQ=255)
                                                                                                                  ggggCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                          >  1:72059/1‑151 (MQ=255)
                                                                                                                      ccAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcc                       <  2:634968/150‑1 (MQ=255)
                                                                                                                        aC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                    >  2:897000/1‑151 (MQ=255)
                                                                                                                              ccGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACGCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg               <  2:124036/151‑1 (MQ=255)
                                                                                                                                    cAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGtt         <  1:696124/151‑1 (MQ=255)
                                                                                                                                      gagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg       <  2:804235/151‑1 (MQ=255)
                                                                                                                                            ttCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGcc  <  2:538162/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAC‑TGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCC  >  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 < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑