Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 790,873 C→G 14.2% intergenic (+119/‑582) ECB_00731 → / → ECB_00732 conserved hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606790,8730C→G14.2% 5.1 ‑2.9 ‑0.034intergenic (+119/‑582)ECB_00731/ECB_00732conserved hypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (14/15):  total (19/15)
Fisher's exact test strand distribution p-value = 5.26e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.90e-01

AAACGGCGTACCAGTGGCGGGTGAAGTATAGGGCATCTATTTATTCATCATTTGAAAAAAGAGTTGTCC  >  REL606/790840‑790908
                                 |                                   
aaaCGGCGTGCTGGTGGCGGGTGAAGTATAGGGCAt                                   <  1:3745374/36‑1 (MQ=21)
  aCGGCGTACCAGTGGCGGGTGAAGTATAGGGCATCt                                 <  1:3427630/36‑1 (MQ=255)
  aCGGCGTACCAGTGGCGGGGGGAGTATAGGGGATCt                                 >  1:901579/1‑36 (MQ=25)
   cGGCGTACCAGTGGCGGGTGAAGTATAGGGCATCta                                <  1:4209813/36‑1 (MQ=255)
   cGGCGTACCAGTGGCGGGTGAAGTATAGGGCATCta                                <  1:2982057/36‑1 (MQ=255)
   cGGCGTACCAGTGGCGGGTGAAGTATAGGGCATCta                                <  1:7064532/36‑1 (MQ=255)
      cGTACCAGTGGCGGGGGGAGGATAGGGGATCtattt                             >  1:420216/1‑36 (MQ=19)
        tACCAGTGGCGGGTGAAGTATAGGGCATCtatttat                           <  1:6716589/36‑1 (MQ=255)
        tACCAGTGGCGGGGGGGGAAAAAGGGATCtattttt                           >  1:2310489/1‑34 (MQ=0)
        tACCAGTGGCGGGGGGAGGATAGGGGCTCtattttt                           >  1:1010836/1‑34 (MQ=255)
         aCCAGTGGCGGGTGAAGTATAGGGCATCtatttatt                          <  1:2385938/36‑1 (MQ=255)
          ccAGTGGCGGGTGAAGTATAGGGCATCTATTTATtc                         <  1:5425542/36‑1 (MQ=255)
           cAGTGGCGGGTGAAGTATAGGGCATCTATTTATtca                        <  1:2582577/36‑1 (MQ=255)
           cAGTGGCGGGGGGGGTATAGGGGATCtttttttttt                        >  1:5042282/1‑34 (MQ=9)
              tagCGGGTTAAGTATAGAGCATCTATTTATtcatca                     <  1:1255940/34‑1 (MQ=37)
               ggCGGGTGAAGTATAGGGCATCTATTTATTcatcat                    <  1:3573321/36‑1 (MQ=255)
                 cGGGTGAAGTATAGGGCATCTATTTATTCATCAttt                  <  1:1271748/36‑1 (MQ=255)
                 cGGGTGAAGTATAGGGCATCTATTTATTCATCAttt                  <  1:7421499/36‑1 (MQ=255)
                   ggTGAAGTATAGGGCATCTATTTATTCATCATTTGa                >  1:306799/1‑36 (MQ=255)
                    gTGAAGTATAGGGCATCTATTTATTCATCATTTGaa               >  1:6552694/1‑36 (MQ=255)
                     tGAAGTATAGGGCATCTATTTATTCATCATTTGaaa              >  1:4464031/1‑36 (MQ=255)
                       aaGTATAGGGCATCTATTTATTCATCATTTGaaaaa            >  1:2852352/1‑36 (MQ=255)
                       aaGTATAGGGCATCTATTTATTCATCATTTGaaaaa            >  1:153588/1‑36 (MQ=255)
                          tataGGGCATCTATTTATTCATCATTTGAAAAAaga         >  1:7374320/1‑36 (MQ=255)
                          tataGGGCATCTATTTATTCATCATTTGAAAAAaga         >  1:7675243/1‑36 (MQ=255)
                           ataGGGCATCTATTTATTCATCATTTGAAAAAagag        >  1:7322317/1‑36 (MQ=255)
                             aGGGCATCTATTTATTCATCATTTGAAAAAAGAGtt      >  1:516087/1‑36 (MQ=255)
                             aGGGCATCTATTTATTCATCATTTGAAAAAAGAGtt      >  1:1898547/1‑36 (MQ=255)
                              gggCATCTATTTATTCATCATTTGAAAAAAGAGTTg     >  1:3895752/1‑36 (MQ=255)
                              gggCATCTATTTATTCATCATTTGAAAAAAGAGTTg     >  1:3597261/1‑36 (MQ=255)
                              gggCATCTATTTATTCATCATTTGAAAAAAGAGTTg     >  1:3125981/1‑36 (MQ=255)
                               ggCATCTATTTATTCATCATTTGAAAAAAGAGTTGt    <  1:309479/36‑1 (MQ=255)
                                gCATCTATTTATTCATCATTTGAAAAAAGAGTTGTc   >  1:4177552/1‑36 (MQ=255)
                                gCATCTATTTATTCATCACTTGAAAAAAGAGTTGTc   <  1:5504973/36‑1 (MQ=255)
                                 cATCTATTTATTCATCATTTGAAAAAAGAGTTGTcc  >  1:275256/1‑36 (MQ=255)
                                 |                                   
AAACGGCGTACCAGTGGCGGGTGAAGTATAGGGCATCTATTTATTCATCATTTGAAAAAAGAGTTGTCC  >  REL606/790840‑790908

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG