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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,992,5780C→G20.7% 5.7 ‑2.9 0.041571ECB_03721conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/1):  ref base (15/16):  total (24/17)
Fisher's exact test strand distribution p-value = 5.36e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TATGCCCGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATACGGGACATGAGAAAACTTCGCA  >  REL606/3992545‑3992612
                                 |                                  
tatGCCCGCCTAACAGGCCTGTAATGGCAGGGGGTg                                  >  1:4477792/1‑36 (MQ=255)
 atGCCCGCCTAACAGGCCTGTAATGGCAGGGGccgg                                 >  1:1827877/1‑33 (MQ=255)
   gCCCGCCTAACAGGCCTGTAATGGCAGGGGCTgcgc                               <  1:3981035/36‑1 (MQ=255)
     ccGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCg                             <  1:804666/36‑1 (MQ=255)
     ccGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCg                             <  1:3448744/36‑1 (MQ=255)
      cGCCTAACAGGCCTGTAATGGCAGGGGGTGCGCCGt                            >  1:3697075/1‑36 (MQ=255)
       gCCTAACAGGCCTGTAATGGCAGGGGGTGCGCCGTc                           >  1:1936834/1‑36 (MQ=255)
       gCCTAACAGGCCTGTAATGGCAGGGGGTGCGCCGTc                           >  1:7360183/1‑36 (MQ=255)
       gCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTc                           <  1:5897683/36‑1 (MQ=255)
       gCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTc                           <  1:6386714/36‑1 (MQ=255)
        ccTAACAGGCCTGTAATGGCAGGGGCTGGGGCGTCa                          >  1:2034990/1‑36 (MQ=37)
        ccTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCa                          <  1:5255616/36‑1 (MQ=255)
        ccTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCa                          <  1:6912454/36‑1 (MQ=255)
         cTAACAGGCCTGTAATGGCAGGGGGTGGGCCGTCAt                         >  1:941309/1‑36 (MQ=37)
            aCAGGCCTGTAATGGCAGGGGGTGCGGCGTCATACg                      >  1:4463546/1‑36 (MQ=38)
            aCAGGCCTGTAATGGCAGGGGCTGCGCCGTCATACg                      <  1:4939460/36‑1 (MQ=255)
             aaGGCCTGTAATGGCAGGGGCTTCGCCGTCATCCgg                     <  1:6996632/35‑1 (MQ=37)
                  ttGTAATGGCAGGGGCTGCGCCGTCATACGGGACAt                <  1:1523007/35‑1 (MQ=255)
                  cTGTAATGGCAGGGGGTGCGCCGGCCTTCGGGGCAt                >  1:4097488/1‑36 (MQ=16)
                  cTGTAATGGCAGGGGCTGCGCCGTCATACGGGACAt                <  1:6779006/36‑1 (MQ=255)
                   tGTAGTGGCGGGGGGTGCGCCGTCATACGGGACATg               <  1:5928881/36‑1 (MQ=25)
                    gTAATGGCAGGGGTTGCGCCCCCATACGGGACAtca              >  1:7743691/1‑34 (MQ=25)
                    gTAATGGCAGGGGGTGCGCCGTCATATGGGGCATgg              >  1:6386792/1‑35 (MQ=255)
                      aaTGGCAGGGGCTGCGCCGTCATACGGGACATggga            >  1:2454448/1‑36 (MQ=255)
                      aaTGGCAGGGGCTGCGCCCACTCACGGGGCATgcga            >  1:518744/1‑36 (MQ=12)
                       aTGGCAGGGGCTGCGCCGTCATACGGGGCATGAGaa           >  1:1274410/1‑36 (MQ=255)
                        tGGCAGGGGCTGCGCCGTCATACGGGACATGAGaaa          <  1:5496232/36‑1 (MQ=255)
                        tGGCAGGGGCTGCGCCGTCATACGGGACATGAGaaa          >  1:5169291/1‑36 (MQ=255)
                         ggCAGGGGCTGCGCCGTCATACGGGACATGAGaaaa         >  1:28996/1‑36 (MQ=255)
                          gCAGGGGCTGCGCCGTCATACGGGACATGAGAAAAc        >  1:3034582/1‑36 (MQ=255)
                           cAGGGGCTGCGCCGTCATACGGGACATGAGAAAACt       >  1:5021678/1‑36 (MQ=255)
                           cAGGGGCTGCGCCGTCATACGGGACATGAGAAAACt       >  1:6912940/1‑36 (MQ=255)
                            aGGGGCTGCGCCGTCATACGGGACATGAGAAAACtt      >  1:2691030/1‑36 (MQ=255)
                             ggggCTGCGCCGTCATACGGGACATGAGAAAACTTc     <  1:4975955/36‑1 (MQ=255)
                             ggggCTGCGCCGTCATACGGGACATGAGAAAACTTc     <  1:4969433/36‑1 (MQ=255)
                             ggggCTGCGCCGTCATACGGGACATGAGAAAACTTc     <  1:4141436/36‑1 (MQ=255)
                              gggCTGCGCCGTCATACGGGACATGAGAAAACTTCg    >  1:4169130/1‑36 (MQ=255)
                              gggCTGCGCCGTCATACGGGACATGAGAAAACTTCg    >  1:2493569/1‑36 (MQ=255)
                               ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   >  1:6071179/1‑36 (MQ=255)
                               ggCTGCGCCGTCATACGGGACATGAGAAAACTTCGc   >  1:1390734/1‑36 (MQ=255)
                                gCTGCGCCGTCATACGGGACATGAGAAAACTTCGCa  <  1:6180842/36‑1 (MQ=255)
                                gCTGCGCCGTCATACGGGACATGAGAAAACTTCGCa  >  1:398344/1‑36 (MQ=255)
                                 |                                  
TATGCCCGCCTAACAGGCCTGTAATGGCAGGGGCTGCGCCGTCATACGGGACATGAGAAAACTTCGCA  >  REL606/3992545‑3992612

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