Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,285,8530GT42.9% 62.8 / 24.1 42A126S (GCC→TCC) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base G (2/22);  new base T (18/0);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 5.37e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285703‑3286002
                                                                                                                                                       |                                                                                                                                                      
gggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTg                                                                                                                                                        <  2:703824/151‑1 (MQ=255)
       taATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag‑cc                                                                                                                                                 <  1:626746/150‑1 (MQ=255)
        aaTTTCGATGGGAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag‑cca                                                                                                                                                <  2:263700/150‑1 (MQ=255)
        aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag‑cca                                                                                                                                                <  2:721008/150‑1 (MQ=255)
                 ggAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACct                                                                                                                                       <  2:942698/150‑1 (MQ=255)
                               aaaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTc                                                                                                                         <  1:926179/150‑1 (MQ=255)
                                  aCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTTTCCTCCTTTTTCCG‑CCAGGGGAAACTGCTGGGGGCGGTTCtttt                                                                                                                     >  2:865970/1‑151 (MQ=255)
                                     gcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTAttttt                                                                                                                  <  2:116135/151‑1 (MQ=255)
                                         cacCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTTTTACCATTTTCCCG‑CCCGGGGAACCCGCTTGGGCCGGTTCTAtttttttcc                                                                                                              >  2:892969/1‑147 (MQ=255)
                                          acCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGTCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCg                                                                                                             <  1:339522/151‑1 (MQ=255)
                                            cATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGGATTTG‑TCCTCCTTTTTCAGCCCAGGAAAACCTCTTGGGGCGGTTCCTtttttttccgtt                                                                                                            >  1:229758/1‑145 (MQ=255)
                                            cATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTGTCCTCCTTTTTCAG‑CCACGGGAAACTGCTGGGGGCGGTTCTtttttt                                                                                                                  >  1:236736/1‑144 (MQ=255)
                                             aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGCTCCTCATTTTCCCG‑CCCGGGGAACCTGCTGGGGGCGGTTTTTTTTTTTGCCGttt                                                                                                          >  1:369461/1‑151 (MQ=255)
                                             aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGttt                                                                                                          <  2:716319/151‑1 (MQ=255)
                                               tAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCTGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTa                                                                                                         <  2:385384/150‑1 (MQ=255)
                                                aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTTTCATCATTTTTCCC‑CCCGGGGAACCCGCTTGGGGCGGTTTTTTTTTTTTCCGTTTTcc                                                                                                       >  2:659540/1‑151 (MQ=255)
                                                 aaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATCATTTTTCAG‑CCAGGGGAAACCGCTGGGGGCGGTTTTATTTTTTTCCGttt                                                                                                          >  2:14858/1‑147 (MQ=255)
                                                                 aCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTg                                                                                      <  1:647323/151‑1 (MQ=255)
                                                                  cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTGGTTCTCCTTTTCCCG‑CCAGGGGAAACCGCTGGGGGCGGTTCTTTTTTTTTCCGCTTTCCTTCttttttttttt                                                                                         >  1:47450/1‑146 (MQ=255)
                                                                  cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCt                                                                                       <  1:129706/149‑1 (MQ=255)
                                                                      ttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTTCCCG‑CCAGGGGAAACCGCTGGGGGCGGTTCTATTTTTTTCCGGTTTCCTGCTTTTCttttt                                                                                          >  2:713886/1‑142 (MQ=255)
                                                                            cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTg                                                                           <  1:46017/151‑1 (MQ=255)
                                                                               tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTTTTCTCCTTTTCCCG‑CCCGGGAAAACCGCTCGGGGCGGTTCTATTTTTTCCCGTTTTCCTGCCTTTTTTTTTTCTTCGCCccggccgg                                                                          >  2:570641/1‑149 (MQ=255)
                                                                                cGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCTTCATTTTTCAG‑CCAGGGGAACCTGCTGGGGGCGGTTTTATTTTTTTCCGCTTTCCTGCTTATTTTTTTTCTTCGCGTCGGCTGGGcc                                                                       >  2:758338/1‑151 (MQ=255)
                                                                                      gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATCCTTTTTCAG‑CCCGGGAAAACTCCTGGGGGCGGTTCTTTTTTTTTCCCTTTTCCTCCTTATTTTTTTTCTTCGCGCCGGCCGGGCCCCCgcc                                                                 >  2:254693/1‑149 (MQ=255)
                                                                                        ttGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGCATTTTGTCCTCATTTTTCAG‑CCAGGGGAAACTTCTGGGGGCGGGTCTTTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTTCGCGTCGGCCgggacccc                                                                    >  1:734663/1‑140 (MQ=255)
                                                                                         tGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCgg                                                               <  2:244462/150‑1 (MQ=255)
                                                                                          ggCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGtt                                                             <  1:484428/151‑1 (MQ=255)
                                                                                              atgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTc                                                         <  1:82474/151‑1 (MQ=255)
                                                                                                gcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCCTTATTTTTCAC‑CCCGGGGAACCCCCTGGGGGCGGTTTTCTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTTCGCGTCGGCTGGGCCCCCGCCGGTTAtt                                                          >  1:331237/1‑148 (MQ=255)
                                                                                                 cgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTTTCATCATTTTCCAG‑CCAGGGGAAACTGCTGGGGCCGGTTCTATTTTTTTCCGCTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCCCGGGTATTTAGt                                                      >  2:58560/1‑151 (MQ=255)
                                                                                                 cgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGt                                                      <  2:510425/151‑1 (MQ=255)
                                                                                                   atTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTTTTCCTCATTTTCCAG‑CCCGGGGAAACCGCTGGGGGCGATTCTTTTTTTTGCCGTTTTTCTGCTTTTTTTTTTTCTTCGCGGCGGGTGGGCCCCCCGCGGGTATTTAGt                                                      >  2:475406/1‑149 (MQ=255)
                                                                                                   atTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCg                                                    <  2:333267/151‑1 (MQ=255)
                                                                                                       tCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg                                                 <  2:684845/150‑1 (MQ=255)
                                                                                                       tCTGGTTTACTTTGGTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCAGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg                                                 <  1:681319/150‑1 (MQ=255)
                                                                                                        cTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCCTCATTTTCCCG‑CCCGGGGAAACTTCTGGGGGCGATTCTTTTTTTTTCCGTTTTCCTGCTTATCTTTTTTCTTCGCGTCGGCCGGGCCCCCCTCGGTTTTTTAGTCGGCGTg                                               >  2:191663/1‑151 (MQ=255)
                                                                                                          ggTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGt                                                <  1:567686/148‑1 (MQ=255)
                                                                                                               aCTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATTATTTTCCAC‑CCAGGGGAAACTGCTGGGGGCGGGTTTTTTTTTTTCCGCTTTTCTGCTTTTTTTTTTTCTGCGCGTCGGCTGGGCCCACGCCGGGTTTTTTGTCGGCGTGAAGGcgg                                        >  1:150338/1‑150 (MQ=255)
                                                                                                                   tgttgCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcg                                     <  2:421165/150‑1 (MQ=255)
                                                                                                                       gCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATaa                                 <  1:87840/150‑1 (MQ=255)
                                                                                                                        ccGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAg                               <  1:808554/151‑1 (MQ=255)
                                                                                                                          gATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTg                             <  2:479704/151‑1 (MQ=255)
                                                                                                                                ggC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgaga                        <  1:892940/150‑1 (MQ=255)
                                                                                                                                          ttGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  1:583960/1‑149 (MQ=255)
                                                                                                                                           tGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  1:561169/1‑148 (MQ=255)
                                                                                                                                                       gccag‑ccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGt  <  2:962582/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285703‑3286002

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