Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,325,522 G→C 22.5% A84G (GCA→GGA)  menD ← 2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,325,5220G→C22.5% 4.0 ‑3.9 0.025251menD2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (11/8):  total (11/14)
Fisher's exact test strand distribution p-value = 1.96e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACC  >  REL606/2325490‑2325551
                                |                             
aCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCa                            >  1:1163429/1‑36 (MQ=255)
aCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCa                            >  1:6721021/1‑36 (MQ=255)
aCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCa                            <  1:6524912/36‑1 (MQ=255)
aCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCa                            <  1:6111574/36‑1 (MQ=255)
  ccGGCTTCAATCAGTGCCGGATAGAGATTTGCCAcc                          <  1:7586341/36‑1 (MQ=255)
  ccGGCTTCAATCAGTGCCGGATAGAGATTTGCCAcc                          <  1:3569150/36‑1 (MQ=255)
   cGGCTTCAATCAGTGCCGGATAGAGATTTGCCAccg                         >  1:3775923/1‑36 (MQ=255)
      ctccAATCAGTCCGGAATAGAAATTTGCCAccgccg                      <  1:6541931/33‑1 (MQ=16)
      cTTCAATCAGTGCCGGATAGAGATTTGCCAccgccg                      >  1:2048741/1‑36 (MQ=255)
       ttCAATCAGTGCCGGATAGAGATTTGCCACCGCCGt                     >  1:229193/1‑36 (MQ=255)
               gTGCCGGAAAGAGTTTTCCCACCGCCGTGCCGGAgg             <  1:8457034/36‑1 (MQ=25)
                tGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGt            >  1:3784033/1‑36 (MQ=255)
                tGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGt            >  1:3604866/1‑36 (MQ=255)
                ccccGGAAAAAAATTTGCCCCCGCCGTGCCGGAGGt            <  1:5789976/34‑1 (MQ=16)
                   cGGATAGAGTTTGCCCCCCGCCGTGCCGGAGGTcac         <  1:7931785/36‑1 (MQ=19)
                    ggATAGAGATTTGCCACCGCCGTGCCGGAGGTcaca        >  1:4146513/1‑36 (MQ=255)
                    ggATAGAGATTTGCCACCGCCGTGCCGGAGGTcaca        >  1:1861727/1‑36 (MQ=255)
                      atggagaTTCGCCACCGCCGTGCCGGAGGTCACAAt      <  1:6165771/33‑1 (MQ=25)
                      aTAAAGTTTTCCCCCCGCCGTGCCGGAGGTCACAAt      <  1:4190004/36‑1 (MQ=19)
                       tAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATc     >  1:6880283/1‑36 (MQ=255)
                       aagagTTTTGCCACCGCCGTGCCGGAGGTCACAATc     <  1:7565731/35‑1 (MQ=255)
                         gagaTTTGCCACCGCCGTGCCGGAGGTCACAATCAc   >  1:2819229/1‑36 (MQ=255)
                         gagaTTTCCCCCCGCCGTGCCGGAGGTCACAATCAc   <  1:1174855/36‑1 (MQ=37)
                         aagaTTTCCCCCCGCCGTGCCGGAGGTCACAATCAc   <  1:7192040/35‑1 (MQ=37)
                          agaTTGCCCACCGCCGTGCCGGAGGTCACAATCAcc  <  1:1174597/36‑1 (MQ=37)
                                |                             
ACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACC  >  REL606/2325490‑2325551

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