Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,014,308 G→T 17.7% 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→T17.7% 7.3 ‑3.5 0.036889ECB_02815Capsule polysaccharide export inner‑membrane protein kpsE
Reads supporting (aligned to +/- strand):  new base (6/1):  ref base (10/19):  total (16/20)
Fisher's exact test strand distribution p-value = 2.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGGGTGACAAGCTAAACCGTATGGCGG.TGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAG  >  REL606/3014273‑3014343
                                    |                                   
aGGGTGACAAGCTAAACCGTATGGCGG.TGGATTTTg                                     <  1:1023590/36‑1 (MQ=255)
aGGGTGACAAGATAAACCGTATGGCGG.TGGATTTTg                                     <  1:6436881/36‑1 (MQ=255)
 gggTGACAAGCTAAACCGTATGGCGG.TGGATTTTga                                    <  1:2268705/36‑1 (MQ=255)
     gACAAGCTAAACCGTATGGCGG.TGGATTTTgaagaa                                <  1:3047032/36‑1 (MQ=255)
       cAAGCTAAACCGTATGGCGGCTGGTTTTTGTAAaaa                               >  1:5523242/1‑36 (MQ=12)
       cAAGCTAAACCGTATGGCGG.GGGGTTTTTTAGaaaa                              >  1:6437428/1‑35 (MQ=19)
       aaaGCTAAACCGTATGGCGG.TGGATTTTTAAGaaat                              <  1:3167832/35‑1 (MQ=255)
        aaGCTAAACCGTATGGCGG.TGGATTTTGAAGaaatc                             <  1:5646143/36‑1 (MQ=255)
         aGCTAAACCGTATGGCGG.AGGGTTTTTAAAaaaaaa                            >  1:893063/1‑33 (MQ=16)
            tAAACCGTATGGCGG.TGGATTTTGAAGaaatcaaat                         <  1:5611623/36‑1 (MQ=255)
            tAAACCGTATGGCGG.GGGGTTTTTTATaaatcaaat                         >  1:1479177/1‑36 (MQ=16)
             aaaCCGTATGGCGG.TGGATTTTGAAGaaatcaaatc                        <  1:5965185/36‑1 (MQ=255)
              aaCCGTATGGCGG.TGGATTTTGAAGAaatcaaatca                       <  1:41939/36‑1 (MQ=255)
               aCCGTATGGCGG.TGGGTTTTTTAGaaatcaaatccc                      >  1:2616391/1‑34 (MQ=25)
                ccGTATGGCGG.TGGGATTTTTAGAAAataaattaaa                     >  1:5019895/1‑36 (MQ=2)
                ccGTATGGCGG.TGGATTTTGAAGAAAtcaaatcaaa                     <  1:7158016/36‑1 (MQ=255)
                  gTATGGCGG.TGGATTTTTAAGAAATCAAATCAAAAg                   >  1:7423678/1‑36 (MQ=255)
                  gTATGGCGG.TGGATTTTGAAGAAATCAAATCAAAAg                   <  1:6096363/36‑1 (MQ=255)
                   tATGTCGG.TGGATTTTGAAGAAATCAAATCAAAAGt                  <  1:2306725/36‑1 (MQ=255)
                    aTGGCGG.TGGATTTTGAAGAAATCAAATCAAAAGTa                 <  1:781235/36‑1 (MQ=255)
                     tGGCGG.TGGATTTTGAAGAAATCAAATCAAAAGTag                <  1:638762/36‑1 (MQ=255)
                      ggCGG.TGGATTTTGAAGAAATCAAATCAAAAGTaga               >  1:522045/1‑36 (MQ=255)
                      ggCGG.TGGATTTTGAAGAAATCAAATCAAAAGTaga               <  1:7618035/36‑1 (MQ=255)
                          g.TGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTc           >  1:7035929/1‑36 (MQ=255)
                          g.TGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTc           <  1:690710/36‑1 (MQ=255)
                          g.TGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTc           >  1:4480968/1‑36 (MQ=255)
                             ggATTTTGAAGAAATCAAATCAAAAGTAGAGTTCaa         <  1:7391952/36‑1 (MQ=255)
                             ggATTTTGAAGAAATCAAATCAAAAGTAGAGTTCaa         >  1:1558301/1‑36 (MQ=255)
                              gATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAac        >  1:7352684/1‑36 (MQ=255)
                              gATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAac        >  1:1893644/1‑36 (MQ=255)
                              gATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAac        <  1:2757673/36‑1 (MQ=255)
                               aTTTTGAAGAAATCAAATCAAAAGTAGAGTTCAaca       <  1:4851471/36‑1 (MQ=255)
                                ttttGAAGAAATCAAATCAAAAGTAGAGTTCAacac      <  1:2564547/36‑1 (MQ=255)
                                ttttGAAGAAATCAAATCAAAAGTAGAGTTCAacac      <  1:3331163/36‑1 (MQ=255)
                                  ttGAAGAAATCAAATCAAAAGTAGAGTTCAacacag    >  1:5907005/1‑34 (MQ=255)
                                  ttGAAGAAATCAAATCAAAAGTAGAGTTCAACACCg    <  1:6755703/36‑1 (MQ=255)
                                  ttGAAGAAATCAAATCAAAAGTAGAGTTCAACACCg    >  1:5710138/1‑36 (MQ=255)
                                  ttGAAGAAATCAAATCAAAAGTAGAGTTCAACACCg    >  1:491722/1‑36 (MQ=255)
                                  ttGAAGAAATCAAATCAAAAGTAGAGTTCAACACCg    <  1:7780469/36‑1 (MQ=255)
                                   tGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGa   <  1:4676304/36‑1 (MQ=255)
                                    gaagaaATCAAATCAAAAGTAGAGTTCAACACCGAg  >  1:1448952/1‑36 (MQ=255)
                                    gaagaaATCAAATCAAAAGTAGAGTTCAACACCGAg  >  1:6624787/1‑36 (MQ=255)
                                    gaagaaATCAAATCAAAAGTAGAGTTCAACACCGAg  <  1:3595232/36‑1 (MQ=255)
                                    gaagaaATCAAATCAAAAGTAGAGGTCAACACCGAg  >  1:4638475/1‑36 (MQ=255)
                                    |                                   
AGGGTGACAAGCTAAACCGTATGGCGG.TGGATTTTGAAGAAATCAAATCAAAAGTAGAGTTCAACACCGAG  >  REL606/3014273‑3014343

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG