Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,062,841 C→A 9.8% L18F (TTG→TTT yccZ ← predicted exopolysaccharide export protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,062,8410C→A 9.8% 2.1 ‑4.3 0.05354yccZpredicted exopolysaccharide export protein
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (34/13):  total (35/18)
Fisher's exact test strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTT  >  REL606/1062807‑1062876
                                  |                                   
aTAAATTCTGACCAGGAACAAGGGTACATGCTGTCa                                    <  1:3436466/36‑1 (MQ=255)
aTAAATTCTGACCAGGAACAAGGGTACATGCTGTCa                                    >  1:2705252/1‑36 (MQ=255)
aTAAATTCTGACCAGGAACAAGGGGACATGCTGTCa                                    >  1:502133/1‑36 (MQ=255)
 tAAATTCTGACCAGGAACAAGGGTACATGCTGTCaa                                   >  1:1394945/1‑36 (MQ=255)
 tAAATTCTGACCAGGAACAAGGGTACATGCTGTCaa                                   <  1:7010760/36‑1 (MQ=255)
  aaaTTCTGACCAGGAACAAGGGTACATGCTGTCaaa                                  <  1:5235082/36‑1 (MQ=255)
  aaaTTCTGACCAGGAACAAGGGTACATGCTGTCaaa                                  <  1:7638640/36‑1 (MQ=255)
    aTTCTGACCAGGAACAAGGGTACATGATGTCAAAGa                                >  1:4177334/1‑36 (MQ=255)
     ttCTGACCAGGAACAAGGGTACATGCTGTCAAAGAc                               >  1:5037258/1‑36 (MQ=255)
     ttCTGACCAGGAACAAGGGTACATGCTGTCAAAGAc                               <  1:2449720/36‑1 (MQ=255)
      tCTGACCAGGAACAAGGGTACATGCTGTCAAAGACa                              >  1:8290063/1‑36 (MQ=255)
      tCTGACCAGGAACAAGGGTACATGCTGTCAAAGACa                              >  1:5574981/1‑36 (MQ=255)
      tCTGACCAGGAACAAGGGTACATGCTGTCAAAGACa                              >  1:8031878/1‑36 (MQ=255)
       cTGACCAGGAACAAGGGTACATGCTGTCAAAGACAg                             <  1:1806062/36‑1 (MQ=255)
        tGACCAGGAACAAGGGTACATGCTGTCAAAGACAGa                            >  1:234146/1‑36 (MQ=255)
         gACCAGGAACAAGGGTACATGCTGTCAAAGACAGaa                           <  1:2998872/36‑1 (MQ=255)
         gACCAGGAACAAGGGTACATGCTGTCAAAGACAGaa                           >  1:8168760/1‑36 (MQ=255)
         gACCAGGAACAAGGGTACATGCTGTCAAAGACAGaa                           >  1:4564125/1‑36 (MQ=255)
         gACCAGGAACAAGGGTACATGCTGTCAAAGACAGaa                           <  1:538735/36‑1 (MQ=255)
         gACCAGGAACAAGGGTACATGCTGTCAAAGACAGaa                           >  1:5795022/1‑36 (MQ=255)
         gACCAGGAACAAGGGTACATGCTGTCAAAGACAGaa                           >  1:7417883/1‑36 (MQ=255)
          aCCAGGAACAAGGGTACATGCTGTCAAAGACAGAAc                          <  1:4149377/36‑1 (MQ=255)
          aCCAGGAACAAGGGTACATGCTGTCAAAGACAGAAc                          >  1:4928730/1‑36 (MQ=255)
          aCCAGGAACAAGGGTACATGCTGTCAAAGACAGAAc                          >  1:1681540/1‑36 (MQ=255)
           ccAGGAACAAGGGTACATGCTGTCAAAGACAGAACa                         <  1:8200615/36‑1 (MQ=255)
            cAGGAACAAGGGTACATGCTGTCAAAGACAGAACAg                        <  1:1807903/36‑1 (MQ=255)
            cAGGAACAAGGGTACATGCTGTCAAAGACAGAACAg                        >  1:7877265/1‑36 (MQ=255)
             aGGAACAAGGGTACATGCTGTCAAAGACAGAACAGc                       <  1:3881185/36‑1 (MQ=255)
              ggAACAAGGGTACATGCTGTCAAAGACAGAACAGcc                      <  1:7016733/36‑1 (MQ=255)
              ggAACAAGGGTACATGCTGTCAAAGACAGAACAGcc                      >  1:679606/1‑36 (MQ=255)
              ggAACAAGGGTACATGCTGTCAAAGACAGAACAGcc                      >  1:1508669/1‑36 (MQ=255)
               gAACAAGGGTACATGCTGTCAAAGACAGAACAGCCa                     >  1:5714973/1‑36 (MQ=255)
                 aCAAGGGTACATGCTGTCAAAGACAGAACAGCCAAt                   >  1:5362168/1‑36 (MQ=255)
                 aCAAGGGTACATGCTGTCAAAGACAGAACAGCCAAt                   >  1:6220713/1‑36 (MQ=255)
                   aaGGGTACATGCTGTCAAAGACAGAACAGCCAATGt                 >  1:2867428/1‑36 (MQ=255)
                   aaGGGTACATGCTGTCAAAGACAGAACAGCCAATGt                 >  1:6985384/1‑36 (MQ=255)
                     gggTACATGCTGAAAAAGAAAGAACAGCCAATGTCa               <  1:212605/36‑1 (MQ=25)
                     gggTACATACTGAAAAAGAAAAACCAGCCAATGTCa               <  1:2497959/36‑1 (MQ=9)
                     gggTAAATGCTGAAAAAGAAAAACCAGCCAATGTCa               <  1:6519962/36‑1 (MQ=9)
                       gTACATGCTGTCAAAGACAGAACAGCCAATGTCAGt             >  1:23755/1‑36 (MQ=255)
                        tACATGCTGTCAAAGACAGAACAGCCAATGTCAGTa            >  1:2113349/1‑36 (MQ=255)
                        tACATGCTGTCAAAGACAGAACAGCCAATGTCAGTa            >  1:5483664/1‑36 (MQ=255)
                         aCATGCTGAAAAAGAAAGACCAGCCAATGTCAGTAc           <  1:1510622/36‑1 (MQ=19)
                          cATGCTGTCAAAGACAGAACAGCCAATGTCAGTAcc          >  1:3261405/1‑36 (MQ=255)
                          cATGCTGTCAAAGACAGAACAGCCAATGTCAGTAcc          >  1:8210541/1‑36 (MQ=255)
                          cATGCTGAAAAAGAAAAAACAGCCAATGTCAGTAcc          <  1:4943600/36‑1 (MQ=16)
                          cATGCTCTCAAAGACAGAACAGCCAATGTCAGTAcc          >  1:7477280/1‑36 (MQ=255)
                           aTGCTGTCAAAGACAGAACAGCCAATGTCAGTACCg         >  1:916328/1‑36 (MQ=255)
                           aTGCTGTCAAAGACAGAACAGCCAATGTCAGTACCg         >  1:1405582/1‑36 (MQ=255)
                           aTGCTGTCAAAGACAGAACAGCCAATGTCAGTACCg         >  1:1321471/1‑36 (MQ=255)
                           aTGCTGTAAAAGACAGAACAGCCAATGTCAGTACCg         >  1:4859747/1‑36 (MQ=255)
                                gTCAAAGACAGAACAGCCAATGTCAGTACCGAAAAc    >  1:2705114/1‑36 (MQ=255)
                                gTCAAAGACAGAACAGCCAATGTCAGTACCGAAAAc    >  1:2021116/1‑36 (MQ=255)
                                  cAAAGACAGAACAGCCAATGTCAGTACCGAAAACtt  >  1:4684915/1‑36 (MQ=255)
                                  cAAAGACAGAACAGCCAATGTCAGTACCGAAAACtt  >  1:4274559/1‑36 (MQ=255)
                                  cAAAGACAGAACAGCCAATGTCAGTACCGAAAACtt  >  1:132802/1‑36 (MQ=255)
                                  |                                   
ATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTT  >  REL606/1062807‑1062876

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG