Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,416,8380TG29.1% 95.1 / 13.7 55intergenic (+1/‑4)ECB_01334/ECB_01335putative Lysozyme/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/16);  new base G (0/16);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 3.44e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCG  >  REL606/1416764‑1416913
                                                                          |                                                                           
gCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAtt                                                                            <  2:288306/76‑1 (MQ=255)
 ccTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTc                                                                           <  2:986721/76‑1 (MQ=255)
  cTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCa                                                                          <  2:319343/76‑1 (MQ=255)
  cTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCa                                                                          <  1:1372085/76‑1 (MQ=255)
  cTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCa                                                                          <  1:547811/76‑1 (MQ=255)
     ttGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATg                                                                       <  2:1286510/76‑1 (MQ=255)
       gACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGAc                                                                     >  1:554768/1‑76 (MQ=255)
       gACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGAc                                                                     >  1:648049/1‑76 (MQ=255)
          gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCa                                                                  <  2:1446614/76‑1 (MQ=255)
          gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCa                                                                  <  1:1542034/76‑1 (MQ=255)
           tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa                                                                 >  2:1077438/1‑76 (MQ=255)
           tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa                                                                 <  1:224471/76‑1 (MQ=255)
           tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa                                                                 >  2:1346497/1‑76 (MQ=255)
           tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa                                                                 >  1:111259/1‑76 (MQ=255)
            aTCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaaa                                                                >  1:1176048/1‑76 (MQ=255)
            aTCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCaaa                                                                <  1:340313/76‑1 (MQ=255)
                gACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTa                                                            >  2:1387664/1‑76 (MQ=255)
                  cAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATg                                                          >  1:1181135/1‑76 (MQ=255)
                      tCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAAca                                                      <  2:307169/76‑2 (MQ=255)
                        tttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACCATGAACTGa                                                    <  2:1411393/76‑1 (MQ=255)
                          tATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAtt                                                  <  2:574185/76‑1 (MQ=255)
                          tATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGAtt                                                  <  1:283256/76‑1 (MQ=255)
                            tttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgt                                                <  2:1109433/76‑1 (MQ=255)
                            tttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgt                                                <  1:121214/76‑1 (MQ=255)
                             ttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgtg                                               <  2:1249633/76‑1 (MQ=255)
                               cAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGaa                                             >  2:1053101/1‑76 (MQ=255)
                                aaGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGaaa                                            >  1:1267939/1‑76 (MQ=255)
                                aaGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGaaa                                            <  1:1243481/76‑1 (MQ=255)
                                   tttAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAggg                                         >  2:99006/2‑76 (MQ=255)
                                    ttAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGa                                        <  2:1583777/76‑1 (MQ=255)
                                     tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt                                       <  1:1490520/76‑1 (MQ=255)
                                     tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt                                       <  2:669597/76‑1 (MQ=255)
                                      aaaGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATc                                      <  2:1569557/76‑1 (MQ=255)
                                      aaaGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATc                                      >  2:1278005/1‑76 (MQ=255)
                                       aagaagATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATcc                                     >  2:378489/1‑76 (MQ=255)
                                          aagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa                                  >  1:742748/1‑76 (MQ=255)
                                             aTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt                               >  2:868983/1‑76 (MQ=255)
                                              tAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAttt                              >  1:1349275/1‑76 (MQ=255)
                                                aGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTg                            <  1:1007940/76‑1 (MQ=255)
                                                aGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTg                            <  1:209980/76‑1 (MQ=255)
                                                  cGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAGAGGGGTCCGCAATTTTgag                          <  2:125507/76‑1 (MQ=255)
                                                  cGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgag                          >  1:1425284/1‑76 (MQ=255)
                                                      aGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATa                      <  1:341129/76‑1 (MQ=255)
                                                       ggCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaa                     >  1:1255993/1‑76 (MQ=255)
                                                       ggCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaa                     <  1:390229/76‑1 (MQ=255)
                                                          tGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGt                  >  1:533337/1‑76 (MQ=255)
                                                          tGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGt                  <  1:1011845/76‑1 (MQ=255)
                                                           gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt                 >  2:1028516/1‑76 (MQ=255)
                                                           gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt                 >  1:426790/1‑76 (MQ=255)
                                                            cACGGGGAGAGTAATTGAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa                <  2:1319828/76‑1 (MQ=255)
                                                            cACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa                <  2:34864/76‑1 (MQ=255)
                                                              cGGGGAGAGTAATTCAATGACTCAAAACTATCAACTGATTGTGGAAGGGATCCGCAATTTTGAGAATAAAGTTACg              <  1:882738/76‑1 (MQ=255)
                                                              cGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACg              <  1:882736/76‑1 (MQ=255)
                                                               ggggAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACgg             <  1:641620/76‑1 (MQ=255)
                                                                  gagaGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTaa          >  1:585979/1‑76 (MQ=255)
                                                                     agTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACtt       >  1:967151/1‑76 (MQ=255)
                                                                          ttCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCg  <  1:619844/76‑1 (MQ=255)
                                                                          |                                                                           
GCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCG  >  REL606/1416764‑1416913

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: ‑