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

GGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTT  >  NC_000913/3285704‑3286003
                                                                                                                                                      |                                                                                                                                                      
ggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgc                                                                                                                                                       <  2:222056/151‑1 (MQ=255)
    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccagc                                                                                                                                                   <  1:654770/151‑1 (MQ=255)
    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag                                                                                                                                                    <  2:301756/150‑1 (MQ=255)
     ttaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccagc                                                                                                                                                   <  2:218634/150‑1 (MQ=255)
       aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccagcca                                                                                                                                                 <  2:1007735/150‑1 (MQ=255)
       aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGccagccag                                                                                                                                                <  1:601560/151‑1 (MQ=255)
        aTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAgg                                                                                                                                               <  1:460738/151‑1 (MQ=255)
          ttCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGCTCATCTTTTTCCAGCCCggg                                                                                                                                              >  1:572658/1‑150 (MQ=255)
                  aagaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAGTTTTGTCTTCCTTTgtcacaccggggaa                                                                                                                                           >  2:235771/1‑135 (MQ=255)
                           gagaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTTTCCTCATTTTTCACCCCGGGGAAACCGCCggggg                                                                                                                                >  2:242465/1‑146 (MQ=255)
                                aaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTGCTTATCATTTTCCCGCCCGGGGAAACTGCTGGGGGCGGtt                                                                                                                           >  2:122778/1‑147 (MQ=255)
                                  ccGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttt                                                                                                                     <  1:747818/151‑1 (MQ=255)
                                     cgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTCTCCTCATTTTTCCGCCCGGGGAAACTGCTGGGGGCGGTTCTttttttt                                                                                                                  >  1:538596/1‑150 (MQ=255)
                                                             ggCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTTTCCTTATTTTTCCGCCCGGGGAAACTCCTGGGGGGGGTTTTTTTTTTTTCCGTTTTCCTCCtttttttttt                                                                                           >  1:852333/1‑150 (MQ=255)
                                                                       ggCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTc                                                                                 <  1:494864/150‑1 (MQ=255)
                                                                          ccAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCt                                                                             <  1:926713/151‑1 (MQ=255)
                                                                                gCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTGTCCTTCTTTTTCCGCCCGGGGAAACC‑GCTGGGGGCGGGTTTTTTTTTTTCCGTTTTTCTGCTTTTTTTTTTTCTTcgcggccggcgggccc                                                                       >  2:368426/1‑139 (MQ=255)
                                                                                     gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTTTTCCTCATTTTTCCGCCCGGGGAAACTGCTGGGGGCGGTTCTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTTTTCGCGGCGGCTGGGccccc                                                                     >  1:882472/1‑146 (MQ=255)
                                                                                           cgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTCTCATCATTTTTCAACCAGGGGAACCCGCTGGGGCCGGTTCTTTTTTTTTCCGCTTTCCTTCTTTTTTTTTTTCTTCGCGGCCGCTGGACCCCCCTCgggttt                                                            >  1:229817/1‑147 (MQ=255)
                                                                                                 gatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTGTCTTCATTTGCCACCCAGGGGAAACCGCTCGGGCGGGTTCTTTTTTTTTCCGGTTACCTTCTTATTTTTTTTCTTCCGGTCGGGTTGGCCCCCGCCgggtttt                                                           >  1:968460/1‑141 (MQ=255)
                                                                                                   ttttCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTCTCATCCTTTTGCCGCCAGGGAAACCCGCTTGGGGCGGTTCTTTTTTTTTCCGGTTTCCTGCTTTTTTTTTTTCTGCGCGGCTGGTGGGCCCCCGTCGGGTTTTCAGTAgg                                                    >  2:932901/1‑151 (MQ=255)
                                                                                                    tttCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTTATCATTTTCCAACCAGGGGAAACTGCTGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCTCGGTtttttcgt                                                       >  1:284587/1‑143 (MQ=255)
                                                                                                          gTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAg                                             <  2:1020376/151‑1 (MQ=255)
                                                                                                           tttACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCCTCATTTTCCCGCCAGGGGAAACTGCTGGGGGCGATTCTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTGCGCGGCGGCTGGGGCCCCGTCGGTTTTTTTGTCGGCGTGAAgg                                            >  1:56235/1‑151 (MQ=255)
                                                                                                              aCTTTGTTGCCGATTATGGC‑GGGAATTTGCTCCTCAATTGCCAGCCAGGGGAAACCGCTGGGGGCGCTTCTTTTTTTTTCCGTTTACCTGCTTTTCTTTTTTCTGCGCGTCGGCTGGGCCCCCGCCGGGTATTTAGTCGGCggggaagcg                                          >  1:296225/1‑142 (MQ=255)
                                                                                                                             aTGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCga                           <  1:687132/150‑1 (MQ=255)
                                                                                                                                 c‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGaa                      <  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)
                                                                                                                                                      |                                                                                                                                                      
GGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTT  >  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: ‑