Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606895,8560GC36.0% 34.7 / 19.5 25A26P (GCC→CCC) ECB_00836putative Regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (14/2);  new base C (0/9);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 2.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGC  >  REL606/895781‑895927
                                                                            |                                                                         
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                           <  1:4301117/76‑1 (MQ=255)
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                           <  1:570242/76‑1 (MQ=255)
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                           <  1:1118428/76‑1 (MQ=255)
 tGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcgc                                                                          <  2:6878055/76‑1 (MQ=255)
  gACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgcc                                                                         <  1:2444987/76‑1 (MQ=255)
   aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgcca                                                                        <  1:2618797/76‑1 (MQ=255)
      cgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccagcc                                                                     <  1:4192337/76‑1 (MQ=255)
       gTGCGCTGGCGGTAGGGGTGGCCCGGGCACTCGTTGCGGTGGGCTGGCAGTGGTGGCGGCTTAACAGCgccagcca                                                                    <  2:2798849/76‑1 (MQ=255)
       gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCGGTCGTGGCGGCTTAACAGCgccagcca                                                                    <  1:1525177/76‑1 (MQ=255)
       gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccagcca                                                                    <  1:7053990/76‑1 (MQ=255)
       gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccagcca                                                                    <  1:4692034/76‑1 (MQ=255)
          cgcTGGCAGTAGTGGTGGCGCTGGCCCTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCacac                                                                 <  1:2239360/76‑1 (MQ=255)
          cgcTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCacac                                                                 <  1:7120021/76‑1 (MQ=255)
            cTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGGGGGGTGAAAGCGCCAGACCCAccc                                                               >  2:2281928/1‑75 (MQ=255)
                    ggtggggGCGCGGGCCCTGGTGGGGTGGGGCGGGGGGTGGGGGCGGCTAACCAGCCCCAGCCACACCATCGAAACg                                                       <  1:6991037/75‑1 (MQ=255)
                    agtGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGGTTACAACGCCCGCCCCACCCTCGCAAcc                                                       >  1:5498246/1‑75 (MQ=255)
                      tggtggCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGGTTAAAACCCCCGCCCCACCCTCGAAAAGCa                                                     >  1:4663695/1‑76 (MQ=255)
                       gtgggCGCGGGCCCTGGTTGCGCGGGGCTGCCAGTGGTGGCGGCTAACCACCCCCAGCCACACCATCGAAACGCAg                                                    <  2:4329536/73‑1 (MQ=255)
                       ggttgCGCGGGCCCTCGTTGCGCGGGGCGGGCAGTCGGGGCGGCTTAACGGCCCCAGCCACACCATCGAAACGCAg                                                    <  1:7814245/76‑1 (MQ=255)
                        ggggCGCGGCCCCTGGTGGCGCGGGGGGGGCGGTGGGGGCGGCTAACCGGCCCCAGCCACACCATCGAAACGCAgc                                                   <  1:4635135/74‑1 (MQ=255)
                          ggCGCGGGCCCTCGTGGCGGGGGGCGGGGGGTGGGGGCGTCTAAACGCCCCCAGCCACACCATCGAAACGCAgcgc                                                 <  1:6190093/76‑1 (MQ=255)
                             ggtGGCCCTGGTTGCGGGGGGCTGGGGGGGGGGGCGGTTAAACAGCCCCAGCCACACCATCGAAACGCAGcgtgcg                                              <  1:1606230/74‑1 (MQ=255)
                             gcTGGCCCTAGTTGGGCGGGGCGGACAGTCGTGGGCGATAAACAGCCCCAGCCACACCATCGAAACGCAGcgcgcg                                              <  1:1777875/76‑1 (MQ=255)
                             gcTGGCCCGGGTTGCGGGGGGCTGGCTGTCGTGGCGCTTAAACGCCCCCAGCCACACCATCGAAACGCAGcgcgcg                                              <  1:321635/76‑1 (MQ=255)
                             gcTGCCACTCGTCGCGCGGGGGGGGCGGTGGGGGCGGCTTACCCCCCCCAGCCACACCATCGAAACGCAGcgcgcg                                              <  1:5264043/76‑1 (MQ=255)
                              cTGGCACTGGTTGGGCGGGGCGGGCGGGGGGGGCGGCTAACCCCCCCCAGCCACACCATCGAAACGCAGCGCgcgg                                             <  2:743918/76‑1 (MQ=255)
                               tGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAGCGGCAGCCACACCATCCAAACGCAGCGCgcggc                                           >  1:7756555/1‑76 (MQ=255)
                                 gCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAGCGCCCGCCACACCATCGAAACGCAGAGCGCGgcgc                                        >  1:1454038/1‑76 (MQ=255)
                                  cACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGGTTACCAGGCCCGCCCCACCCTCGCAACGCAGCGCGCGgggcg                                       >  2:6412490/1‑75 (MQ=255)
                                  cACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAACGCCCGCCACACCATCCAAACGCCGCGCCCGGGGCt                                       >  1:2720577/1‑76 (MQ=255)
                                   aCTCGTTGCGCTGGGCTGGCAGTCGTGGGGGGCTTACACCGCCAGCCCCACCACCGAAACGCAGCGCGCGgcgccg                                      >  1:6718542/1‑74 (MQ=255)
                                    cTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTTAAAGCGCCAGCCACAACCTCCAAAAGCAGCGCGCGGCGCTgg                                     >  2:4658897/1‑75 (MQ=255)
                                     tCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTTACAGCGCCAGCCCCCCCATCGAAACGCAGCGCGCGGGGCTgca                                    >  1:2018252/1‑74 (MQ=255)
                                     tCGGTGCGGGGGGGTGGCAGGGGGGGCGGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGaa                                    <  1:3707732/76‑1 (MQ=255)
                                      cGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTTAAAAGGCCCACCACCACCACCCAAAGCAGCGCGCGGCGCTGgaa                                   >  2:6264791/1‑76 (MQ=255)
                                       gTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGGGCTGgaaa                                  >  1:3423658/1‑76 (MQ=255)
                                        ttGCGCTGGGCTGGCAGTCGTGGCGGCTTTAAAACGCCCACCACACCCTTCAAAAGCACCGGGCGGGGCCGgaaaa                                 >  1:1502060/1‑75 (MQ=255)
                                        gtGCGCTGGGCGGGCAGGGGGGGCGGCTAAACACCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAg                                 <  1:5970545/75‑1 (MQ=255)
                                         tGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCCGCCACACCCTCGAAACCCAGGGCGCGGCGCCGCAAAgg                                >  1:3249693/1‑75 (MQ=255)
                                         tGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGGAAAGt                                >  1:5921732/1‑76 (MQ=255)
                                         tGCGCTGGGCTGGCAGTCGTGGCGGCTTAACACCGCCCGCCACACCAACTAAACGCAGCGCGCGGCCGCGGAAAAg                                 >  1:6231680/1‑76 (MQ=255)
                                          gcgcTGGGCTGGCAGTCGTGGCGGGTTAACAGCCCCCGCCCCACCACCGCAACGCAGCGCGCGGCGCGGACAAgaa                               >  1:121938/1‑74 (MQ=255)
                                             ctgggctggCAGTCGTGGCGGGTTTAAAACGCCAGCCCCCCCCCCCAAAAGCAGCGCGCCGGGCCGCAAAGaaaaa                            >  2:5070884/1‑75 (MQ=255)
                                              tgggctggCAGTCGTGGCGGGTTTACAAAGCCAGCCCCACCCCAGAAACGCAGCGCGCGGGGCCGGTAAGTAAAgc                           >  2:5778958/1‑76 (MQ=255)
                                              tgggctggCAGTCGTGGCGGGCTTACAGGGCCCGCCACCCCAACGACAAGCCACGCGCGGCGCCGgaaaagaaaac                            >  2:5501881/1‑75 (MQ=255)
                                              ggggctggCTGTGGGGGGGGCTAACCCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgc                           <  1:2919305/75‑1 (MQ=255)
                                               gggctggCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCCTCCAAAAGCAGCGCGCGGCGCGGGAAAGGAAAacg                          >  1:3019022/1‑76 (MQ=255)
                                               gggcgggcAGTGGGGGCGGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgcg                          <  1:6961813/76‑1 (MQ=255)
                                                 gctggCGGTGGGGGCGGCTTAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCa                        <  2:7366380/76‑1 (MQ=255)
                                                  cgggCAGGCGGGGCGGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAg                       <  1:5734875/74‑1 (MQ=255)
                                                   gggCGGTCGGGGCGCCTTACCCGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAgg                      <  1:1598272/75‑1 (MQ=255)
                                                     gCAGTCGTGGCGGGTTAAAAACCCCAGCCCCACCCACGGAAAGCAGCGCGCGGCGCCGGAAAATAAAAAGCCGGaa                    >  2:5042339/1‑76 (MQ=255)
                                                      cAGTCGTGGCGGGTTTAAAGCGCCCGCCACACCACCCGAAAGCAGCGCGCCGCGCCGGAAAAAAAAACGCAGGaaa                   >  2:41124/1‑75 (MQ=255)
                                                      cAGTCGTGGCGGCTTAACAGGGCCCCCCCCACCAACGCAAAGCAGCGCGCGGGCCTGTAAAGTAAAACGCCGGAAc                   >  2:7977004/1‑76 (MQ=255)
                                                      cAGTCGTGGCGGCTTAACAGCGCCAGCCCCACCACCGAAACGCAGCGCGCGGCGCTGAAAAGAACAAAGCAGGaga                   >  2:4302842/1‑74 (MQ=255)
                                                       aGTCGTGGCGGGTTTAAAGCGCCCGCCCCACCCCCGAAACGCAGCGCGCGGCGCTGAAAACAAAAGAGCAGGaaat                  >  2:2031722/1‑74 (MQ=255)
                                                        gTCGTTGCGGCTTTTTAGCGCCAGCCACACCATCGAAAAGCAGCGCGCGGCGCTGGAAAGTAAAACGCAGGGAccg                 >  2:173394/1‑74 (MQ=255)
                                                        gTCGTGGCGGGTTTACAGCGCCCGCCCCCCCATCCAAAAGCAGCGCGCGGCGCTGAAAAGAAAAACGCAGGAAcag                 >  2:259980/1‑74 (MQ=255)
                                                        gTCGTGGCGGGTTTACAACGCCAGCCCCACCATCGAAACACAGCGCGCGGGGCCGGAAAGGAAAAAGCAGGAAcag                 >  2:4054601/1‑74 (MQ=255)
                                                        gTCGTGGCGGGTTAACAGAGCCAGCCCCACCCTCCAAACGCAGCGCGCGGCGCTGGAAAGTAAAGCGGAGGGACTg                 >  2:3369061/1‑76 (MQ=255)
                                                          cGTGGCGGGTTAACCGCGCCAGCCACACCACCGAAAAGCAGCGCGCGGCGCTGAAAAGCAAAGAGCAGGAACTggc               >  2:2894511/1‑74 (MQ=255)
                                                          cGTGGCGGCTTTACAGCGCCCGCCCCACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGAc               >  2:7354989/1‑76 (MQ=255)
                                                          cGTGGCGGCTTAACAGCGCCCGCCACACCATCGAAACGCAGCGCGCGGCGCTGGAAAGGAAAACGCAGGAAATGAc               >  2:3922755/1‑76 (MQ=255)
                                                           ggggcggcTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACg              <  1:3961676/74‑1 (MQ=255)
                                                           gTGGCGGGTTTACAGCGCCCGCCACACCATCGAAAACCAGCGCGCCGCGGTGGAAAAAAAAACGCAGGGACTGACg              >  2:3458548/1‑76 (MQ=255)
                                                            tGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACga             >  1:907277/1‑76 (MQ=255)
                                                            tGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACga             >  2:2017357/1‑76 (MQ=255)
                                                            tGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCCGCGCTGAAAAGTAAAGCGCAGGAACTGACga             >  2:1979081/1‑76 (MQ=255)
                                                             ggcggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa            >  2:3929213/1‑76 (MQ=255)
                                                              gcggcTAACCACCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaag           <  1:7066826/76‑1 (MQ=255)
                                                               cggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaaga          >  1:3873238/1‑76 (MQ=255)
                                                                  ccTTACCACCGCAAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaa        <  1:474720/75‑1 (MQ=255)
                                                                   ttaaaCAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa       <  2:5348196/73‑1 (MQ=255)
                                                                   cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAACAGAAAAGCGCAGGAACTGACGAAGaaaa       >  2:1939545/1‑76 (MQ=255)
                                                                   cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa       >  1:7613698/1‑76 (MQ=255)
                                                                    ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa      >  1:75942/1‑76 (MQ=255)
                                                                    ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa      >  2:4743263/1‑76 (MQ=255)
                                                                    ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa      >  2:6639940/1‑76 (MQ=255)
                                                                     tAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc     >  1:1591238/1‑76 (MQ=255)
                                                                      aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca    >  2:2713753/1‑76 (MQ=255)
                                                                      aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca    >  1:1943529/1‑76 (MQ=255)
                                                                      aaCACCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca    <  2:2851941/76‑1 (MQ=255)
                                                                      aaCACCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca    <  2:2918345/76‑1 (MQ=255)
                                                                        gcccgcCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagc  <  1:3183742/73‑1 (MQ=255)
                                                                            |                                                                         
ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGC  >  REL606/895781‑895927

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