Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,459,924 A→C 19.8% V12V (GTT→GTG zipA ← cell division protein ZipA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,459,9240A→C19.8% 2.0 ‑4.2 0.04036zipAcell division protein ZipA
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (16/14):  total (26/14)
Fisher's exact test strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCATGTACCAGTAAAGCGATTATGGCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCAT  >  REL606/2459889‑2459959
                                   |                                   
ccATGTACCAGTAAAGCGATTATGGCGATCGCGCca                                     <  1:7794725/36‑1 (MQ=255)
ccATGTACCAGTAAAGCGATTATGGCGATCGCGCca                                     <  1:6962728/36‑1 (MQ=255)
 cATGTACCAGTAAAGCGATTATGGCGATCGCGCcaa                                    <  1:5583874/36‑1 (MQ=255)
  aTGTACCAGTAAAGCGATTATGGCGATCGCGCcaac                                   <  1:7001964/36‑1 (MQ=255)
  aTGTACCAGTAAAGCGATTATGGCGATCGCGCcaac                                   <  1:5384491/36‑1 (MQ=255)
  aTGTACCAGTAAAGCGATTATGGCGATCGCGCcaac                                   <  1:1099474/36‑1 (MQ=255)
   tGTACCAGTAAAGCGATTATGGCGATCGCGCcaaca                                  <  1:7378590/36‑1 (MQ=255)
   tGTACCAGTAAAGCGATTATGGCGATCGCGCcaaca                                  >  1:292136/1‑36 (MQ=255)
   tGTACCAGTAAAGCGATTATGGCGATCGCGCcaaca                                  <  1:295430/36‑1 (MQ=255)
    gTACCAGTAAAGCGATTATGGCGATCTCGTccacaa                                 >  1:5206612/1‑36 (MQ=25)
     tACCAGTAAAGCTATTATGGCGATCGCGCCAACAAt                                <  1:6993861/36‑1 (MQ=255)
     tACCAGTAAAGCGATTATGGCGATTGCGCcaacaaa                                >  1:2583219/1‑35 (MQ=255)
     tACCAGTAAAGCGATTATGGCGATCGCGCccacaaa                                >  1:4458468/1‑35 (MQ=255)
     tACCAGTAAAGCGATTATGGCGATCGCGCCAACAAt                                <  1:3124658/36‑1 (MQ=255)
      aCCAGTAAATCGATTATGGCGATCGCGCCAACAATg                               <  1:149063/36‑1 (MQ=255)
      aCCAGTAAAGCGATTATGGCGATCGCGCCACCAATg                               >  1:2344213/1‑36 (MQ=255)
      aCCAGTAAAGCGATTATGGCGACCGCGCCAACAATg                               >  1:752874/1‑36 (MQ=255)
       ccAGTAAAGCGATTATGGCGATCCCGCCCACAATGa                              >  1:7135672/1‑36 (MQ=37)
        cAGTAAAGCGATTATGGCGATCGCGCCAACAATGAt                             <  1:6216842/36‑1 (MQ=255)
         aGTAAAGCGATTATGGCGATTCGGCCCACAATTAtt                            >  1:6633306/1‑36 (MQ=12)
          gTAAAGCGATTATGGCGATCGGGCCCACAAAAATTa                           >  1:4684138/1‑36 (MQ=16)
          gTAAAGCGATTATGGAGAAAACGCAAAGAATGAGTa                           >  1:3537347/1‑36 (MQ=0)
           tAAAGCGATTATGGCGATCGCGCCAACAATGATTaa                          <  1:6511589/36‑1 (MQ=255)
            aaaGCGATTATGGCGATCGCGCCCACAAAGAtttaa                         >  1:7834088/1‑33 (MQ=25)
             aaGCGATTATGGCGATCGCGCCAACAATGATTAata                        <  1:871393/36‑1 (MQ=255)
             aaGCGATTATGGCGATCGCGCCAACAAAGATTTata                        >  1:4490436/1‑36 (MQ=37)
              aGCGATTATGGCGATGCCGGCAACAAAGATTAataa                       >  1:8249885/1‑35 (MQ=16)
              aGCGATTATGGCGATCGCGCCAACAATGATTAatat                       <  1:174422/36‑1 (MQ=255)
               gCGATTATGGCGATCGCGCCCACAATGATTAATATc                      >  1:5458842/1‑36 (MQ=255)
               gCGATTATGGCGATCGCGCCAACAATGATTAATATc                      >  1:8354723/1‑36 (MQ=255)
               gCGATTATGGCGATCGCGCCAACAATGATTAATATc                      <  1:6145854/36‑1 (MQ=255)
               gCGATTATGGCGATCCCGCCCCCAAAAATTTatata                      >  1:5726318/1‑35 (MQ=2)
                cGATTATGGCGATGACGACGAGAATCAGTAatatta                     >  1:6958171/1‑34 (MQ=0)
                cGATTATGGCGATCTCGCCCACAAAGATTTatatta                     >  1:4463445/1‑34 (MQ=16)
                cGATTATGGCGATCGCGCCCACAATAATTAATAACa                     >  1:4991708/1‑36 (MQ=25)
                cGATTATGGCGATCGCGCCAACAATGATTAATATCa                     <  1:2413710/36‑1 (MQ=255)
                cGATTATGGCGATCGCGCCAACAATGAATTatatta                     >  1:7896565/1‑34 (MQ=21)
                 gATTATGGCGATCGCGCCAACAATGATTAATATCAg                    >  1:7668333/1‑36 (MQ=255)
                     aTGGCGATCGCGCCCACAATGATTAATATCAGACGc                >  1:2354717/1‑36 (MQ=255)
                     aTGGCGATCGCGCCAACAATGATTAATATCAGACGc                >  1:5705208/1‑36 (MQ=255)
                      tGGCGATCGCGCCAACAATGATTAATATCAGACGCa               >  1:3939338/1‑36 (MQ=255)
                      tGGCGATCGCGCCAACAATGATTAATATCAGACGCa               <  1:2654737/36‑1 (MQ=255)
                       ggCGATCGCGCCAACAATGATTAATATCAGACGCaa              >  1:811964/1‑36 (MQ=255)
                       ggCGATCGCGCCAACAATGATTAATATCAGACGCaa              <  1:4616958/36‑1 (MQ=255)
                         cGATCGCGCCAACAATGATTAATATCAGACGCAAAt            >  1:6284890/1‑36 (MQ=255)
                          gATCGCGCCAACAATGATTAATATCAGACGCAAATc           >  1:8342997/1‑36 (MQ=255)
                          gATCGCGCCAACAATGATTAATATCAGACGCAAATc           >  1:8268964/1‑36 (MQ=255)
                          gATCGCGCCAACAATGATTAATATCAGACGCAAATc           <  1:6508046/36‑1 (MQ=255)
                            tCGCGCCANCAATGATTAATATCAGACGCAAATCCt         <  1:6563074/36‑1 (MQ=255)
                            tCGCGCCAACAATGATTAATATCAGACGCAAATCCt         <  1:915496/36‑1 (MQ=255)
                              gcgcCAACAATGATTAATATCAGACGCAAATCCTGc       >  1:2968785/1‑36 (MQ=255)
                                 ccAACAATGATTAATATCAGACGCAAATCCTGcatc    >  1:1124797/1‑36 (MQ=255)
                                  caacaaTGATTAATATCAGACGCAAATCCTGcatca   <  1:1632280/36‑1 (MQ=255)
                                  caacaaTGATTAATATCAGACGCAAATCCTGcatca   >  1:818454/1‑36 (MQ=255)
                                  caacaaTGATTAATATCAGACGCAAATCCTGcatca   >  1:3643412/1‑36 (MQ=255)
                                  caacaaTGATTAATATCAGACGCAAATCCTGcatca   >  1:3014256/1‑36 (MQ=255)
                                   aacaaTGATTAATATCAGACGCAAATCCTGcatcat  <  1:4498704/36‑1 (MQ=255)
                                   aacaaTGATTAATATCAGACGCAAATCCTGcatcat  <  1:8327089/36‑1 (MQ=255)
                                   aacaaTGATTAATATCAGACGCAAATCCTGcatcat  <  1:1872613/36‑1 (MQ=255)
                                   aacaaTGATTAATATCAGACGCAAATCCTGcatcat  >  1:1839846/1‑36 (MQ=255)
                                   |                                   
CCATGTACCAGTAAAGCGATTATGGCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCAT  >  REL606/2459889‑2459959

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