Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,094,789 C→T 7.6% G21G (GGG→GGA)  tdcG ← L‑serine dehydratase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,094,7890CT 7.6% 154.4 / 2.9 53G21G (GGG→GGA) tdcGL‑serine dehydratase 3
Reads supporting (aligned to +/- strand):  ref base C (30/18);  new base T (2/2);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCA  >  minE/2094726‑2094856
                                                               |                                                                   
cgtcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACgg                                                              >  1:2343337/1‑71 (MQ=255)
  tcgCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACg                                                               <  1:3547851/68‑1 (MQ=255)
    gcgGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTAt                                                           <  1:2195380/70‑1 (MQ=255)
          aataaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAAc                                                     <  1:3417285/70‑1 (MQ=255)
          aataaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAAc                                                     <  1:1578825/70‑1 (MQ=255)
           ataaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACt                                                    <  1:1976398/70‑1 (MQ=255)
           ataaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACt                                                    <  1:82307/70‑1 (MQ=255)
            taaGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTg                                                   >  1:3668564/1‑70 (MQ=255)
              aGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                 <  1:3557410/70‑1 (MQ=255)
               gCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAgg                                               <  1:2442444/71‑1 (MQ=255)
                ccGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAggg                                              <  1:164692/71‑1 (MQ=255)
                 cGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGc                                                                     <  1:788024/47‑1 (MQ=255)
                  gCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGAc                                            >  1:1656447/1‑71 (MQ=255)
                   cTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGt                                                             >  1:2865411/1‑53 (MQ=255)
                   cTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATg                                                          <  1:499652/56‑1 (MQ=255)
                   cTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGa                                                 >  1:3069114/1‑65 (MQ=255)
                     aCTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGAc                                            >  1:2220794/1‑68 (MQ=255)
                     aCTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCaa                                         >  1:228068/1‑71 (MQ=255)
                      cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTg                                                   <  1:916623/60‑1 (MQ=255)
                      cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAgg                                               >  1:582941/1‑64 (MQ=255)
                         tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                      >  1:844217/1‑70 (MQ=255)
                          ccAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGc                                                         >  1:3270821/1‑50 (MQ=255)
                            aGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                  <  1:2269486/71‑1 (MQ=255)
                              ccGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAttt                                <  1:2829431/71‑1 (MQ=255)
                                  acAATAAAACTTTTTCCGCCATTCATTGG‑CCCACGGTATGCGAACTg                                                   <  1:3177870/46‑1 (MQ=255)
                                   cAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTAt                                                           >  1:1637726/1‑39 (MQ=255)
                                   cAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCa                                    >  1:2769049/1‑62 (MQ=255)
                                      tAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGaa                                                      <  1:3332576/41‑1 (MQ=255)
                                        aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                      >  1:1161045/1‑71 (MQ=255)
                                         aaCTTTTTCCGGCATTCATTGGCCCCAAGGTATGCGAACTg                                                   <  1:806014/41‑1 (MQ=255)
                                          ttgtttaCCCGCTTTCATTGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCg                     >  1:716263/4‑70 (MQ=255)
                                           cTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                   >  1:1677234/1‑71 (MQ=255)
                                           cTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                   >  1:1424773/1‑71 (MQ=255)
                                            tttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCa                                    <  1:666668/53‑1 (MQ=255)
                                             ttttCCGGCATTCATTGGCCCCACGGTATGCGAACt                                                    >  1:3869745/1‑36 (MQ=255)
                                             ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCa                                    >  1:1399229/1‑52 (MQ=255)
                                             ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:2781429/1‑71 (MQ=255)
                                             ttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCACAATTTTGAAAATATCGAATg                 >  1:2286589/1‑71 (MQ=255)
                                              tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCaa                                   >  1:3692426/1‑52 (MQ=255)
                                              tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:2389924/1‑70 (MQ=255)
                                              tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:691851/1‑70 (MQ=255)
                                                tCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGAcc                                           >  1:1734414/1‑42 (MQ=255)
                                                    gctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAAt          >  1:2976114/4‑71 (MQ=255)
                                                     cattcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                 >  1:2490450/1‑63 (MQ=255)
                                                      tttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGaa                   <  1:3824543/59‑1 (MQ=255)
                                                      attcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGa                             >  1:1959186/1‑50 (MQ=255)
                                                       ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                          >  1:324940/1‑52 (MQ=255)
                                                        tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata      >  1:3131267/1‑71 (MQ=255)
                                                         cattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat     >  1:919355/1‑71 (MQ=255)
                                                          attGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAtac                <  1:2428345/59‑3 (MQ=255)
                                                          attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                     <  1:430453/54‑1 (MQ=255)
                                                           ttGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                         <  1:1707496/49‑1 (MQ=255)
                                                            tGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAata                        >  1:2180987/1‑49 (MQ=255)
                                                            tGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCca  >  1:1812294/1‑71 (MQ=255)
                                                               ccccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                          <  1:1680818/44‑1 (MQ=255)
                                                               |                                                                   
CGTCGCGGTTAATAAGCCGCTACTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCA  >  minE/2094726‑2094856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑