Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,205 C→A 10.4% S234S (TCC→TCA) ‡ ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,2050CA10.4% 239.3 / 8.6 76S234S (TCC→TCA) ‡ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base A (2/6);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATG  >  minE/189136‑189272
                                                                     |                                                                   
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:2002719/71‑1 (MQ=255)
    agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                                                >  1:2208361/1‑71 (MQ=255)
     gTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGt                                                               <  1:3892865/71‑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: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: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:632539/1‑71 (MQ=255)
                                     aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                               <  1:3667322/71‑1 (MQ=255)
                                     aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc                               <  1:3354540/71‑1 (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)
                                                                     |                                                                   
ACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATG  >  minE/189136‑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: