Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,020,046 A→G 15.9% I328V (ATC→GTC)  ECB_02819 → KpsS protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,020,0460A→G15.9% 4.7 ‑3.9 0.066982ECB_02819KpsS protein
Reads supporting (aligned to +/- strand):  new base (11/1):  ref base (28/26):  total (39/27)
Fisher's exact test strand distribution p-value = 2.01e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CAAAAGCGGTGGTGACAATTAACAGTACGGCGGGGATCTCTGCGCTGATTCATAACAAACCACTCAAAGT  >  REL606/3020011‑3020080
                                   |                                  
aaaaaGCGGTGGTGACAATTAACAGTACGGCGGGGa                                    <  1:1478836/35‑1 (MQ=255)
  aaaGCGGTGGTGACAATTAACAGTACGGCGGGGGtc                                  >  1:2466976/1‑36 (MQ=255)
  aaaGCGGTGGTGACAATTAACAGTACGGCGGGGGtc                                  >  1:6043032/1‑36 (MQ=255)
  aaaGCGGTGGTGACAATTAACAGTACGGCGGGGAtt                                  >  1:2546159/1‑35 (MQ=255)
   aaGCGGTGGTGACAATTAACAGTACGGCGGGGGtct                                 >  1:3476519/1‑36 (MQ=255)
   aaGCGGTGGTGACAATTAACAGTACGGCGGGGGtct                                 >  1:7184717/1‑36 (MQ=255)
    aGCGGTGGTTACAATTAACAGTACGGCGGGGAtctc                                <  1:3422581/36‑1 (MQ=255)
    aGCGGTGGTGACAATTAACAGTACGGCGGGGAtctc                                <  1:6968579/36‑1 (MQ=255)
     gCGGTGGTGACAATTAACAGTACGGCGGGGATctct                               <  1:7350856/36‑1 (MQ=255)
      cGGTGGTGACAATTAACAGTACGGCGGGGATCTCTg                              <  1:1655836/36‑1 (MQ=255)
       ggtggtGACAATTAACAGTACTGCGGGGATCTCTgc                             <  1:1937994/36‑1 (MQ=255)
       ggtggtGACAATTAACAGTACGGCGGGGGTctcttc                             >  1:3172790/1‑34 (MQ=38)
       ggtggtGACAATTAACAGTACGGCGGGGATCTCTgc                             >  1:1605067/1‑36 (MQ=255)
       ggtggtGACAATTAACAGTACGGCGGGGATCTCTgc                             <  1:6799563/36‑1 (MQ=255)
        gtggtgACAATTAACAGTACGGCGGGGGTCTCTtcg                            >  1:6886056/1‑36 (MQ=37)
         tggtgACAATTAACAGTACGGCGGGGATCTCTgcgc                           >  1:6896395/1‑36 (MQ=255)
         tggtgACAATTAACAGTACGGCGGGGATCTCTgcgc                           >  1:7219054/1‑36 (MQ=255)
         tggtgACAATTAACAGTACGGCGGGGATCTCTgcgc                           <  1:1266656/36‑1 (MQ=255)
          ggtgACAATTAACAGTACGGCGGGGATCTCTGCGCt                          <  1:6863436/36‑1 (MQ=255)
           gtgACAATTAACAGTACGGCGGGGGTATTTTCGCTg                         >  1:4488341/1‑36 (MQ=19)
           gtgACAATTAACAGTACGGCGGGGATCTTTGCGCTg                         >  1:2492075/1‑36 (MQ=255)
           gtgACAATTAACAGTACGGCGGGGATCTCTGCGCTg                         <  1:6354697/36‑1 (MQ=255)
           gtgACAATTAACAGTACGGCGGGGATCTCTGCGCTg                         >  1:7818875/1‑36 (MQ=255)
            tgACAATTAACAGTACGGCGGGGATCTCTGCGCTGa                        <  1:6666179/36‑1 (MQ=255)
            tgACAATTAACAGTACGGCGGGGATCTCTGCGCTGa                        <  1:6656880/36‑1 (MQ=255)
              aCAATTAACAGTACGGCGGGGGTCTCTTCGCTGGtt                      >  1:5790852/1‑36 (MQ=25)
              aCAATTAACAGTACGGCGGGGATCTCTGCGCTGAtt                      <  1:2147798/36‑1 (MQ=255)
              aCAATTAACAGTACGGCGGGGATCTCTGCGCTGAtt                      <  1:4032408/36‑1 (MQ=255)
              aCAATTAACAGTACGGCGGGGATCTCTGCGCTGAtt                      <  1:1035654/36‑1 (MQ=255)
               cAATTAACAGTACGGCGGGGATCTCTGCGCTGATTc                     <  1:7037579/36‑1 (MQ=255)
                caTTAACAGTACGGCGGGGATTTCTGCGCTGATTCa                    >  1:5689900/2‑36 (MQ=255)
                aaTTAACAGTACGGCGGGGATCTCTGCGCTGATTCa                    >  1:5624103/1‑36 (MQ=255)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTTATTTAt                   >  1:5035119/1‑36 (MQ=37)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTGATTTAt                   >  1:4880444/1‑36 (MQ=255)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTGATTCAt                   <  1:5857520/36‑1 (MQ=255)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTGATTCAt                   <  1:1981599/36‑1 (MQ=255)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTGATTCAt                   <  1:2531767/36‑1 (MQ=255)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTGATTCAt                   <  1:7641888/36‑1 (MQ=255)
                 aTTAACAGTACGGCGGGGATCTCTGCGCTGATTCAt                   <  1:3243171/36‑1 (MQ=255)
                  ttAACAGTACGGCGGGGATCTCTGCGCTGATTCATa                  <  1:47364/36‑1 (MQ=255)
                   tAACAGTACGGCGGGGATCTCTGCGCTGCTTCATaa                 >  1:712688/1‑36 (MQ=255)
                    aaCAGTACGGCGGGGATCTCTGCGCTGATTTATaaa                >  1:7492729/1‑35 (MQ=255)
                     aCAGTACGGCGGGGGTCTCTGCGCTGCTTCATAACa               <  1:121225/36‑1 (MQ=255)
                     aCAGTACGGCGGGGGTCTCTGCGCTGATTCTTaaaa               >  1:2931612/1‑34 (MQ=37)
                     aCAGTACGGCGGGGGTATCTGCGTTGGTTTATAACa               >  1:1034221/1‑36 (MQ=255)
                     aCAGTACGGCGGGGATTTTTGGGCTTGTTCATAACa               >  1:4626060/1‑36 (MQ=16)
                     aCAGTACGGCGGGGATCTCTGCGCTGGTTCATAACa               >  1:2848466/1‑36 (MQ=255)
                     aCAGTACGGCGGGGATCTATGCGCTGGTTTATaaaa               >  1:898063/1‑34 (MQ=25)
                      cAGTACGGCGGGGGTCTCTGCGCTGATTCATAACaa              >  1:6830605/1‑36 (MQ=255)
                      cAGTACGGCGGGGATCTCTGCGCTGATTCATAACaa              >  1:5842054/1‑36 (MQ=255)
                      cAGTACGGCGGGGATCTCTGCGCTGATTCATAACaa              <  1:4721952/36‑1 (MQ=255)
                      cAGTACGGCGGGGATCTCTGCGCTGATTCATAACaa              >  1:3695187/1‑36 (MQ=255)
                       aGTACGGCGGGGATCTCTGCGCTGGTTCATAACaaa             >  1:7496070/1‑36 (MQ=255)
                        gTACGGCGGGGATCTCTGCGCTGATTCATAACAAAc            >  1:7446115/1‑36 (MQ=255)
                          aCGGCGGGGATCTCTTCGCTGATTCATAACAAACCa          >  1:6868621/1‑36 (MQ=255)
                           cggcggGTATCTCTGCGCTTATTCATAACAAACCAc         <  1:4486284/36‑1 (MQ=37)
                           cggcggGGATCTCTGCGCTGATTCATAACAAACCAc         <  1:4981925/36‑1 (MQ=255)
                             gcggGGATCTCTGCGCTGATTCATAACAAACCACTc       >  1:2100520/1‑36 (MQ=255)
                             gcggGGATCTCTGCGCTGATTCATAACAAACCACTc       >  1:4310482/1‑36 (MQ=255)
                              cggGGATCTCTGCGCTGATTCATAACAAACCACTCa      >  1:649060/1‑36 (MQ=255)
                               ggggATCTCTGCGCTGATTCATAACAAACCACTCaa     <  1:2343960/36‑1 (MQ=255)
                               ggggATCTCTGCGCTGATTCATAACAAACCACTCaa     <  1:7937102/36‑1 (MQ=255)
                               ggggATCTCTGCGCTGATTCATAACAAACCACTCaa     >  1:6792750/1‑36 (MQ=255)
                               ggggATCTCTGAGCTGATTCATAACAAACCACTCaa     >  1:3815329/1‑36 (MQ=255)
                                gggATCTCTGCGCTGATTCATAACAAACCACTCaaa    >  1:7349390/1‑36 (MQ=255)
                                 ggATCTCTGCGCTGATTCATAACAAACCACTCAAAg   >  1:3348031/1‑36 (MQ=255)
                                  gATCTCTGCGCTGATTCATAACAAACCACTCAAAGt  >  1:7423012/1‑36 (MQ=255)
                                   |                                  
CAAAAGCGGTGGTGACAATTAACAGTACGGCGGGGATCTCTGCGCTGATTCATAACAAACCACTCAAAGT  >  REL606/3020011‑3020080

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