Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,014,308 G→T 33.3% E297* (GAA→TAA)  ECB_02815 → Capsule polysaccharide export inner‑membrane protein kpsE

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,014,3080G→T33.3% 6.1 ‑4.6 0.012889ECB_02815Capsule polysaccharide export inner‑membrane protein kpsE
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (1/7):  total (5/7)
Fisher's exact test strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGGGTGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACC  >  REL606/3014273‑3014340
                                   |                                
aGGGTGACAAGCTAAACCGTATGGCGGTGGATTTTg                                  <  1:3101928/36‑1 (MQ=255)
      aCAAGCTAAACCGTATGGCGGTGGGTTTTTTAGaaa                            >  1:6806691/1‑36 (MQ=37)
      aCAAGCTAAACCGTATGGCGGTGGATTTTTAAGaaa                            >  1:413170/1‑36 (MQ=255)
        aaGCTAAACCGTATGGCGGTGGATTTTGAAGaaatc                          <  1:1341003/36‑1 (MQ=255)
         aGCTAAACCGTATGGCGGTGGATTTTGAAGaaatca                         <  1:1881705/36‑1 (MQ=255)
             aaaCCGTATGGCGGTGGATTTTGAAGaaatcaaatc                     <  1:2124405/36‑1 (MQ=255)
              aaCCGTATGGCGGTGGATTTTGAAGAaatcaaatca                    <  1:6700235/36‑1 (MQ=255)
               aCCGTATGGCGGTGGGTTTTTAAGAAataaaatcaa                   >  1:1404639/1‑36 (MQ=25)
                ccGTATGGCGGTGGATTTTTAAGAAAtcaaatcaaa                  >  1:7558724/1‑36 (MQ=255)
                         ggTGGATTTTGAAGAAATCAAATCAAAAGTAGAGtt         <  1:846310/36‑1 (MQ=255)
                         ggTGGATTTTGAAGAAATCAAATCAAAAGTAGAGtg         <  1:961991/36‑2 (MQ=255)
                            ggATTTTGAAGAAATCAAATCAAAAGTAGAGTTCaa      >  1:8350930/1‑36 (MQ=255)
                                tttGAAGAAATCAAATCAAAAGTAGAGTTCAACAcc  <  1:2279738/36‑1 (MQ=255)
                                   |                                
AGGGTGACAAGCTAAACCGTATGGCGGTGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACC  >  REL606/3014273‑3014340

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG