Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,950,0381.→C 3.9% 0.6 ‑3.1 0.051coding (1424/1662 nt)tarmethyl‑accepting chemotaxis protein II
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (40/9):  total (40/11)
Fisher's exact test strand distribution p-value = 4.31e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGCGTTCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGCCAATGCGACTTGATCGATGCCACGGCTCTGT  >  REL606/1950004‑1950073
                                   |                                   
tGCGTTCTGTTGCGTGACGCGATCCATTTCCGAAA.c                                    <  1:7519719/36‑1 (MQ=255)
 gCGTTCTGTTGCGTGACGCGATCCATTTCCGAAA.cc                                   >  1:3050043/1‑36 (MQ=255)
   gTTCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGc                                 >  1:6799697/1‑36 (MQ=255)
   gTTCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGc                                 >  1:3929522/1‑36 (MQ=255)
    ttCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGcc                                <  1:4523764/36‑1 (MQ=255)
    ttCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGcc                                <  1:7746935/36‑1 (MQ=255)
     tCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGCCa                               >  1:6591873/1‑36 (MQ=255)
        gTTGCGTGACGCGATCCATTTCCGAAA.CCGCCAATg                            >  1:3763419/1‑36 (MQ=255)
         ttGCGTGACGCGATCCATTTCCGAAA.CCGCCAATGc                           >  1:3545167/1‑36 (MQ=255)
          tGCGTGACGCGATCCATTTCCGAAA.CCGCCAATGCg                          >  1:3240984/1‑36 (MQ=255)
          tGCGTGACGCGATCCATTTCCGAAA.CCGCCAATGCg                          >  1:7688469/1‑36 (MQ=255)
           gCGTGACGCGATCCATTTCCGAAA.CCGCCAATGCGa                         >  1:3148459/1‑36 (MQ=255)
           gCGCGACGCGA.CCAATTCCAAAACCCGCCAATGCGa                         <  1:3638578/36‑1 (MQ=2)
            gttGACGCGATCCATTCCATAAC.CCGCCAATGCGAc                        <  1:3931553/34‑1 (MQ=12)
            cGTGACGCGATCCATTTCCGAAA.CCGCCAATGCGAc                        >  1:4205953/1‑36 (MQ=255)
            cGTGACGCGATCCATTTCCGAAA.CCGCCAATGCGAc                        >  1:424083/1‑36 (MQ=255)
            cGTGACGCGATCCATTTCCGAAA.CCGCCAATGCGAc                        >  1:7923445/1‑36 (MQ=255)
             gTGACGCGATCCATTTCCGAAA.CCGCCAATGCGagt                       >  1:1470596/1‑34 (MQ=255)
             gTGACGCGATCCATTTCCGAAA.CCGCCAATGCGACt                       >  1:3899148/1‑36 (MQ=255)
             gTGACGCGATCCATTTCCGAAA.CCGCCAATGCGACt                       >  1:6646247/1‑36 (MQ=255)
               gccgcgATCCATTCCCAAAA.CCGCCAATGCGACTTg                     <  1:1567619/34‑1 (MQ=25)
               gACGCGATCCATTTCCGAAA.CCGCCAATGCGACTTg                     >  1:2651962/1‑36 (MQ=255)
                 cgcgATCCATTTCCGAAA.CCGCCAATGCGACTTGAt                   >  1:4223885/1‑36 (MQ=255)
                  gcgATCCATTTCCGAAA.CCGCCAATGCGACTTGATc                  >  1:6574187/1‑36 (MQ=255)
                  gcgATCCATTTCCGAAA.CCGCCAATGCGACTTGATc                  >  1:2449003/1‑36 (MQ=255)
                  gcgATCCATTTCCGAAA.CCGCCAATGCGACTTGATc                  >  1:2414850/1‑36 (MQ=255)
                  gcgATCCATTTCCGAAA.CCGCCAAGGCGACTTGATc                  >  1:2728438/1‑36 (MQ=255)
                  gcgATCCATTAACAAAC.CCGCCAATGCGACTTGATc                  <  1:5917081/36‑1 (MQ=16)
                  gaaaTCAATTACCAAAA.CCGCCAATGCGACTTGATc                  <  1:7201441/33‑1 (MQ=16)
                   cgATCCATTTCCGAAA.CCGCCAATGCGACTTGATCg                 >  1:4633571/1‑36 (MQ=255)
                    gATCCATTTCCGAAA.CCGCCAATGCGACTTGATCGa                >  1:4104024/1‑36 (MQ=255)
                    gATCCATTTCCGAAA.CCGCCAATGCGACTTGATCGa                >  1:833388/1‑36 (MQ=255)
                    gATCCATTTCCGAAA.CCGCCAATGCGACTTGATCGa                >  1:3683562/1‑36 (MQ=255)
                    gATCCATTTCCGAAA.CCGCCAATGCGACTTGATCGa                >  1:305122/1‑36 (MQ=255)
                     aTCCATTTCCGAAA.CCGCCAATGCGACTTGATCGAt               >  1:2223949/1‑36 (MQ=255)
                     aTCCATTTCCGAAA.CCGCCAATGCGACTTGATCGAt               >  1:1756319/1‑36 (MQ=255)
                      tCCATTCCCAAAA.CCGCCAATGCGACTTGATCGATg              <  1:496985/36‑1 (MQ=37)
                      tCCATTCCCAAAA.CCGCCAATGCGACTTGATCGATg              <  1:7070692/36‑1 (MQ=37)
                       ccATTTCCGAAA.CCGCCAATGCGACTTGATCGATGc             >  1:7381180/1‑36 (MQ=255)
                       ccATTTCCGAAA.CCGCCAATGCGACTTGATCGATGc             >  1:123699/1‑36 (MQ=255)
                        cctttCCAAAACCCGCCAATGCGACTTGATCGATGc             <  1:402091/34‑1 (MQ=21)
                        cATTTCCGAAA.CCGCCAATGCGACTTGATCGATGcc            >  1:3400998/1‑36 (MQ=255)
                        cATTTCCGAAA.CCGCCAATGCGACTTGATCGATGcc            >  1:7535001/1‑36 (MQ=255)
                         aTTTCCGAAA.CCGCCAATGCGACTTGATCGATGCCa           >  1:2861858/1‑36 (MQ=255)
                         aTTTCCGAAA.CCGCCAATGCGACTTGATCGATGCCa           >  1:2624280/1‑36 (MQ=255)
                           ttCCGAAA.CCGCCAATGCGACTTGATCGATGCCACg         >  1:271325/1‑36 (MQ=255)
                             ccaaaa.CCGCCAATGCGACTTGATCGATGCCACGGc       <  1:6175730/33‑1 (MQ=255)
                             ccGAAA.CCGCCAATGCGACTTGATCGATGCCACGGc       <  1:7837919/36‑1 (MQ=255)
                             ccGAAA.CCGCCAATGCGACTTGATCGATGCCACGGc       >  1:4124293/1‑36 (MQ=255)
                              cGAAA.CCGCCAATGCGACTTGATCGATGCCACGGct      >  1:6349501/1‑36 (MQ=255)
                              cGAAA.CCGCCAATGCGACTTGATCGACGCCACGGct      >  1:1184025/1‑36 (MQ=255)
                               gAAA.CCGCCAATGCGACTTGATCGATGCCACGGctc     >  1:6574283/1‑36 (MQ=255)
                               gAAA.CCGCCAATGCGACTTGATCGATGCCACGGctc     >  1:7744949/1‑36 (MQ=255)
                               gAAA.CCGCCAATGCGACTTGATCGATGCCACGGctc     >  1:2397177/1‑36 (MQ=255)
                               aaaa.CCGCCAATGCGACTTGATCGATGCCACGGctc     <  1:403232/35‑1 (MQ=255)
                                aaa.CCGCCAATGCGACTTGATCGATGCCACGGctct    >  1:6049987/1‑36 (MQ=255)
                                aaa.CCGCCAATGCGACTTGATCGATGCCACGGctct    >  1:5393895/1‑36 (MQ=255)
                                aaa.CCGCCAATGCGACTTGATCGATGCCACGGctct    >  1:975886/1‑36 (MQ=255)
                                 aa.CCGCCAATGCGACTTGATCGATGACACGGCTCTg   >  1:978247/1‑36 (MQ=255)
                                  a.CCGCCAATGCGACTTGATCGATGCCACGGCTCTGt  >  1:4901159/1‑36 (MQ=255)
                                   |                                   
TGCGTTCTGTTGCGTGACGCGATCCATTTCCGAAA.CCGCCAATGCGACTTGATCGATGCCACGGCTCTGT  >  REL606/1950004‑1950073

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