Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,992,578 C→G 16.8% A191P (GCC→CCC)  ECB_03721 ← conserved hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,992,5780C→G16.8% 3.1 ‑1.9 0.060571ECB_03721conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/1):  ref base (30/17):  total (41/19)
Fisher's exact test strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AATATGCCCGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATACGGGACATGAGAAAACTTCGCA  >  REL606/3992543‑3992612
                                   |                                  
aaTATGCCCGCCTAACAGGCCTGTAATGGCAGGGGc                                    >  1:1127687/1‑36 (MQ=255)
aaTATGCCCGCCTAACAGGCCTGTAATGGCAGGGGc                                    <  1:3644172/36‑1 (MQ=255)
 atatGCCCGCCTAACAGGCCTGTAATGGCAGGGGCt                                   <  1:1049320/36‑1 (MQ=255)
 atatGCCCGCCTAACAGGCCTGTAATGGCAGGGGCt                                   >  1:5728468/1‑36 (MQ=255)
  tatGCCCGCCTAACAGGCCTGTAATGGCAGGGGGTg                                  >  1:394686/1‑36 (MQ=255)
  tatGCCCGCCTAACAGGCCTGTAATGGCAGGGGGTg                                  >  1:3671416/1‑36 (MQ=255)
  tatGCCCGCCTAACAGGCCTGTAATGGCAGGGGCTg                                  >  1:3110520/1‑36 (MQ=255)
       ccGCCTAACAGGCCTGTAATGGCAGGGGCTGCGcct                             <  1:6901977/36‑2 (MQ=255)
        cGCCTAACAGGCCTGTAATGGCAGGGGCTTGGCCGt                            >  1:1307258/1‑36 (MQ=37)
        cGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGt                            >  1:1032713/1‑36 (MQ=255)
        cGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGt                            <  1:1396884/36‑1 (MQ=255)
         gCCTAACAGGCCTGTAATGGCAGGGGGTGGGCCggg                           >  1:3907025/1‑34 (MQ=37)
         gCCTAACAGGCCTGTAATGGCAGGGGGTGCGCCGTc                           >  1:6712144/1‑36 (MQ=255)
         gCCTAACAGGCCTGTAATGGCAGGGGATTCGCCGTc                           >  1:7400858/1‑36 (MQ=37)
          ccTAACAGGCCTGTAATGGCAGGGGGTTCGCCGTct                          >  1:7203557/1‑35 (MQ=37)
           cTAACAGGCCTGTCATGGCAGGGGGTGCGCCGTCAt                         >  1:4340849/1‑36 (MQ=38)
           cTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCAt                         <  1:6956557/36‑1 (MQ=255)
           cTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCAt                         <  1:1778347/36‑1 (MQ=255)
            tAACAGGCCTGTAATGGCAGGGGGTGCGCCGTTATa                        >  1:3163418/1‑36 (MQ=37)
            tAACAGGCCTGTAATGGCAGGGGGTGCGCCGTCATa                        <  1:1387531/36‑1 (MQ=255)
            tAACAGGCCTGTAATGGCAGGGGGTGCGCCGGCATa                        >  1:4580792/1‑36 (MQ=37)
            tAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATa                        >  1:6576484/1‑36 (MQ=255)
            tAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATa                        <  1:4729011/36‑1 (MQ=255)
            tAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATa                        <  1:1532555/36‑1 (MQ=255)
               cAGGCCTGTAATGGCAGGGGGTGCGGCGTCATATgg                     >  1:6980184/1‑36 (MQ=37)
               cAGGCCTGTAATGGCAGGGGCTTCGCCGTCCTTCgg                     >  1:4914027/1‑36 (MQ=25)
                aGGCGTGTAATGGCAGGGGCTGCGCCGTCATACggg                    <  1:4346568/36‑1 (MQ=255)
                aGGCCTGTAATGGCAGGGGCTGGGCCGTCAGACcgg                    >  1:7444523/1‑36 (MQ=25)
                aGGCCTGTAATGGCAGGGGCTGCGGCGTCATACggg                    >  1:5982359/1‑36 (MQ=255)
                aGGCCTGTAATGGCAGGGGCTGCGCCGTCATACggg                    >  1:7294062/1‑36 (MQ=255)
                aGGCCTGTAATGGCAGGGGCTGCGCCGTCATACggg                    >  1:4926963/1‑36 (MQ=255)
                 ggCCTGTAATGGCAGGGGCTGCGCCGTCCTACgggg                   >  1:5770968/1‑35 (MQ=255)
                 ggCCTGTAATGGCAGGGGCTGCGCCGTCATACGGGa                   <  1:3604622/36‑1 (MQ=255)
                 ggCCTGTAATGGCAGGGGCGGCGCCGTCCTACcggg                   >  1:7175357/1‑35 (MQ=25)
                   ccTGTAATGGCAGGGGCTGCGCCGTCATACGGGGCa                 >  1:1308534/1‑36 (MQ=255)
                   ccTGTAATGGCAGGGGCTGCGCCGTCATACGGGACa                 <  1:959748/36‑1 (MQ=255)
                    cTGTAATGGCAGGNGCTGCGCCGTCATACGGGACAt                >  1:3247952/1‑36 (MQ=255)
                     tGTAATGGCAGGGGGTGCGGCGTCATACGGGACCTg               >  1:6818063/1‑36 (MQ=255)
                       tAATGGCAGGCGCAGCGCCGTCATACGGGGCATgag             >  1:4425913/1‑36 (MQ=18)
                          tGGCAGGGGCTGCGCCGTCATACGGGACATGAGaaa          <  1:5523604/36‑1 (MQ=255)
                          tGGCAGGGGCTGCGCCGTCATACGGGACATGAGaaa          >  1:3602937/1‑36 (MQ=255)
                          tGGCAGGGGCTGCGCCGTCATACGGGACATGAGaaa          >  1:1739113/1‑36 (MQ=255)
                          tGGCAGGGGATGCGCCGTCATACGGGACATGAGaaa          <  1:6807599/36‑1 (MQ=255)
                           ggCAGGGGCTGCGCCGTCATACGGGACATGAGaaaa         <  1:3510456/36‑1 (MQ=255)
                            gCAGGGGCTGCGCCGTCATACGGGACATGAGAAAAc        >  1:5057301/1‑36 (MQ=255)
                            gCAGGGGCTGCGCCGTCATACGGGACATGAGAAAAc        >  1:3221058/1‑36 (MQ=255)
                            gCAGGGGCTGCGCCGTCATACGGGACATGAGAAAAc        >  1:5275041/1‑36 (MQ=255)
                            gCAGGGGCTGCGCCGTCATACGGGACATGAGAAAAc        >  1:2191258/1‑36 (MQ=255)
                             cAGGGGCTGCGCCGTCATACGGGACATGAGAAAACt       >  1:5871179/1‑36 (MQ=255)
                             cAGGGGCTGCGCCGTCATACGGGACATGAGAAAACt       >  1:7438633/1‑36 (MQ=255)
                              aGGGGCTGCGCCGTCATACGGGACATGAGAAAACtt      >  1:5106779/1‑36 (MQ=255)
                              aGGGGCTGCGCCGTCATACGGGACATGAGAAAACtt      >  1:6220685/1‑36 (MQ=255)
                              aGGGGCTGCGCCGTCATACGGGACATGAGAAAACtt      <  1:5019318/36‑1 (MQ=255)
                               ggggCTGCGCCGTCATACGGGACATGAGAAAACTTc     >  1:5117981/1‑36 (MQ=255)
                               ggggCTGCGCCGTCATACGGGACATGAGAAAACTTc     <  1:4874246/36‑1 (MQ=255)
                                gggCTGCGCCGTCATACGGGACATGAGAAAACTTCg    >  1:2238399/1‑36 (MQ=255)
                                 ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   <  1:707479/36‑1 (MQ=255)
                                 ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   >  1:3046585/1‑36 (MQ=255)
                                 ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   <  1:6621121/36‑1 (MQ=255)
                                 ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   <  1:4219908/36‑1 (MQ=255)
                                 ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   >  1:757598/1‑36 (MQ=255)
                                  gCTGCGCCGTCATACGGGACATGAGAAAACTTCGCa  >  1:5445492/1‑36 (MQ=255)
                                   |                                  
AATATGCCCGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATACGGGACATGAGAAAACTTCGCA  >  REL606/3992543‑3992612

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG