Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,285,8530GT44.0% 26.6 / 13.9 25A126S (GCC→TCC) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base G (5/9);  new base T (11/0);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTT  >  NC_000913/3285704‑3286003
                                                                                                                                                      |                                                                                                                                                      
ggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgc                                                                                                                                                       <  2:222056/151‑1 (MQ=255)
    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagc                                                                                                                                                   <  1:654770/151‑1 (MQ=255)
    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccag                                                                                                                                                    <  2:301756/150‑1 (MQ=255)
     ttaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagc                                                                                                                                                   <  2:218634/150‑1 (MQ=255)
       aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTgccagcca                                                                                                                                                 <  2:1007735/150‑1 (MQ=255)
       aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGccagccag                                                                                                                                                <  1:601560/151‑1 (MQ=255)
        aTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAgg                                                                                                                                               <  1:460738/151‑1 (MQ=255)
          ttCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCATCTTTTTCCAGCCCggg                                                                                                                                              >  1:572658/1‑150 (MQ=255)
                  aagaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAGTTTTGTCTTCCTTTgtcacaccggggaa                                                                                                                                           >  2:235771/1‑135 (MQ=255)
                           gagaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTTTCCTCATTTTTCACCCCGGGGAAACCGCCggggg                                                                                                                                >  2:242465/1‑146 (MQ=255)
                                aaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTGCTTATCATTTTCCCGCCCGGGGAAACTGCTGGGGGCGGtt                                                                                                                           >  2:122778/1‑147 (MQ=255)
                                  ccGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttt                                                                                                                     <  1:747818/151‑1 (MQ=255)
                                     cgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTCTCCTCATTTTTCCGCCCGGGGAAACTGCTGGGGGCGGTTCTttttttt                                                                                                                  >  1:538596/1‑150 (MQ=255)
                                                             ggCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTTTCCTTATTTTTCCGCCCGGGGAAACTCCTGGGGGGGGTTTTTTTTTTTTCCGTTTTCCTCCtttttttttt                                                                                           >  1:852333/1‑150 (MQ=255)
                                                                       ggCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTc                                                                                 <  1:494864/150‑1 (MQ=255)
                                                                          ccAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCt                                                                             <  1:926713/151‑1 (MQ=255)
                                                                                gCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTGTCCTTCTTTTTCCGCCCGGGGAAACCGCTGGGGGCGGGTTTTTTTTTTTCCGTTTTTCTGCTTTTTTTTTTTCTTcgcggccggcgggccc                                                                       >  2:368426/1‑139 (MQ=255)
                                                                                     gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTTTCCTCATTTTTCCGCCCGGGGAAACTGCTGGGGGCGGTTCTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTTTTCGCGGCGGCTGGGccccc                                                                     >  1:882472/1‑146 (MQ=255)
                                                                                           cgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTCTCATCATTTTTCAACCAGGGGAACCCGCTGGGGCCGGTTCTTTTTTTTTCCGCTTTCCTTCTTTTTTTTTTTCTTCGCGGCCGCTGGACCCCCCTCgggttt                                                            >  1:229817/1‑147 (MQ=255)
                                                                                                 gatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTTGTCTTCATTTGCCACCCAGGGGAAACCGCTCGGGCGGGTTCTTTTTTTTTCCGGTTACCTTCTTATTTTTTTTCTTCCGGTCGGGTTGGCCCCCGCCgggtttt                                                           >  1:968460/1‑141 (MQ=255)
                                                                                                   ttttCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTTCTCATCCTTTTGCCGCCAGGGAAACCCGCTTGGGGCGGTTCTTTTTTTTTCCGGTTTCCTGCTTTTTTTTTTTCTGCGCGGCTGGTGGGCCCCCGTCGGGTTTTCAGTAgg                                                    >  2:932901/1‑151 (MQ=255)
                                                                                                    tttCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTTATCATTTTCCAACCAGGGGAAACTGCTGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCTCGGTtttttcgt                                                       >  1:284587/1‑143 (MQ=255)
                                                                                                          gTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAg                                             <  2:1020376/151‑1 (MQ=255)
                                                                                                           tttACTTTGTTGCCGATTATGGCGGGAATTTTCTCCTCATTTTCCCGCCAGGGGAAACTGCTGGGGGCGATTCTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTGCGCGGCGGCTGGGGCCCCGTCGGTTTTTTTGTCGGCGTGAAgg                                            >  1:56235/1‑151 (MQ=255)
                                                                                                              aCTTTGTTGCCGATTATGGCGGGAATTTGCTCCTCAATTGCCAGCCAGGGGAAACCGCTGGGGGCGCTTCTTTTTTTTTCCGTTTACCTGCTTTTCTTTTTTCTGCGCGTCGGCTGGGCCCCCGCCGGGTATTTAGTCGGCggggaagcg                                          >  1:296225/1‑142 (MQ=255)
                                                                                                                             aTGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCga                           <  1:687132/150‑1 (MQ=255)
                                                                                                                                 cGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGaa                      <  1:923297/151‑1 (MQ=255)
                                                                                                                                          tGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAg               <  2:95816/150‑1 (MQ=255)
                                                                                                                                                tcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTg        >  2:266419/1‑151 (MQ=255)
                                                                                                                                                tcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTg        >  2:911251/1‑151 (MQ=255)
                                                                                                                                                  atTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGcc      >  2:909249/1‑151 (MQ=255)
                                                                                                                                                      gccagccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGtt  <  1:301355/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTT  >  NC_000913/3285704‑3286003

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