Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 135,246 +G 5.6% coding (828/2598 nt) acnB → aconitate hydratase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606135,2461.→G 5.6% 2.2 ‑2.4 ‑0.036coding (828/2598 nt)acnBaconitate hydratase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (9/25):  total (11/25)
Fisher's exact test strand distribution p-value = 8.73e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.71e-01

GAACAAACGCGGCGGTGGTTTGTGCCTCGGCGGT.AAAATTGC.ACCAATCTTCTTTAACACAATGGAAGA  >  REL606/135213‑135281
                                  |                                    
gAACAAACGCGGCGGTGGTTTTTGCCTCGGCGGT.aa                                    >  1:7540354/1‑36 (MQ=255)
gAACAAACGCGGCGGTGGTTTGTGCCTCGGCGGT.aa                                    <  1:1160461/36‑1 (MQ=255)
caaCAAACGCGGCGGTGGTTTTTGCCTCGGCGGT.aa                                    <  1:2858173/35‑1 (MQ=255)
    aaaCGCGGCGGTGGTTTGTGCCTCGGCGGT.AAAAtt                                <  1:6028558/36‑1 (MQ=255)
    aaaCGCGGCGGTGGTTTGTGCCTCGGCGGT.AAAAtt                                <  1:3393065/36‑1 (MQ=255)
     aaCGCGGCGGTGGTTTTTGCCTCGGCGGTGAAAAtt                                >  1:3630554/1‑36 (MQ=25)
      aCGCGGCGGTGGTTTGTGCCTCGGCGGT.TAAAtttc                              >  1:3895755/1‑34 (MQ=37)
            cggTGGTTTGTGCCTCGGCGGT.AAAATTGC.ACCAAt                       <  1:1929169/36‑1 (MQ=255)
            cggTGGTTTGTGCCTCGGCGGT.AAAATTGC.ACCAAt                       <  1:3231282/36‑1 (MQ=255)
             ggtggtTTGTGCCTCGGCGGT.AAAATTGT.ACCAAtt                      >  1:6474944/1‑35 (MQ=255)
             ggtggtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtc                      <  1:1369241/36‑1 (MQ=255)
             ggtggtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtc                      <  1:5739952/36‑1 (MQ=255)
              gtggtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtct                     <  1:3517853/36‑1 (MQ=255)
                ggtTTGTGCCTCGGCGGT.CAAATTGCAACCAAttgt                    >  1:5632091/1‑33 (MQ=16)
                ggtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtcttc                   <  1:3190467/36‑1 (MQ=255)
                ggtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtcttc                   >  1:231437/1‑36 (MQ=255)
                 gtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtcttct                  <  1:1918880/36‑1 (MQ=255)
                 gtTTGTGCCTCGGCGGT.AAAATTGC.ACCAAtcttct                  <  1:4964643/36‑1 (MQ=255)
                  tttGTGCCTCGGCGGTGAAAATTTC..CCAAActtctt                 >  1:2021891/1‑36 (MQ=9)
                   ttGTGCCTCGGCGGT.AAAATTGC.ACCAATCTTCttt                <  1:3306922/36‑1 (MQ=255)
                     gtgCCTCGGCGGT.AAAATTGC.ACCAATCTTCTTTaa              <  1:716826/36‑1 (MQ=255)
                     gtgCCTCGGCGGT.AAAATTGC.ACCAATCTTCTTTaa              <  1:1729523/36‑1 (MQ=255)
                      tgCCTCGGCGGT.AAAATTTC.ACCAATCTTCTTTTac             >  1:127766/1‑36 (MQ=37)
                      tgCCTCGGCGGT.AAAATTGC.CCCAATCTTCTTTAac             <  1:4387660/36‑1 (MQ=255)
                      tgCCTCGGCGGT.AAAATTGC.ACCAATCTTCTTTAac             <  1:1447050/36‑1 (MQ=255)
                        ccTCGGCGGT.AAAATTTC.ACCAATCTTCTTTAacac           >  1:1257325/1‑36 (MQ=255)
                         cgcggcggT.AAAATGGC.ACCGAGCTTCTTTAAcaca          <  1:3597096/34‑1 (MQ=21)
                           cggcggT.AAAATTGC.ACCAATCTTCTTTAACACAAt        <  1:6401215/36‑1 (MQ=255)
                           cggcggT.AAAAGTGC.ACCACTCTTCTTTAACACAAt        <  1:3089453/36‑1 (MQ=37)
                           cggcggT.AAAACTGC.ACCAATCTTCTTTAACACAAt        <  1:7508062/36‑1 (MQ=255)
                            ggcggT.AAAATTGC.CCCAATCTTCTTTAACACAATg       <  1:44930/36‑1 (MQ=255)
                            ggcggT.AAAATTGC.ACCAATCTTCTTTAACACAATg       <  1:5021725/36‑1 (MQ=255)
                             gcggT.AAAATTGC.ACCAATCTTCTTTAACACAATgg      <  1:5669399/36‑1 (MQ=255)
                             gcggT.AAAATTGC.ACCAATCTTCTTTAACACAATgg      >  1:3103117/1‑36 (MQ=255)
                               ggT.AAAATTGC.ACCAATCTTCTTTAACACAATGGaa    <  1:1556119/36‑1 (MQ=255)
                               ggT.AAAATTGC.ACCAATCTTCTTTAACACAATGGaa    <  1:1564923/36‑1 (MQ=255)
                               ggT.AAAATTGC.ACCAATCTTCTTTAACACAATGGaa    >  1:2329346/1‑36 (MQ=255)
                                gT.AAAATTGC.ACCAATCTTCTTTAACACAATGGAAg   <  1:861369/36‑1 (MQ=255)
                                gT.AAAATTGA.ACCAATATTCTTTAACACAATGGAAg   >  1:7052885/1‑36 (MQ=37)
                                 t.AAAATTGC.ACCAATCTTCTTTAACACAATGGAAGa  <  1:6414438/36‑1 (MQ=255)
                                 t.AAAATTGC.ACCAATCTTCTTTAACACAATGGAAGa  >  1:7140906/1‑36 (MQ=255)
                                  |                                    
GAACAAACGCGGCGGTGGTTTGTGCCTCGGCGGT.AAAATTGC.ACCAATCTTCTTTAACACAATGGAAGA  >  REL606/135213‑135281

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG