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

GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  NC_000913/3285703‑3286002
                                                                                                                                                       |                                                                                                                                                      
gggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTg                                                                                                                                                        <  2:703824/151‑1 (MQ=255)
       taATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagcc                                                                                                                                                 <  1:626746/150‑1 (MQ=255)
        aaTTTCGATGGGAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagcca                                                                                                                                                <  2:263700/150‑1 (MQ=255)
        aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagcca                                                                                                                                                <  2:721008/150‑1 (MQ=255)
                 ggAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACct                                                                                                                                       <  2:942698/150‑1 (MQ=255)
                               aaaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTc                                                                                                                         <  1:926179/150‑1 (MQ=255)
                                  aCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTTTCCTCCTTTTTCCGCCAGGGGAAACTGCTGGGGGCGGTTCtttt                                                                                                                     >  2:865970/1‑151 (MQ=255)
                                     gcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttttt                                                                                                                  <  2:116135/151‑1 (MQ=255)
                                         cacCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTTTTTACCATTTTCCCGCCCGGGGAACCCGCTTGGGCCGGTTCTAtttttttcc                                                                                                              >  2:892969/1‑147 (MQ=255)
                                          acCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGTCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCg                                                                                                             <  1:339522/151‑1 (MQ=255)
                                            cATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGGATTTGTCCTCCTTTTTCAGCCCAGGAAAACCTCTTGGGGCGGTTCCTtttttttccgtt                                                                                                            >  1:229758/1‑145 (MQ=255)
                                            cATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTGTCCTCCTTTTTCAGCCACGGGAAACTGCTGGGGGCGGTTCTtttttt                                                                                                                  >  1:236736/1‑144 (MQ=255)
                                             aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCCTCATTTTCCCGCCCGGGGAACCTGCTGGGGGCGGTTTTTTTTTTTGCCGttt                                                                                                          >  1:369461/1‑151 (MQ=255)
                                             aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGttt                                                                                                          <  2:716319/151‑1 (MQ=255)
                                               tAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCTGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTa                                                                                                         <  2:385384/150‑1 (MQ=255)
                                                aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTTTCATCATTTTTCCCCCCGGGGAACCCGCTTGGGGCGGTTTTTTTTTTTTCCGTTTTcc                                                                                                       >  2:659540/1‑151 (MQ=255)
                                                 aaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCATCATTTTTCAGCCAGGGGAAACCGCTGGGGGCGGTTTTATTTTTTTCCGttt                                                                                                          >  2:14858/1‑147 (MQ=255)
                                                                 aCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTg                                                                                      <  1:647323/151‑1 (MQ=255)
                                                                  cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTGGTTCTCCTTTTCCCGCCAGGGGAAACCGCTGGGGGCGGTTCTTTTTTTTTCCGCTTTCCTTCttttttttttt                                                                                         >  1:47450/1‑146 (MQ=255)
                                                                  cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCt                                                                                       <  1:129706/149‑1 (MQ=255)
                                                                      ttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTTCCCGCCAGGGGAAACCGCTGGGGGCGGTTCTATTTTTTTCCGGTTTCCTGCTTTTCttttt                                                                                          >  2:713886/1‑142 (MQ=255)
                                                                            cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTg                                                                           <  1:46017/151‑1 (MQ=255)
                                                                               tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTTTTTCTCCTTTTCCCGCCCGGGAAAACCGCTCGGGGCGGTTCTATTTTTTCCCGTTTTCCTGCCTTTTTTTTTTCTTCGCCccggccgg                                                                          >  2:570641/1‑149 (MQ=255)
                                                                                cGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCTTCATTTTTCAGCCAGGGGAACCTGCTGGGGGCGGTTTTATTTTTTTCCGCTTTCCTGCTTATTTTTTTTCTTCGCGTCGGCTGGGcc                                                                       >  2:758338/1‑151 (MQ=255)
                                                                                      gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCATCCTTTTTCAGCCCGGGAAAACTCCTGGGGGCGGTTCTTTTTTTTTCCCTTTTCCTCCTTATTTTTTTTCTTCGCGCCGGCCGGGCCCCCgcc                                                                 >  2:254693/1‑149 (MQ=255)
                                                                                        ttGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGCATTTTGTCCTCATTTTTCAGCCAGGGGAAACTTCTGGGGGCGGGTCTTTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTTCGCGTCGGCCgggacccc                                                                    >  1:734663/1‑140 (MQ=255)
                                                                                         tGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCgg                                                               <  2:244462/150‑1 (MQ=255)
                                                                                          ggCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGtt                                                             <  1:484428/151‑1 (MQ=255)
                                                                                              atgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTc                                                         <  1:82474/151‑1 (MQ=255)
                                                                                                gcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCCTTATTTTTCACCCCGGGGAACCCCCTGGGGGCGGTTTTCTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTTCGCGTCGGCTGGGCCCCCGCCGGTTAtt                                                          >  1:331237/1‑148 (MQ=255)
                                                                                                 cgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTTTCATCATTTTCCAGCCAGGGGAAACTGCTGGGGCCGGTTCTATTTTTTTCCGCTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCCCGGGTATTTAGt                                                      >  2:58560/1‑151 (MQ=255)
                                                                                                 cgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGt                                                      <  2:510425/151‑1 (MQ=255)
                                                                                                   atTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTTTCCTCATTTTCCAGCCCGGGGAAACCGCTGGGGGCGATTCTTTTTTTTGCCGTTTTTCTGCTTTTTTTTTTTCTTCGCGGCGGGTGGGCCCCCCGCGGGTATTTAGt                                                      >  2:475406/1‑149 (MQ=255)
                                                                                                   atTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCg                                                    <  2:333267/151‑1 (MQ=255)
                                                                                                       tCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg                                                 <  2:684845/150‑1 (MQ=255)
                                                                                                       tCTGGTTTACTTTGGTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCAGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg                                                 <  1:681319/150‑1 (MQ=255)
                                                                                                        cTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCCTCATTTTCCCGCCCGGGGAAACTTCTGGGGGCGATTCTTTTTTTTTCCGTTTTCCTGCTTATCTTTTTTCTTCGCGTCGGCCGGGCCCCCCTCGGTTTTTTAGTCGGCGTg                                               >  2:191663/1‑151 (MQ=255)
                                                                                                          ggTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGt                                                <  1:567686/148‑1 (MQ=255)
                                                                                                               aCTTTGTTGCCGATTATGGCGGGAATTTTCTCATTATTTTCCACCCAGGGGAAACTGCTGGGGGCGGGTTTTTTTTTTTCCGCTTTTCTGCTTTTTTTTTTTCTGCGCGTCGGCTGGGCCCACGCCGGGTTTTTTGTCGGCGTGAAGGcgg                                        >  1:150338/1‑150 (MQ=255)
                                                                                                                   tgttgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcg                                     <  2:421165/150‑1 (MQ=255)
                                                                                                                       gCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATaa                                 <  1:87840/150‑1 (MQ=255)
                                                                                                                        ccGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAg                               <  1:808554/151‑1 (MQ=255)
                                                                                                                          gATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTg                             <  2:479704/151‑1 (MQ=255)
                                                                                                                                ggCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgaga                        <  1:892940/150‑1 (MQ=255)
                                                                                                                                          ttGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  1:583960/1‑149 (MQ=255)
                                                                                                                                           tGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                >  1:561169/1‑148 (MQ=255)
                                                                                                                                                       gccagccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGt  <  2:962582/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT  >  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: