Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,196 C→T 12.8% A231A (GCC→GCT) ‡ ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,1960CT12.8% 181.9 / 15.5 70A231A (GCC→GCT) ‡ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (2/7);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

TTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATG  >  minE/189126‑189265
                                                                      |                                                                     
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgcc                                                                       >  1:2834480/1‑71 (MQ=255)
ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgcc                                                                       >  1:877985/1‑71 (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)
          aCAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCg                                                             <  1:3200051/71‑1 (MQ=255)
           cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCgg                                                            <  1:951995/71‑1 (MQ=255)
           cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCgg                                                            >  1:3109812/1‑71 (MQ=255)
              agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                         >  1:3779182/1‑71 (MQ=255)
              agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                         >  1:2208361/1‑71 (MQ=255)
              agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                         <  1:2002719/71‑1 (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:2600727/1‑71 (MQ=255)
                     tCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTg                                                  >  1:1964065/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:2743344/71‑1 (MQ=255)
                            cAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCg                                           <  1:3005263/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: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:3354540/71‑1 (MQ=255)
                                               aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                        <  1:3667322/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: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:2441252/71‑1 (MQ=255)
                                                           ggTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCt            <  1:1883741/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:2850883/70‑1 (MQ=255)
                                                               tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt         <  1:213731/70‑1 (MQ=255)
                                                               tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt         >  1:224687/1‑70 (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)
                                                                      |                                                                     
TTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATG  >  minE/189126‑189265

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: