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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,2030TG10.4% 178.8 / 17.5 77S234A (TCC→GCC) ‡ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (37/32);  new base G (2/6);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATG  >  minE/189133‑189272
                                                                      |                                                                     
ggAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCAt                                                                       <  1:3144273/71‑1 (MQ=255)
 gAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATc                                                                      >  1:2717643/1‑71 (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:2208361/1‑71 (MQ=255)
       agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                                <  1:2002719/71‑1 (MQ=255)
       agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                                >  1:3779182/1‑71 (MQ=255)
        gTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGt                                                               <  1:3892865/71‑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:3141094/71‑1 (MQ=255)
                 ttACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTg                                                      <  1:1473330/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:820214/1‑68 (MQ=255)
                       aCGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                                   >  1:3720346/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:6910/1‑71 (MQ=255)
                                cGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGt                                       >  1:233731/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:2661353/43‑1 (MQ=255)
                                    tgcGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCAc                                                             <  1:2711827/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:3667322/71‑1 (MQ=255)
                                        aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                               <  1:3354540/71‑1 (MQ=255)
                                        aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                               >  1:632539/1‑71 (MQ=255)
                                         aCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCg                              >  1:1985098/1‑71 (MQ=255)
                                         aCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCg                              >  1:2615817/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:780675/1‑45 (MQ=255)
                                             aaTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGAt                                                    >  1:1246270/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:3402489/45‑1 (MQ=255)
                                                       aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAt                                          <  1:1532994/45‑1 (MQ=255)
                                                       aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                >  1:3883740/1‑71 (MQ=255)
                                                        tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                <  1:2850883/70‑1 (MQ=255)
                                                        tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                >  1:224687/1‑70 (MQ=255)
                                                        tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                <  1:213731/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)
                                                              ccgcGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATa                                   <  1:771326/45‑1 (MQ=255)
                                                              ccgcGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                               <  1:3501113/49‑1 (MQ=255)
                                                               cgcgCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTtgatgatg         >  1:2218557/1‑70 (MQ=255)
                                                                gcgcCATCCGGCAGTACGCTGTTGATCGCCCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGaa       >  1:3147988/1‑71 (MQ=255)
                                                                 cgcCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAAc      <  1:1653186/71‑1 (MQ=255)
                                                                  gcCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAAc      >  1:581932/1‑70 (MQ=255)
                                                                  gcCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACg     <  1:1890757/71‑1 (MQ=255)
                                                                    cATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGAt   >  1:3187404/1‑71 (MQ=255)
                                                                     aTCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATg  >  1:1249521/1‑71 (MQ=255)
                                                                      |                                                                     
GGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATG  >  minE/189133‑189272

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: