Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,194 G→T 12.6% A231S (GCC→TCC) ‡ ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,1940GT12.6% 222.4 / 11.0 71A231S (GCC→TCC) ‡ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base T (2/7);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGA  >  minE/189126‑189263
                                                                    |                                                                     
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgcc                                                                     >  1:877985/1‑71 (MQ=255)
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgcc                                                                     >  1:2834480/1‑71 (MQ=255)
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgc                                                                      >  1:2669286/1‑70 (MQ=255)
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgc                                                                      >  1:400174/1‑70 (MQ=255)
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgc                                                                      <  1:2147232/70‑1 (MQ=255)
  ggCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgccg                                                                    <  1:941026/70‑1 (MQ=255)
   gcgcGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCcgcg                                                                  >  1:1330996/1‑71 (MQ=255)
      cgGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCa                                                               <  1:775655/71‑1 (MQ=255)
       ggAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCAt                                                              <  1:3144273/71‑1 (MQ=255)
        gAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATc                                                             >  1:2717643/1‑71 (MQ=255)
         aaCAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgc                                                                      <  1:2429123/61‑1 (MQ=255)
          aCAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCg                                                           <  1:3200051/71‑1 (MQ=255)
           cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCgg                                                          <  1:951995/71‑1 (MQ=255)
           cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCgg                                                          >  1:3109812/1‑71 (MQ=255)
              agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                       <  1:2002719/71‑1 (MQ=255)
              agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                       >  1:2208361/1‑71 (MQ=255)
              agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                       >  1:3779182/1‑71 (MQ=255)
               gTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGt                                                      <  1:3892865/71‑1 (MQ=255)
                  atatCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGcc                                                                <  1:3010056/58‑1 (MQ=255)
                   tatCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGc                                                  <  1:1843947/71‑1 (MQ=255)
                     tCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTg                                                >  1:1964065/1‑71 (MQ=255)
                     tCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTg                                                >  1:2600727/1‑71 (MQ=255)
                        ttACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTg                                             <  1:1473330/71‑1 (MQ=255)
                        ttACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTg                                             <  1:3141094/71‑1 (MQ=255)
                        tgaCCAACGGTACTTCCACTGCGAACAAAATTGTTGGTATGTACTCTGCTCCAgcaggcag                                                       <  1:542284/59‑1 (MQ=255)
                          aCCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGAt                                           >  1:690367/1‑71 (MQ=255)
                           ccAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGa                                            <  1:235535/69‑1 (MQ=255)
                            cAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCg                                         <  1:3005263/71‑1 (MQ=255)
                            cAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCg                                         <  1:2743344/71‑1 (MQ=255)
                              aCGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                          >  1:3720346/1‑68 (MQ=255)
                              aCGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                          >  1:820214/1‑68 (MQ=255)
                                ggAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTg                                                <  1:1628883/60‑1 (MQ=255)
                                    cATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCa                                                        >  1:3336808/1‑48 (MQ=255)
                                    cATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCa                                                        >  1:642726/1‑48 (MQ=255)
                                       cGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGt                              >  1:233731/1‑71 (MQ=255)
                                       cGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGt                              >  1:6910/1‑71 (MQ=255)
                                       cGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGt                              >  1:930947/1‑71 (MQ=255)
                                        gACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTg                                             <  1:1803808/55‑1 (MQ=255)
                                         aCGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCa                            >  1:1704588/1‑71 (MQ=255)
                                          ctgcGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTaa                                  <  1:248430/61‑1 (MQ=255)
                                          cGTCGAACAAAATCGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGa                                        >  1:1360648/1‑58 (MQ=255)
                                           tgcGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCAc                                                    <  1:2711827/43‑1 (MQ=255)
                                           tgcGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCAc                                                    <  1:2661353/43‑1 (MQ=255)
                                           gTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCaa                                  <  1:3052458/63‑1 (MQ=255)
                                            tCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGc                                                  <  1:3079407/46‑1 (MQ=255)
                                            tCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATaa                         <  1:3226810/71‑1 (MQ=255)
                                            gcGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTaa                                  >  1:3658149/2‑62 (MQ=255)
                                             cGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCACaaa                        >  1:1595306/1‑71 (MQ=255)
                                             cGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCACaa                         <  1:958532/70‑1 (MQ=255)
                                             cGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCACaa                         <  1:1529292/70‑1 (MQ=255)
                                              gAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCACAAAt                       <  1:1642960/71‑1 (MQ=255)
                                              gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCa                                                        >  1:3643118/1‑38 (MQ=255)
                                              gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATaaa                        <  1:3724611/70‑1 (MQ=255)
                                              gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAAt                       >  1:3053044/1‑71 (MQ=255)
                                               aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                      <  1:3354540/71‑1 (MQ=255)
                                               aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                      >  1:632539/1‑71 (MQ=255)
                                               aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                      <  1:3667322/71‑1 (MQ=255)
                                                aCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCg                     >  1:2615817/1‑71 (MQ=255)
                                                aCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCg                     >  1:1985098/1‑71 (MQ=255)
                                                 cAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGc                    >  1:1673698/1‑71 (MQ=255)
                                                   aaaTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTg                  >  1:3691979/1‑71 (MQ=255)
                                                    aaTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGc                                                  <  1:1870698/38‑1 (MQ=255)
                                                    aaTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGAt                                           >  1:1246270/1‑45 (MQ=255)
                                                    aaTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGAt                                           >  1:780675/1‑45 (MQ=255)
                                                     aTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAATGTCATAAATCGCTgg                 <  1:845810/70‑1 (MQ=255)
                                                       tgtgGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGgcg               >  1:1209976/1‑70 (MQ=255)
                                                        gtgGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCa             >  1:146827/1‑71 (MQ=255)
                                                           ggTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCt          <  1:1883741/71‑1 (MQ=255)
                                                           ggTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCt          <  1:2441252/71‑1 (MQ=255)
                                                            gTATGTACGCCGCGCCATCCGGGAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCt          <  1:3763025/70‑1 (MQ=255)
                                                            gTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCaa                                  >  1:715104/1‑46 (MQ=255)
                                                              aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAt                                 <  1:1532994/45‑1 (MQ=255)
                                                              aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAt                                 <  1:3402489/45‑1 (MQ=255)
                                                              aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt       >  1:3883740/1‑71 (MQ=255)
                                                               tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt       <  1:213731/70‑1 (MQ=255)
                                                               tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt       >  1:224687/1‑70 (MQ=255)
                                                               tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt       <  1:2850883/70‑1 (MQ=255)
                                                                 tACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTtgat    >  1:3370770/1‑71 (MQ=255)
                                                                   cgccgcGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTtgatga  >  1:3658826/1‑71 (MQ=255)
                                                                   cgccgcGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTtgatg   <  1:2470767/70‑1 (MQ=255)
                                                                    |                                                                     
TTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGA  >  minE/189126‑189263

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