Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,882,5260TG22.8% 175.7 / 10.9 79I21S (ATT→AGT) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (44/17);  new base G (0/18);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.65e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCT  >  REL606/1882455‑1882600
                                                                       |                                                                          
ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt                                                                        >  2:105776/1‑76 (MQ=255)
 gACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTg                                                                       <  1:1819492/76‑1 (MQ=255)
 gACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTg                                                                       <  1:1500463/76‑1 (MQ=255)
  aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg                                                                      <  1:896149/76‑1 (MQ=255)
  aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg                                                                      >  1:1668745/1‑76 (MQ=255)
   cTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGc                                                                     <  1:1207933/76‑1 (MQ=255)
    tAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCg                                                                    >  1:1625997/1‑76 (MQ=255)
     aaCAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATGCGTGGCGt                                                                   >  1:1320020/1‑76 (MQ=255)
       cAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCt                                                                 >  1:1015409/1‑76 (MQ=255)
        aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCtt                                                                >  1:400955/1‑76 (MQ=255)
        aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCtt                                                                <  2:1343877/76‑1 (MQ=255)
         ggTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTc                                                               <  2:520335/76‑1 (MQ=255)
         ggTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTc                                                               <  2:432578/76‑1 (MQ=255)
          gTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTcc                                                              >  1:986400/1‑76 (MQ=255)
           tGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCt                                                             >  1:476333/1‑76 (MQ=255)
            gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTg                                                            <  2:262029/76‑1 (MQ=255)
              gCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCg                                                          >  2:700332/1‑76 (MQ=255)
               cGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGt                                                         >  1:1248171/1‑76 (MQ=255)
               cGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGt                                                         >  2:1504443/1‑76 (MQ=255)
               cGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGt                                                         >  2:1696595/1‑76 (MQ=255)
                gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt                                                        >  1:776961/1‑76 (MQ=255)
                gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt                                                        >  1:748520/1‑76 (MQ=255)
                  aaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCt                                                      <  2:1019868/76‑1 (MQ=255)
                   aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa                                                     <  2:601763/76‑1 (MQ=255)
                   aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa                                                     <  2:2059063/76‑1 (MQ=255)
                    tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa                                                    >  1:1829302/1‑76 (MQ=255)
                    tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTaa                                                    <  1:452905/76‑1 (MQ=255)
                    tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTaa                                                    <  1:1576788/76‑1 (MQ=255)
                       ttGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACct                                                 >  2:515334/1‑76 (MQ=255)
                       ttGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACct                                                 <  1:626121/76‑1 (MQ=255)
                        tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAAcccc                                                >  1:1278530/1‑74 (MQ=255)
                        tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                >  2:1085553/1‑76 (MQ=255)
                        tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                >  1:910600/1‑76 (MQ=255)
                        tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                >  2:638588/1‑76 (MQ=255)
                        tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                >  2:1044198/1‑76 (MQ=255)
                         gATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctct                                               >  1:1223217/1‑76 (MQ=255)
                         gATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctct                                               >  1:1554432/1‑76 (MQ=255)
                          aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt                                              >  1:414495/1‑76 (MQ=255)
                          aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt                                              >  2:1227914/1‑76 (MQ=255)
                          aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCtt                                              <  2:84224/76‑1 (MQ=255)
                           tACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTc                                             <  2:909599/76‑1 (MQ=255)
                            aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc                                            >  2:1862781/1‑76 (MQ=255)
                            aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc                                            >  2:70610/1‑76 (MQ=255)
                            aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc                                            <  2:2009466/76‑1 (MQ=255)
                             cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa                                           <  2:27449/76‑1 (MQ=255)
                             cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa                                           <  2:1389481/76‑1 (MQ=255)
                             cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa                                           <  2:388093/76‑1 (MQ=255)
                              aaCTCAAACTGCCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAg                                          <  2:1887147/76‑1 (MQ=255)
                               aCTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAgg                                         >  2:968032/1‑76 (MQ=255)
                                cTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAggg                                        >  1:1626544/1‑76 (MQ=255)
                                cTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAggg                                        <  2:2016697/76‑1 (MQ=255)
                                  cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt                                      >  2:64147/1‑76 (MQ=255)
                                   aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc                                     >  1:1206061/1‑76 (MQ=255)
                                   aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc                                     >  2:2026451/1‑76 (MQ=255)
                                   aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc                                     <  2:1914727/76‑1 (MQ=255)
                                   aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc                                     <  1:457554/76‑1 (MQ=255)
                                     aCTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCAt                                   <  1:1040578/76‑1 (MQ=255)
                                      cTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATg                                  <  2:1233846/76‑1 (MQ=255)
                                       tACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATgg                                 >  2:49313/1‑76 (MQ=255)
                                         ccaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGaa                               >  2:1079636/1‑76 (MQ=255)
                                          caccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATgggcc                              >  1:1306093/1‑73 (MQ=255)
                                          caccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAAc                              >  2:1721996/1‑76 (MQ=255)
                                            ccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACtt                            <  2:670647/76‑1 (MQ=255)
                                            ccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACtt                            <  2:1110708/76‑1 (MQ=255)
                                                gaAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGaa                        <  1:1833649/76‑1 (MQ=255)
                                                 aaaaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc                       >  2:849737/1‑76 (MQ=255)
                                                   aaaaCTCACTCAAAGTGAGATTCGGGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGt                     <  2:307128/76‑1 (MQ=255)
                                                     aaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTAt                   <  2:234630/76‑1 (MQ=255)
                                                        tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa                >  1:980341/1‑76 (MQ=255)
                                                        tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa                <  2:887697/76‑1 (MQ=255)
                                                        tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa                <  1:1584378/76‑1 (MQ=255)
                                                         cactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAg               >  1:804372/1‑76 (MQ=255)
                                                         cactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATCGAACTTCGAACGTATGCAg               >  2:1352867/1‑76 (MQ=255)
                                                         cactcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAg               <  2:730700/76‑1 (MQ=255)
                                                             caAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCAc           <  2:2086847/76‑1 (MQ=255)
                                                               aaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTg         >  1:360412/1‑76 (MQ=255)
                                                                 gTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTggg       >  2:916460/1‑76 (MQ=255)
                                                                  tGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGt      >  1:1737746/1‑76 (MQ=255)
                                                                   gATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGtt     >  1:1942931/1‑76 (MQ=255)
                                                                    atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt    >  2:1101135/1‑76 (MQ=255)
                                                                    atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt    >  2:128452/1‑76 (MQ=255)
                                                                    atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt    <  2:2212649/76‑1 (MQ=255)
                                                                    atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt    <  2:2212668/76‑1 (MQ=255)
                                                                     tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc   <  2:2161165/76‑1 (MQ=255)
                                                                     tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc   >  2:1637294/1‑76 (MQ=255)
                                                                      atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt  <  2:1657364/76‑1 (MQ=255)
                                                                      atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt  >  1:1150778/1‑76 (MQ=255)
                                                                       |                                                                          
GGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCT  >  REL606/1882455‑1882600

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