Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606716,8320TG22.1% 171.4 / 15.1 68A31A (GCA→GCC) ECB_00663hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/36);  new base G (15/0);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.28e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAT  >  REL606/716758‑716907
                                                                          |                                                                           
aaaaGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTg                                                                            <  2:893130/76‑1 (MQ=255)
 aaaGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGc                                                                           <  1:831789/76‑1 (MQ=255)
    gAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCAtt                                                                        <  1:1806723/76‑1 (MQ=255)
    gAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCAtt                                                                        <  2:535453/76‑1 (MQ=255)
      aTGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTAc                                                                      <  2:1944624/76‑1 (MQ=255)
      aTGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTAc                                                                      <  2:540294/76‑1 (MQ=255)
       tGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGAtttac                                                                     >  2:1704849/1‑73 (MQ=255)
        gAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCa                                                                    <  1:166539/76‑1 (MQ=255)
        gAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCa                                                                    <  2:874600/76‑1 (MQ=255)
        gAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCa                                                                    <  1:344047/76‑1 (MQ=255)
            aTTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTAACAACgg                                                                >  2:1243139/1‑75 (MQ=255)
             ttAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTg                                                               >  1:1881383/1‑76 (MQ=255)
              tAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTTCCAACGTGa                                                              >  2:1051870/1‑76 (MQ=255)
               aGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGaa                                                             <  1:1628285/76‑1 (MQ=255)
               aGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGATTACCATCGTgga                                                             >  2:215366/1‑74 (MQ=255)
                 aTACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGAATACCAACCTGaaaa                                                           >  1:1149570/1‑74 (MQ=255)
                  tACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACaa                                                          <  2:1637078/76‑1 (MQ=255)
                  tACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGCTTACCATCGTGAaaaa                                                          >  1:1410965/1‑76 (MQ=255)
                   aCGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAAt                                                         <  2:1368662/76‑1 (MQ=255)
                   aCGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAAt                                                         <  1:1024950/76‑1 (MQ=255)
                      cgATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAgg                                                      <  1:1995054/76‑1 (MQ=255)
                       gaggcggcgGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGa                                                     <  1:1060300/73‑1 (MQ=255)
                       gATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGa                                                     <  2:1449111/76‑1 (MQ=255)
                       gATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGa                                                     <  2:777417/76‑1 (MQ=255)
                        aTGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGAc                                                    <  1:276488/76‑1 (MQ=255)
                        aTGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGAc                                                    <  1:1817418/76‑1 (MQ=255)
                            ggcggcAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCg                                                >  2:60398/1‑76 (MQ=255)
                              cggcAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTa                                              <  2:1957171/76‑1 (MQ=255)
                                gcAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTAct                                            >  1:623586/1‑76 (MQ=255)
                                 cAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTActg                                           >  1:825823/1‑76 (MQ=255)
                                  aGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTActgc                                          <  1:391928/76‑1 (MQ=255)
                                   ggTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCGTTACCATCGTGAACAATAGGACGCCGTActgct                                         >  1:518767/1‑76 (MQ=255)
                                   ggTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTActgct                                         <  1:1278597/76‑1 (MQ=255)
                                     ttACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATTGGACGCCGCACTGCTgg                                       >  1:302867/1‑76 (MQ=255)
                                     ttACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAAAAATAGGACGCCGTAATGCTgg                                       >  1:1698819/1‑76 (MQ=255)
                                     ttACACGGTTATCACTACGATAGCCAATACTGGCTGGGGGCTTACCAATGTGAAAAAAAGGACACCCGGCTGCTgg                                       >  1:803859/1‑76 (MQ=255)
                                     tgacacGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTgg                                       <  2:1102509/74‑1 (MQ=255)
                                      tACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGa                                      <  1:18845/76‑1 (MQ=255)
                                      tACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGa                                      <  1:451988/76‑1 (MQ=255)
                                        cacGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTTCCCTCGTGGACAAAAGGGGGCCGCACTGCTGGgaa                                    >  2:1471339/1‑76 (MQ=255)
                                          cGGTTATCACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAAtt                                  >  2:619424/1‑76 (MQ=255)
                                              tATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCtt                              >  1:1934259/1‑76 (MQ=255)
                                                tCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAAtttttttc                            >  2:1541305/1‑74 (MQ=255)
                                                tCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTAc                            <  2:928526/76‑1 (MQ=255)
                                                 cACTACGATAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGGACTGCTGGAAATTTCTTTcc                           >  1:109444/1‑76 (MQ=255)
                                                  aCTACGATAGCCAATACTGGCTGGGGGATTACCATCGTGAAAAATAGGACGCCGGACTGCTGGAAATTTCTTTCCt                          >  2:328573/1‑76 (MQ=255)
                                                   cTACGATAGCCAATACTGGCTGGTGGATTACCATCGTGAACAATAGGAAGCCGTACTGCTGGAAATTTCTTTCCtt                         >  2:1319415/1‑76 (MQ=255)
                                                   cTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTTTTACCtt                         >  1:242002/1‑76 (MQ=255)
                                                   cTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCtt                         >  1:1859697/1‑76 (MQ=255)
                                                     aCGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTa                       >  1:1551688/1‑76 (MQ=255)
                                                     aCGATAGCCAATACTGGCTGGTGCATTACCATCGTCAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTa                       <  1:448341/76‑1 (MQ=255)
                                                      cGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTAt                      <  1:1973828/76‑1 (MQ=255)
                                                       gATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATg                     <  1:788429/76‑1 (MQ=255)
                                                         tAGCCAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGAc                   >  2:612445/1‑76 (MQ=255)
                                                          aGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACt                  <  1:985176/76‑1 (MQ=255)
                                                          aGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACt                  <  1:375360/76‑1 (MQ=255)
                                                             cAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTccc               >  2:934548/1‑76 (MQ=255)
                                                             cAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTCTTACCTTTATTACTccc               >  1:974364/1‑76 (MQ=255)
                                                             cAATACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTccc               >  2:87591/1‑76 (MQ=255)
                                                               aTACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCcg             <  1:888227/76‑1 (MQ=255)
                                                               aTACTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTTTTTCCTTTTTGACTCCCcg             >  2:711181/1‑76 (MQ=255)
                                                                tACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTCTTACCTTTATGACTcccccc            >  2:171987/1‑74 (MQ=255)
                                                                tACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCcgc            <  2:1684447/76‑1 (MQ=255)
                                                                 aCTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCcgcg           <  1:269288/76‑1 (MQ=255)
                                                                 aCTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCcgcc           >  2:12288/1‑75 (MQ=255)
                                                                 aCTGGCTGGGGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGGAATTTTTTTCCTTTTTGACTCCCcccg           >  2:849579/1‑76 (MQ=255)
                                                                  ctggctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCgcgc          >  2:176222/1‑76 (MQ=255)
                                                                  ctggctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCgcgc          >  1:1610/1‑76 (MQ=255)
                                                                    ggctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGcgccg        >  2:1634742/1‑76 (MQ=255)
                                                                    ggctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGcgccg        >  2:87139/1‑76 (MQ=255)
                                                                      ctggTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCgccgcc      <  1:796615/76‑1 (MQ=255)
                                                                         gTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTa   >  2:571511/1‑76 (MQ=255)
                                                                         gTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTa   <  2:1808539/76‑1 (MQ=255)
                                                                          tGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAt  <  2:1763844/76‑1 (MQ=255)
                                                                          |                                                                           
AAAAGAATGAACATTAGATACGCGATGCGGCGGCAGGTTACACGGTTATCACTACGATAGCCAATACTGGCTGGTGCATTACCATCGTGAACAATAGGACGCCGTACTGCTGGAAATTTCTTACCTTTATGACTCCCCGCGCCGCCGTAT  >  REL606/716758‑716907

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