Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913766,7160GT40.0% 7.5 / 10.6 25A245S (GCG→TCG) mngA2‑O‑alpha‑mannosyl‑D‑glycerate specific PTS enzyme II
Reads supporting (aligned to +/- strand):  ref base G (9/6);  new base T (10/0);  total (19/6)
Fisher's exact test for biased strand distribution p-value = 5.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGTGCT  >  NC_000913/766474‑766966
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                
tGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAATGCAAACCCCTGGATTTTTTCGGCTgg                                                                                                                                                                                                                                                           >  1:168121/1‑249 (MQ=255)
                                                     aaCGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGGTTAACGGCGGAACAACTCCAAAATGCAAACCCCTGGTTTTTTTCCGCTGGAGTGGCCCCCAAAGAGAGTGAGGGTTTTTATgggcttccccccgctttactgac                                                                                                                                                                                                    >  1:389865/1‑230 (MQ=255)
                                                        cgccaatctcccGTCTGCGGAAGGGCCTGGCGGGGGGGTGTTGCTCACTCTTATGTGGCGACGGAATATCTGAAAAAAGCCGGACGAAGCCGCGGTTTAAAGTTTGAGTTGAAAAAACAAGGCGGAAGCGTTTTGGAAGGGTGTTTATGGGGGGAGCAATTAAATGGTGCAGCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGTAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGt                                                                                                                                                                                                <  1:620071/240‑1 (MQ=255)
                                                                            aaCGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGAACAACTCAACAATGGAACCGCCTGTATTTTTTCGGCTGAAAGAGCCATCAAGGAGAGTGGGGGGTTTAATTGCATTTCCCCGCCTTTCTGGCCTGTTGCCCGGCCGAATCGCCAtt                                                                                                                                                                           >  1:381602/1‑250 (MQ=255)
                                                                                           gggTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTTCGGCTGGAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCCGAAGCGTTGGTcc                                                                                                                                                           >  1:210800/1‑251 (MQ=255)
                                                                                                 tGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTACCGGCATTGAAGGGCGTTTAACGGCGGATTAAATCAATAGTGCCAACCCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGCGAGAGCGTTTTAATGGAATTTCCGCCCTTTCAGTGACTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCcaacc                                                                                                                                                       >  1:660975/1‑248 (MQ=255)
                                                                                                                                 aaaaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGAACAAATCAATAAAGCAACCGCCCGTATTTTTTCGGGTGAAGCCGCCACCCAGGGGAGGGAGGGGTTTAAAGGGATTCCCGCGCTTTTAGTGCGTGTTTTCGCACCGATTACGCATGCAGGAGCGGTTGTTCCACCAACGCTTACCCTcccgacgaacgaggtgacccgc                                                                                                                    >  2:547835/1‑230 (MQ=255)
                                                                                                                                   aaaGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGCTCAACTCAACAGTCCCACCGCCTGTATTTTTTCGGCTGAAGTCGCCACCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCCTGCCAAAGCGTTGATCCCACAAGCGCTTACCCTCAAGCGTAGCGGTGAGACGcgtac                                                                                                                  >  1:343774/1‑251 (MQ=255)
                                                                                                                                       ccGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACATGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCCACCGCCTGTATTTTTGCGGCTGAAGGCGCCATCAAGGAGAGTGAGCGGTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATTCAACAAGCGCTTACCCTCCAGCGTAGCGATGAGACGcgtaccgt                                                                                                              >  1:67245/1‑250 (MQ=255)
                                                                                                                                         ggACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATCGTGCAACCGCCTGTATTTTTTCGGCTGCAGGCGCCCTCAAGGAGAGTGAGCGGTTTTATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCCATTCGCCATGCCGGAGCGTTGATCCAACAAGCGCTTCCCCCCCAGCGGAGCGGTGAGACGcgtaccgtac                                                                                                            >  1:304539/1‑250 (MQ=255)
                                                                                                                                                 aCTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCCACCGCCTGTATTTTTTCGGCTGAAGTCGCCATCAAGGGGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGGCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGAt                                                                                                    >  1:143615/1‑250 (MQ=255)
                                                                                                                                                          aaaTGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGGGCAACCGCCTGTATTTTTGCGGCTGAAGGCGCCATCAAGGGGAGTGAGCGTTTTTATGGCATTCCCGCGCTTTCCGTGCCTGTTGCCGAGCCGATTCGCCCTGCAGAAGCGTTGATCCAACAAGCGCTTACCCCCAAGCGGAGCGGTGAGACGCGTACCGTAAAGCCAGATAcagatcggg                                                                                          >  1:10536/1‑243 (MQ=255)
                                                                                                                                                                            aCAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCCACTCAATAGTGCAACCGCCCGTATTTTTTCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTTATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCCTACCCTCAAGCGGAGCGGTGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACgg                                                                        >  1:323110/1‑251 (MQ=255)
                                                                                                                                                                                    gcTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCCATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCCTTTCCGCGCTTTCAGTGCCTGTTGCCGCGCCGATTCGCCATGCAGAAGCGTTGGTCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTCCCGTACAGCAAGATACGCAACCGGGGAAAAGTGGCAAAAAGGAGCTGa                                                                  >  1:22737/1‑250 (MQ=255)
                                                                                                                                                                                          ggCATTGAAGGGCGTTTAACGGCGGATCAACTTAATAGTGCAACCGCCTGTATTTTTTCGGCTGAAGGCGCCATCCAGGGGAGGGAGCGGTTTAATTGGCTTCCCCCGCCTTCCGTGCCTGTTGCCGAGCCGGTTCGCCCTGCAGAAGCGGTTGTCCAACAAGCGCCTTACCTCAACCGGAGCGGTTGGACGCGGACCCGTAAGCCAGAGACGCAACCGGGGCAAAGTGTCAAAAAGGAGGTGAAAAAggg                                                           >  2:510124/1‑250 (MQ=255)
                                                                                                                                                                                           gtgtaggggggatttgagcgcaggATCAACTCAATAGCGCACCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCa                                                          <  2:168121/229‑1 (MQ=255)
                                                                                                                                                                                              ttGAAGGGCGTTTAACGGCGGCTCAACTCAAAAGTGCAAACACCTCTATTTTTGCGGCTGAAGACGCCATCCAGGAGAGGGAGCGCTTTTATGGCATTCCCGCGCCTTCAGTGCCTGGTGCCGAGCCGATTCTCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGCTGAGACGCGTACCGTACAGCAAGAAACGCAAACGGTGAAAAGTGTCAAAAAGGAGCTGAAACAGGCACTg                                                       >  1:414153/1‑251 (MQ=255)
                                                                                                                                                                                                              gggcgaaaaaccaaaaaGTGCAATCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATGCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCGGAAACAGGCACTGTTGAGCGGAATctct                                        <  1:570395/235‑1 (MQ=255)
                                                                                                                                                                                                                        cTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATTAAGGAGAGTGAGCGTATTAATGGCATTCCCGCGATTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGAGCAAAGTGTCAACACGGCGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTAcc                             >  1:81902/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTAGTGCAACAGCCTGTATTTTTGCGGCTGAAGTCGACATCAAGGCGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGCGCCGATTCGCCATGCAGACGCGTTGATCCAACAAGCGCTTACCCACAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACGGTTg                         >  2:662099/1‑251 (MQ=255)
                                                                                                                                                                                                                              aGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGCGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGCGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGAt                       >  1:613523/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCCGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGCGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGCATCCCTTTTGCCGTACCGTTGATTGTcgg                 >  2:298389/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aTTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCa         <  1:2620/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aTTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCAc        <  1:674359/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttGCGGCTGAAGTCGCCATCAAGGAGAGTGTGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCACCCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGt     <  2:381602/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTACTGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGTGCt  >  1:39650/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAAc                                                               >  2:361407/1‑190 (MQ=255)
                                                                                                                                                                                                                                                     gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTggtggcggcagcacagaa    >  1:38457/1‑232 (MQ=255)
                                                                                                                                                                                                                                                     gCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGAGGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAAc                                                               <  1:361407/190‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                
TGCTCTCGGCGCTGGATGACAAGGGAGGCACGCAACCTTCTGCCTCTTTTTCCAACGCGCCAACTATCGTCTGCGTAACGGCCTGTCCGGCGGGTATTGCTCACACCTATATGGCTGCGGAATATCTGGAAAAAGCCGGACGCAAACTCGGCGTAAATGTTTACGTTGAAAAACAAGGCGCTAACGGCATTGAAGGGCGTTTAACGGCGGATCAACTCAATAGTGCAACCGCCTGTATTTTTGCGGCTGAAGTCGCCATCAAGGAGAGTGAGCGTTTTAATGGCATTCCCGCGCTTTCAGTGCCTGTTGCCGAGCCGATTCGCCATGCAGAAGCGTTGATCCAACAAGCGCTTACCCTCAAGCGTAGCGATGAGACGCGTACCGTACAGCAAGATACGCAACCGGTGAAAAGTGTCAAAACGGAGCTGAAACAGGCACTGTTGAGCGGAATCTCTTTTGCCGTACCGTTGATTGTCGCGGGGGGCACGGTGCT  >  NC_000913/766474‑766966

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