Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,734,196 T→G 8.0% L411R (CTG→CGG)  pykF → pyruvate kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,734,1960T→G 8.0% 6.2 ‑0.4 ‑0.8501232pykFpyruvate kinase
Reads supporting (aligned to +/- strand):  new base (1/3):  ref base (19/27):  total (20/30)
Fisher's exact test strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.59e-01

GACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGA  >  REL606/1734161‑1734231
                                   |                                   
gACCACCAACGAAAAAACGGCTCATCAGTTGGTACt                                     >  1:923753/1‑36 (MQ=255)
gACCACCAACGAAAAAACGGCTCATCAGTTGGTACt                                     >  1:1002872/1‑36 (MQ=255)
 accaccAACGAAAAAACGGCTCATCAGTTGGTACTg                                    >  1:7313068/1‑36 (MQ=255)
  ccaccaACGAAAAAACGGCTCATCAGTTGGTACTGa                                   >  1:5150099/1‑36 (MQ=255)
   caccaACGAAAAAACGGCTCATCAGTTGGTACTGAg                                  >  1:6747648/1‑36 (MQ=255)
    accaACGAAAAAACGGCTCATCAGTTGGTACTGAGc                                 >  1:1323348/1‑36 (MQ=255)
     ccaACGAAAAAACGGCTCATCAGTTGGTACTTAGCa                                <  1:5385121/36‑1 (MQ=255)
     ccaACGAAAAAACGGCTCATCAGTTGGTACTGAGCa                                >  1:3857304/1‑36 (MQ=255)
      caACGAAAAAACGGCTCATCAGTTGGTACTGAGCaa                               <  1:4002408/36‑1 (MQ=255)
        aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAg                             <  1:5633273/36‑1 (MQ=255)
        aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAg                             >  1:4466882/1‑36 (MQ=255)
        aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAg                             >  1:3660620/1‑36 (MQ=255)
        aCGAAAAAACGGCTCATCAGTTGGTACGGAGCAAAg                             <  1:4511143/36‑1 (MQ=255)
          gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGc                           <  1:6081951/36‑1 (MQ=255)
          gAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGc                           <  1:3143016/36‑1 (MQ=255)
          gAAAAAACGGCTCATCAGCTGGTACTGAGCAAAGGc                           <  1:237058/36‑1 (MQ=255)
           aaaaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCg                          >  1:1621766/1‑36 (MQ=255)
              aaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTtg                       >  1:273243/1‑36 (MQ=255)
               aaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTtgt                      <  1:5886544/36‑1 (MQ=255)
               aaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTtgt                      <  1:6284735/36‑1 (MQ=255)
               aaCGGCTCATCAGTTGGTACGGAGCAAAGGCGTtgt                      <  1:418471/36‑1 (MQ=255)
                 cGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGc                    >  1:6741384/1‑36 (MQ=255)
                   gCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCg                  <  1:7039671/36‑1 (MQ=255)
                   gCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCg                  >  1:6156935/1‑36 (MQ=255)
                    cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGc                 >  1:2677298/1‑36 (MQ=255)
                     tcatcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCa                >  1:3313235/1‑36 (MQ=255)
                      catcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAg               <  1:5158993/36‑1 (MQ=255)
                        tcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCt             <  1:5742025/36‑1 (MQ=255)
                        tcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCt             >  1:5864972/1‑36 (MQ=255)
                         caGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCtt            >  1:2248064/1‑36 (MQ=255)
                         caGTTGGTACGGAGCAAAGGCGTTGTGCCGCAGCtt            >  1:6027547/1‑36 (MQ=255)
                         caGGTGGTACGGAGCAAAGGCGTTGTGCCGCAGCtt            <  1:2504935/36‑1 (MQ=37)
                          aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg           <  1:7425373/36‑1 (MQ=255)
                          aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg           <  1:195354/36‑1 (MQ=255)
                          aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg           <  1:1826069/36‑1 (MQ=255)
                          aGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg           <  1:1313646/36‑1 (MQ=255)
                           gTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGt          <  1:84156/36‑1 (MQ=255)
                            ttGGTACTGAGCAAAGNCGTTGTGCCGCAGCTTGtt         <  1:2570395/36‑1 (MQ=255)
                            ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGtt         <  1:6220725/36‑1 (MQ=255)
                            ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGtt         <  1:5661131/36‑1 (MQ=255)
                            ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGtt         >  1:5136807/1‑36 (MQ=255)
                            ttGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGtt         <  1:2058557/36‑1 (MQ=255)
                             tGGTACTGAGCATAGGCTTTGTGCCGCAGCTTGTTa        <  1:7122742/36‑1 (MQ=37)
                             tGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTa        >  1:7029656/1‑36 (MQ=255)
                              ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaa       <  1:5705124/36‑1 (MQ=255)
                              ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaa       <  1:6513499/36‑1 (MQ=255)
                              ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaa       >  1:2661616/1‑36 (MQ=255)
                               gTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaaa      <  1:7187930/36‑1 (MQ=255)
                                 aCTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAaga    <  1:6645385/36‑1 (MQ=255)
                                 aCTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAaga    <  1:311395/36‑1 (MQ=255)
                                 aCTGAGCAAAGGCGCTGTGCCGCAGCTTGTTAAaga    >  1:68750/1‑36 (MQ=255)
                                  cTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag   <  1:4578258/36‑1 (MQ=255)
                                   tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAgaga  >  1:5609545/1‑36 (MQ=255)
                                   tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAgaga  <  1:1952078/36‑1 (MQ=255)
                                   tGAGCAAAGGCGCTGTGCCGCAGCTTGTTAAAgaga  >  1:6332651/1‑36 (MQ=255)
                                   |                                   
GACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGA  >  REL606/1734161‑1734231

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