Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606988,1800G→. 4.0% 0.8 0.0 ‑0.150coding (237/747 nt)kdsB3‑deoxy‑manno‑octulosonate cytidylyltransferase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (34/14):  total (36/14)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.17e-01

ACGCGCGCCGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGA  >  REL606/988145‑988215
                                   |                                   
aCGCGCGCCGATCATCAGTCAGGAACAGAACGTCTg                                     >  1:1759419/1‑36 (MQ=255)
 cgcgcgCCGATCATCAGTCAGGAACAGAACGTCTgg                                    <  1:2682451/36‑1 (MQ=255)
  gcgcgcCGATCATCAGTCAGGAACAGAACGTCTGGc                                   >  1:2558471/1‑36 (MQ=255)
   cgcgcCGATCATCAGTCAGGAACAGAACGTCTGGCg                                  >  1:1724914/1‑36 (MQ=255)
      gcCGATCATCAGTCAGGAACAGAACGTCTGGCGGaa                               <  1:423380/36‑1 (MQ=255)
        cGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGt                             <  1:4107024/36‑1 (MQ=255)
        cGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGt                             <  1:5499232/36‑1 (MQ=255)
         gATCATCAGTCAGGAACAGAACGTCTGGCGGAAGtt                            >  1:5861426/1‑36 (MQ=255)
         gATCATCAGTCAGGAACAGAACGTCTGGCGGAAGtt                            <  1:1359642/36‑1 (MQ=255)
          atcatcAGTCAGGAACAGACCGTCTGGCGGAAGTTg                           <  1:3840765/36‑1 (MQ=255)
          atcatcAGTCAGGAACAGAACGTCTGGCGGAAGTTg                           >  1:66887/1‑36 (MQ=255)
          atcatcAGTCAGGAAAAGAACGTATGGCGGAAGTTg                           >  1:3569643/1‑36 (MQ=37)
           tcatcaGTCAGGAACAGAACGTCTGGCGGAAGTTGt                          >  1:2601322/1‑36 (MQ=255)
            catcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTc                         <  1:5210107/36‑1 (MQ=255)
            catcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTc                         <  1:1465179/36‑1 (MQ=255)
            catcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTc                         >  1:1056368/1‑36 (MQ=255)
            catcaGTCAGAAACAGAACGTCTGGCGGAAGTTGTc                         <  1:2970938/36‑1 (MQ=255)
             atcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTCg                        >  1:5427402/1‑36 (MQ=255)
             atcaGTCAGGAACAGAACGTCTGGCGGAAGTTGTCg                        >  1:1574694/1‑36 (MQ=255)
              taagtcagGAACAGAACGTCTGGCGGGAGTTGTCGa                       >  1:4383253/3‑36 (MQ=37)
               cagtcagGAACAGAACGTCTGGCGGAAGTTGTCGaa                      <  1:2929734/36‑1 (MQ=255)
                agtcagGAACAGAACGTCTGGAGGAAGTTGTCGaaa                     >  1:2408547/1‑36 (MQ=255)
                 gtaagGAACAGAACGTCT.GCGGAAGTTGTCGaaaaa                   >  1:3160544/4‑36 (MQ=37)
                   cagGAACAGAACGTCTGGCGGAAGTTGTCGAAAAAt                  >  1:539511/1‑36 (MQ=255)
                   cagAAACAGAACGTCTGGCGGAAGTTGTCGAAAAAt                  <  1:2561545/36‑1 (MQ=255)
                    agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg                 >  1:4305993/1‑36 (MQ=255)
                    agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg                 >  1:3055619/1‑36 (MQ=255)
                    agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg                 >  1:1486589/1‑36 (MQ=255)
                    agGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATg                 >  1:1155624/1‑36 (MQ=255)
                     ggAACAGAACTTCTGGCGGAAGTTGTCGAAAAATgc                <  1:5588551/36‑1 (MQ=255)
                     ggAACAGAACGTCTGGCGGAAGTTGTCGAAAAATgc                >  1:634907/1‑36 (MQ=255)
                      gAACAGAACGTCTGGCGGAAGTTGTCGAAAAATgcg               >  1:4633207/1‑36 (MQ=255)
                      gAACAGAACGTCTGGCGGAAGTTGTCGAAAAATgcg               >  1:832661/1‑36 (MQ=255)
                      gAACAGAACGTCTGGCGGAAGTTGCCGAAAAATggg               >  1:1320810/1‑34 (MQ=37)
                        aCAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCa             >  1:1819073/1‑36 (MQ=255)
                          aGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCCtt           >  1:1044367/1‑36 (MQ=255)
                          aGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCAtt           >  1:3173130/1‑36 (MQ=255)
                          aGAACGTCTGGCGGAAGTTGTCGAAAAATGAGAAtt           >  1:5177819/1‑36 (MQ=37)
                          aGAACGTCTGGAGGAAGTTGTCGAAAAATGCGCAtt           >  1:1963048/1‑36 (MQ=255)
                          aGAACGTCT.GCGGAAGATGTCGACAAATGCGCATTc          >  1:292966/1‑36 (MQ=21)
                           gAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTc          >  1:4481000/1‑36 (MQ=255)
                           gAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTc          >  1:4575860/1‑36 (MQ=255)
                           gAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTc          >  1:1939066/1‑36 (MQ=255)
                            aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCCTTCa         >  1:130747/1‑36 (MQ=255)
                            aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCa         >  1:4832614/1‑36 (MQ=255)
                            aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCa         <  1:4797747/36‑1 (MQ=255)
                            aaCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCa         >  1:4096314/1‑36 (MQ=255)
                             aCGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAg        >  1:3983975/1‑36 (MQ=255)
                              cGTCTGGCGGGAGTTGTCGAAAAATGCGCATTCAGc       <  1:679638/36‑1 (MQ=255)
                               gTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGcg      >  1:4576926/1‑36 (MQ=255)
                               gTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGcg      <  1:1005327/36‑1 (MQ=255)
                                  tGGCGGAAGTTGTCGAAAAATGCGCATTCAGcgacg   >  1:5524317/1‑36 (MQ=255)
                                   ggCGGAAGTTGTCGAAAAATGCGCATTCAGcgacga  >  1:5052528/1‑36 (MQ=255)
                                   |                                   
ACGCGCGCCGATCATCAGTCAGGAACAGAACGTCTGGCGGAAGTTGTCGAAAAATGCGCATTCAGCGACGA  >  REL606/988145‑988215

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG