Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,329,291 A→T 20.6% A22A (GCA→GCT elaD → predicted enzyme

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,329,2910A→T20.6% 17.1 ‑0.9 ‑0.23466elaDpredicted enzyme
Reads supporting (aligned to +/- strand):  new base (5/2):  ref base (13/14):  total (18/16)
Fisher's exact test strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.68e-01

CAATTATCTCAAACGCAACTCAGTCAGACATTTGCAGAAAAATTTACTGTGACCGAGGAATTACTGCAGT  >  REL606/2329256‑2329325
                                   |                                  
cAATTATCTCAAACGCAACTCAGTCAGACATTTGCa                                    <  1:5938835/36‑1 (MQ=255)
cAATTATCTCAAACGCAACTCAGTCAGACATTTGCa                                    <  1:6345379/36‑1 (MQ=255)
 aaTTATCTCAAACGCAACTCAGTCAGACATTTGCAg                                   >  1:4991893/1‑36 (MQ=255)
 aaTTATCTCAAACGCAACTCAGTCAGACATTTGCAg                                   >  1:7553793/1‑36 (MQ=255)
  aTTATCTCAAACGCAACTCAGTCAGACATTTGCTGa                                  >  1:3584125/1‑36 (MQ=255)
  aTTATCTCAAACGCAACTCAGTCAGACATTTGCAGa                                  <  1:5855590/36‑1 (MQ=255)
   ttATCTCAAACGCAACTCAGTCAGACATTTGCTGaa                                 >  1:4196293/1‑36 (MQ=255)
   atATCTCAAACGCAACTCAGTCAGACATTTGCAGaa                                 <  1:6903303/35‑1 (MQ=255)
    tATCTCAAACGCAACTCAGTCAGACATTTGCAGaaa                                >  1:7373535/1‑36 (MQ=255)
     aTCTCAAACGCAACTCAGTCAGACATTTGCTGaaaa                               <  1:2526777/36‑1 (MQ=255)
      tctcAAACGCAACTCAGTCAGACATTTGCAGaaaaa                              <  1:7611684/36‑1 (MQ=255)
         cAAACGCAACTCAGTCAGACATTTGCTGAAAAAttt                           >  1:7218158/1‑36 (MQ=255)
          aaaCGCAACTCAGTCAGACATTTGCAGAAAAATTTa                          <  1:7365007/36‑1 (MQ=255)
          aaaCGCAACTCAGTCAGACATTTGCAGAAAAATTTa                          >  1:1661883/1‑36 (MQ=255)
           aaCGCAACTCAGTCAGACATTTGCAGAAAAATTTAc                         >  1:1896989/1‑36 (MQ=255)
           aaCGCAACTCAGTCAGACATTTGCAGAAAAATTTAc                         >  1:3091717/1‑36 (MQ=255)
            aCGCAACACAGTCAGACATTTGCAGAAAAATTTACt                        <  1:5404199/36‑1 (MQ=255)
             cGCAACTCAGTCAGACATTTGCAGAAAAATTTACtg                       <  1:5764011/36‑1 (MQ=255)
             cGCAACTCAGTCAGACATTTGCAGAAAAATTTACtg                       <  1:134464/36‑1 (MQ=255)
             cGCAACTCAGTCAGACATTTGCAGAAAAATTTACtg                       <  1:2263961/36‑1 (MQ=255)
             cGCAACTCAGTCAGACATTTGCAGAAAAATTTACtg                       <  1:1716199/36‑1 (MQ=255)
              gCAACTCAGTCAGACATTTGCAGAAAAATTTACtgt                      >  1:2404762/1‑36 (MQ=255)
               cAACTCAGTCAGACATTTGCTGAAAAATTTACtgtg                     <  1:6669768/36‑1 (MQ=255)
                  cTCAGTCAGACATTTGCAGAAAAATTTACTGTGAcc                  >  1:1474617/1‑36 (MQ=255)
                     aggcagACATTTGCTGAAAAATTTACTGTGACCGAg               >  1:1664394/4‑36 (MQ=37)
                       tcagACATTTGCAGAAAAATTTACTGTGACCGAGGa             >  1:2577423/1‑36 (MQ=255)
                        cagACATTTGCTGAAAAATTTACTGTGACCGAGGaa            >  1:6982828/1‑36 (MQ=255)
                          gACATTTGCAGAAAAATTTACTGTGACCGAGGAAtt          >  1:5830833/1‑36 (MQ=255)
                          gACATTTGCAGAAAAATTTACTGTGACCGAGGAAtt          <  1:7401577/36‑1 (MQ=255)
                           aCATTTGCAGAAAAATTTACTGTGACCGAGGAATTa         >  1:7425698/1‑36 (MQ=255)
                             aTTTGCAGAAAAATTTACTGTGACCGAGGAATTACt       >  1:3735913/1‑36 (MQ=255)
                             aTTTGCAGAAAAATTTACTGTGACCGAGGAATTACt       <  1:2141074/36‑1 (MQ=255)
                               ttGCAGAAAAATTTACTGTGACCGAGGAATTACTGc     <  1:7659933/36‑1 (MQ=255)
                               ttGCAGAAAAATTTACTGTGACCGAGGAATTACTGc     >  1:1962039/1‑36 (MQ=255)
                                 gCAGAAAAATTTACTGTGACCGAGGAATTACTGCAg   <  1:4453197/36‑1 (MQ=255)
                                  cAGAAAAATTTACTGTGACCGAGGAATTACTGCAGt  <  1:2544380/36‑1 (MQ=255)
                                   |                                  
CAATTATCTCAAACGCAACTCAGTCAGACATTTGCAGAAAAATTTACTGTGACCGAGGAATTACTGCAGT  >  REL606/2329256‑2329325

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG