Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,211 T→C 11.9% S236S (AGT→AGC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,2110TC11.9% 176.8 / 12.6 67S236S (AGT→AGCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base C (2/6);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATGTAGTGC  >  minE/189141‑189278
                                                                      |                                                                   
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)
           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:820214/1‑68 (MQ=255)
               aCGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                                         >  1:3720346/1‑68 (MQ=255)
                 ggAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTg                                                               <  1:1628883/60‑1 (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)
                               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: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:3402489/45‑1 (MQ=255)
                                               aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAt                                                <  1:1532994/45‑1 (MQ=255)
                                               aTGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                      >  1:3883740/1‑71 (MQ=255)
                                                tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                      >  1:224687/1‑70 (MQ=255)
                                                tGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGtt                      <  1:213731/70‑1 (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)
                                                      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)
                                                                     gTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATGTAGTGc  >  1:3388078/1‑69 (MQ=255)
                                                                      tACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATGt       <  1:2320742/63‑1 (MQ=255)
                                                                      |                                                                   
GTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCATCTGTTGATGATGAACGATGTAGTGC  >  minE/189141‑189278

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: