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

AACTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTG  >  REL606/1882451‑1882601
                                                                           |                                                                           
aaCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGatat                                                                             <  1:431330/76‑1 (MQ=255)
 aCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAtt                                                                            <  2:1026901/76‑1 (MQ=255)
 aCTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAtt                                                                            >  2:887935/1‑76 (MQ=255)
    ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt                                                                         >  2:982113/1‑76 (MQ=255)
    ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt                                                                         <  1:844714/76‑1 (MQ=255)
    ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGt                                                                         >  1:1842773/1‑76 (MQ=255)
    ggACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGt                                                                         <  1:844687/76‑1 (MQ=255)
     gACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTg                                                                        <  1:652933/76‑1 (MQ=255)
      aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg                                                                       >  2:16038/1‑76 (MQ=255)
      aCTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTgg                                                                       >  2:1249446/1‑76 (MQ=255)
       cTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGc                                                                      >  1:1159499/1‑76 (MQ=255)
        tAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCg                                                                     >  2:1849247/1‑76 (MQ=255)
         aaCAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGt                                                                    >  1:1018417/1‑76 (MQ=255)
           cAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCt                                                                  >  2:1225040/1‑76 (MQ=255)
           cAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCt                                                                  <  2:640549/76‑1 (MQ=255)
            aGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCtt                                                                 >  1:972738/1‑76 (MQ=255)
              gTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTcc                                                               >  1:454856/1‑76 (MQ=255)
               tGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCt                                                              >  1:1209880/1‑76 (MQ=255)
                gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTg                                                             >  1:845825/1‑76 (MQ=255)
                gAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTg                                                             >  1:1698250/1‑76 (MQ=255)
                  gCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCg                                                           <  2:416271/76‑1 (MQ=255)
                  gCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCg                                                           <  2:1897194/76‑1 (MQ=255)
                    gAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGtt                                                         >  1:231037/1‑76 (MQ=255)
                     aaaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTc                                                        >  1:687739/1‑76 (MQ=255)
                     aaaTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTc                                                        >  1:1634096/1‑76 (MQ=255)
                       aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa                                                      >  1:1226424/1‑76 (MQ=255)
                       aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTa                                                      <  1:109367/76‑1 (MQ=255)
                       aTGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTa                                                      <  1:1288862/76‑1 (MQ=255)
                        tGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTaa                                                     <  2:1733937/76‑1 (MQ=255)
                         ggTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAAc                                                    <  1:1661421/76‑1 (MQ=255)
                          gTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAAcc                                                   >  1:1802107/1‑76 (MQ=255)
                            tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                 >  1:494038/1‑76 (MQ=255)
                            tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                 >  2:79111/1‑76 (MQ=255)
                            tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                 >  2:919971/1‑76 (MQ=255)
                            tGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACctc                                                 >  1:1982616/1‑76 (MQ=255)
                             gATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCCGCGTTCTAACctct                                                >  2:306254/1‑76 (MQ=255)
                              aTACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCtt                                               <  2:1344398/76‑1 (MQ=255)
                               tACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTc                                              >  2:1058289/1‑76 (MQ=255)
                                aCAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTcc                                             <  1:919971/76‑1 (MQ=255)
                                 cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa                                            >  2:410303/1‑76 (MQ=255)
                                 cAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCa                                            <  1:571291/76‑1 (MQ=255)
                                  aaCTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAg                                           >  2:1932669/1‑76 (MQ=255)
                                    cTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAggg                                         <  2:692493/76‑1 (MQ=255)
                                      cAAACTACCCCCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt                                       <  2:579014/76‑1 (MQ=255)
                                      cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt                                       >  1:1226692/1‑76 (MQ=255)
                                      cAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGtt                                       <  2:517130/76‑1 (MQ=255)
                                       aaaCTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTc                                      >  2:1348238/1‑76 (MQ=255)
                                          cTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATg                                   <  1:1671465/76‑1 (MQ=255)
                                          cTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCCCTTCCACAGTTCATg                                   >  1:1366120/1‑76 (MQ=255)
                                           tACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATgg                                  >  1:504081/1‑76 (MQ=255)
                                            accaccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa                                 >  1:845725/1‑76 (MQ=255)
                                            accaccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGa                                 <  1:651575/76‑1 (MQ=255)
                                              caccGAGAAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAAc                               <  1:1953359/76‑1 (MQ=255)
                                               accGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACt                              <  2:1223515/76‑1 (MQ=255)
                                                ccGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACtt                             <  2:1603705/76‑1 (MQ=255)
                                                 cGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTc                            >  1:309317/1‑76 (MQ=255)
                                                 cGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTc                            <  2:1620280/76‑1 (MQ=255)
                                                  gagaAAAAACTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCg                           <  2:59895/76‑1 (MQ=255)
                                                   agaAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGa                          >  2:1930766/1‑76 (MQ=255)
                                                     aaaaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc                        <  2:1149699/76‑1 (MQ=255)
                                                     aaaaaaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc                        <  2:1761548/76‑1 (MQ=255)
                                                     aaaaaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAAc                        <  1:1356404/76‑1 (MQ=255)
                                                       aaaaCTCACTCAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGt                      <  2:1774056/76‑1 (MQ=255)
                                                         aaCTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGATCATGGAACTTCGAACGTAt                    >  1:1024656/1‑76 (MQ=255)
                                                          actcactcAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATg                   >  2:561673/1‑76 (MQ=255)
                                                          actcactcAAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATg                   <  1:104011/76‑1 (MQ=255)
                                                           ctcactcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGc                  <  1:1275679/76‑1 (MQ=255)
                                                           ctcactcaAAGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGc                  <  1:1379399/76‑1 (MQ=255)
                                                            tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa                 >  1:1241401/1‑76 (MQ=255)
                                                            tcactcaAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCa                 >  1:1752949/1‑76 (MQ=255)
                                                                   aaGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTg          >  1:1563569/1‑76 (MQ=255)
                                                                    aGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTgg         <  2:1358766/76‑1 (MQ=255)
                                                                    aGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTgg         <  1:1052505/76‑1 (MQ=255)
                                                                    aGTGATAGTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTgg         <  2:112055/76‑1 (MQ=255)
                                                                      tGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGt       >  1:163501/1‑76 (MQ=255)
                                                                        atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt     <  2:859984/76‑1 (MQ=255)
                                                                        atatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt     >  2:1905695/1‑76 (MQ=255)
                                                                        atagTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGttt     <  2:360565/76‑1 (MQ=255)
                                                                         tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc    >  1:552928/1‑76 (MQ=255)
                                                                         tatTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTc    <  2:1408659/76‑1 (MQ=255)
                                                                          atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCCGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt   >  2:1295391/1‑76 (MQ=255)
                                                                          atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt   <  1:1722501/76‑1 (MQ=255)
                                                                          atTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCt   >  1:1793106/1‑76 (MQ=255)
                                                                           ttCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTg  >  1:564632/1‑76 (MQ=255)
                                                                           |                                                                           
AACTGGACTAACAGGTGAGCGAAATGGTTGATACAACTCAAACTACCACCGAGAAAAAACTCACTCAAAGTGATATTCGTGGCGTCTTCCTGCGTTCTAACCTCTTCCAGGGTTCATGGAACTTCGAACGTATGCAGGCACTGGGTTTCTG  >  REL606/1882451‑1882601

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