Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,184 G→T 11.6% V227V (GTG→GTT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,1840GT11.6% 239.5 / 9.4 76V227V (GTG→GTTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base T (2/7);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGCGGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCA  >  minE/189114‑189252
                                                                      |                                                                    
tCGCGCGGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATtgtg                                                                      <  1:1837532/71‑1 (MQ=255)
 cgcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                     >  1:3800370/1‑71 (MQ=255)
 cgcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                     >  1:3379104/1‑71 (MQ=255)
 cgcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                     >  1:2504122/1‑71 (MQ=255)
 cgcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                     >  1:1808578/1‑71 (MQ=255)
  gcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTggg                                                                    >  1:2470260/1‑71 (MQ=255)
  gcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                     <  1:2648124/70‑1 (MQ=255)
   cgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGt                                                                   <  1:414653/71‑1 (MQ=255)
    gcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGt                                                                   <  1:3888113/70‑1 (MQ=255)
     cgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTAt                                                                 <  1:3473156/71‑1 (MQ=255)
     cgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTAt                                                                 <  1:1016025/71‑1 (MQ=255)
      ggACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTAt                                                                 <  1:1148431/70‑1 (MQ=255)
      ggACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                <  1:3318694/71‑1 (MQ=255)
       gACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                >  1:2993479/1‑70 (MQ=255)
       gACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                >  1:3458608/1‑70 (MQ=255)
        aCTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTa                                                              <  1:2389298/71‑1 (MQ=255)
            ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgcc                                                          >  1:2834480/1‑71 (MQ=255)
            ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgcc                                                          >  1:877985/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)
                 gcgGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATtgtg                                                                      <  1:1769666/54‑1 (MQ=255)
                  cgGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCa                                                    <  1:775655/71‑1 (MQ=255)
                   ggAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCAt                                                   <  1:3144273/71‑1 (MQ=255)
                    gAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATc                                                  >  1:2717643/1‑71 (MQ=255)
                     aaCAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                <  1:1880775/56‑1 (MQ=255)
                     aaCAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgc                                                           <  1:2429123/61‑1 (MQ=255)
                      aCAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCg                                                <  1:3200051/71‑1 (MQ=255)
                       cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCgg                                               >  1:3109812/1‑71 (MQ=255)
                       cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCgg                                               <  1:951995/71‑1 (MQ=255)
                          agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                            >  1:2208361/1‑71 (MQ=255)
                          agTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAg                                            >  1:3779182/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)
                                 tCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                     <  1:2508023/39‑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)
                                     tACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGt                                                               <  1:2689639/41‑1 (MQ=255)
                                      aCCAACGGAACATCGACGTCGAACAAAATTGTggg                                                                    >  1:3894720/1‑35 (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:930947/1‑71 (MQ=255)
                                                   cGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGt                   >  1:233731/1‑71 (MQ=255)
                                                   cGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGt                   >  1:6910/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)
                                                                      |                                                                    
TCGCGCGGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCA  >  minE/189114‑189252

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: