Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,887,4650G→T 5.4% 0.8 ‑0.5 ‑0.155intergenic (‑121/‑38)yobG/ECB_01797hypothetical protein/hypothetical protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/1):  ref base (25/27):  total (27/28)
Fisher's exact test strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.88e-01

TCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTG  >  REL606/1887431‑1887500
                                  |                                   
tCGGAATGAGTATCTTGCTTCAGCGCATATTGCCGc                                    <  1:4690149/36‑1 (MQ=255)
tCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGc                                    <  1:133686/36‑1 (MQ=255)
tCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGc                                    <  1:3555087/36‑1 (MQ=255)
  ggAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTa                                  >  1:2375611/1‑36 (MQ=255)
  ggAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTa                                  <  1:2146516/36‑1 (MQ=255)
   gcaTGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc                                 <  1:2819168/34‑1 (MQ=255)
   gAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc                                 <  1:2895358/36‑1 (MQ=255)
   gAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc                                 >  1:957271/1‑36 (MQ=255)
   gAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTAc                                 <  1:4890631/36‑1 (MQ=255)
    aaTGAGTAGCTTGCTTCAGCGCATATTGCCGCTACg                                <  1:5581277/36‑1 (MQ=255)
    aaTGAGTAGCTTGCTTCAGCGCATATTGCCGCTACg                                >  1:909462/1‑36 (MQ=255)
     aTGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGa                               >  1:2455569/1‑36 (MQ=255)
     aTGAGGAGCTTGCTTCAGCGCATATTGCCGCTACGa                               >  1:556885/1‑36 (MQ=255)
      tGAGTAGCTTGCTTCAGCGCATATTGCCTCTACGAt                              >  1:3639372/1‑36 (MQ=255)
      tGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAt                              <  1:964170/36‑1 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             <  1:3873355/36‑1 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             <  1:3254669/36‑1 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             <  1:4368491/36‑1 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             <  1:3117778/36‑1 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             <  1:964796/36‑1 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             >  1:2314780/1‑36 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             >  1:1567971/1‑36 (MQ=255)
       gAGTAGCTTGCTTCAGCGCATATTGCCGCTACGAtt                             <  1:1763877/36‑1 (MQ=255)
         gTAGCTTGCTTCAGCGCATATTGCCTCTACGATTaa                           <  1:1474219/36‑1 (MQ=255)
         gTAGCTTGCTTCAGCGCATATTGCCGCTACGATTaa                           <  1:4139745/36‑1 (MQ=255)
         gTAGCTTGCTTCAGCGCATATTGCCGCTACGATTaa                           <  1:3451106/36‑1 (MQ=255)
         gTAGCTTGCTTCAGCGCATATTGCCGCTACGATTaa                           >  1:2688278/1‑36 (MQ=255)
          tAGCTTGCTTCAGCGCATATTGCCTCTACGATTAAg                          >  1:934551/1‑36 (MQ=255)
           aGCTTGCTTCAGCGCATATTTCCGCTACGATTAAGc                         <  1:4636858/36‑1 (MQ=255)
            gcttgcttCAGCGCATATTGCCGCTACGATTAAGCg                        >  1:5465423/1‑36 (MQ=255)
            gcttgcttCAGCGCATATTGCCGCTACGATTAAGCg                        >  1:4254077/1‑36 (MQ=255)
            gcttgcttCAGCGCATATTGCCGCTACGATTAAGCg                        >  1:4946310/1‑36 (MQ=255)
             cttgcttCAGCGCATATTGCCGCTACGCTTAAGCGa                       >  1:1423096/1‑36 (MQ=255)
              ttgcttCAGCGCATATTGCCGCTACGATTAAGCGaa                      >  1:1281292/1‑36 (MQ=255)
               tgcttCAGCGCATATTGCCGCTACGATTAAGCGAAc                     <  1:2234931/36‑1 (MQ=255)
                  ttCAGCGCATATTGCCGCTACGATTAAGCGAACATa                  <  1:3442359/36‑1 (MQ=255)
                  ttCAGCGCATATGGCCGCTACGATTAAGCGAACATa                  <  1:5022686/36‑1 (MQ=255)
                   tCAGCGCATATTGCCGCTACGATTAAGCGAACATaa                 >  1:4953442/1‑36 (MQ=255)
                    cAGCGCATATTGCCGCTACGATTAAGCGAACATaaa                <  1:526084/36‑1 (MQ=255)
                      gcgcATATTGCCGCTACGATTAAGCGAACATaaaaa              >  1:4703504/1‑36 (MQ=255)
                      gcgaATATTGCCGCTACGATTAAGCGAACATaaaaa              >  1:1088256/1‑36 (MQ=255)
                       cgcATATTGCCGCTACGATTAAGCGAACATAAAaag             >  1:5352913/1‑36 (MQ=255)
                        gcATATTGCCGCTACGATTAAGCGAACATAAAaaga            >  1:4825229/1‑36 (MQ=255)
                        gcATATTGCCGCTACGATTAAGCGAACATAAAaaga            >  1:3804430/1‑36 (MQ=255)
                          atatTGCCGCTACGATTAAGCGAACATAAAaagaga          >  1:4466694/1‑36 (MQ=255)
                          atatTGCCGCTACGATTAAGCGAACATAAAaagaga          >  1:5428939/1‑36 (MQ=255)
                          atatTGCCGCTACGATTAAGCGAACATAAAaagaga          <  1:3948888/36‑1 (MQ=255)
                          atatTGCCGCTACGATTAAGCGAACATAAAaagaga          >  1:2070096/1‑36 (MQ=255)
                           tatTGCCGCTACGATTAAGCGAACATAAAaagagaa         <  1:4438366/36‑1 (MQ=255)
                           tatTGCCGCTACGATTAAGCGAACATAAAaagagaa         >  1:813416/1‑36 (MQ=255)
                           tatTGCCGCTACGATTAAGCGAACATAAAaagagaa         <  1:2206669/36‑1 (MQ=255)
                           tatTGCCGCTACGATTAAGCGAACATAAAaagagaa         >  1:1415270/1‑36 (MQ=255)
                             ttGCCGCTACGATTAAGCGAACATAAAaagagaaga       >  1:5392211/1‑36 (MQ=255)
                             ttGCCGCTACGATTAAGCGAACATAAAaagagaaga       <  1:3575979/36‑1 (MQ=255)
                             ttGCCGCTACGATTAAGCGAACATAAAaagagaaga       <  1:805664/36‑1 (MQ=255)
                                  gCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTg  <  1:2218940/36‑1 (MQ=255)
                                  |                                   
TCGGAATGAGTAGCTTGCTTCAGCGCATATTGCCGCTACGATTAAGCGAACATAAAAAGAGAAGAGGTTG  >  REL606/1887431‑1887500

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG