Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,184 G→T 9.2% 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,1840GT 9.2% 209.6 / 5.0 65V227V (GTG→GTTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base T (2/4);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CGCGCGGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCAT  >  minE/189115‑189253
                                                                     |                                                                     
cgcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                      >  1:2237304/1‑71 (MQ=255)
 gcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTggg                                                                     <  1:936190/71‑1 (MQ=255)
 gcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTggg                                                                     <  1:351132/71‑1 (MQ=255)
 gcgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTggg                                                                     >  1:369627/1‑71 (MQ=255)
  cgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGt                                                                    <  1:3860261/71‑1 (MQ=255)
  cgcgGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGt                                                                    <  1:1815588/71‑1 (MQ=255)
     ggACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                 <  1:3694745/71‑1 (MQ=255)
        cTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAc                                                              >  1:1635857/1‑71 (MQ=255)
         ttttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcg                                                             >  1:427410/1‑71 (MQ=255)
         ttttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcg                                                             >  1:358159/1‑71 (MQ=255)
           ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCTACGTCGAACAAAATTGTGGGTATGTAcgcc                                                           >  1:3890852/1‑71 (MQ=255)
           ttGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgc                                                            >  1:1267283/1‑70 (MQ=255)
              gcgcGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgccgc                                                         <  1:672161/70‑1 (MQ=255)
              gcgcGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgccgc                                                         <  1:2816228/70‑1 (MQ=255)
                 cgGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                 <  1:867325/59‑1 (MQ=255)
                   gAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                 <  1:2392419/57‑1 (MQ=255)
                   gAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                 <  1:2837987/57‑1 (MQ=255)
                   gAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                 <  1:1445493/57‑1 (MQ=255)
                      cAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCa                                                     >  1:2087194/1‑66 (MQ=255)
                          gTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCAt                                                    >  1:767364/1‑63 (MQ=255)
                          gTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCAt                                                    >  1:706291/1‑63 (MQ=255)
                          gTTATATCGTTACCAACGGAACATCGACGTCGAACAAAAT‑‑TGGGTATGTACGCCGCGCCAt                                                    >  1:2869515/1‑61 (MQ=255)
                           ttATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                      <  1:2498662/44‑1 (MQ=255)
                           ttATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTa                                                                   <  1:4077573/47‑1 (MQ=255)
                            tataTCGTTACCAACGGAACATCGACGTCGAACAAAATtgtg                                                                       >  1:4060671/1‑42 (MQ=255)
                             atatCGTTACCAACGGAACATCGACGTCGAACAAAATTGTgg                                                                      >  1:747902/1‑42 (MQ=255)
                             atatCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTAcgccgc                                                         >  1:3692682/1‑55 (MQ=255)
                             atatCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACg                                         >  1:3513510/1‑71 (MQ=255)
                              tatCGTTACCAACGTAACATCGACGTCGAACAAAATTGTGGGTa                                                                   <  1:3219928/44‑1 (MQ=255)
                              tatCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGt                                                                    <  1:2192168/43‑1 (MQ=255)
                               atCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATg                                                                 <  1:401639/45‑1 (MQ=255)
                                 cGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGt                                     >  1:2555742/1‑71 (MQ=255)
                                 cGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGt                                     >  1:653708/1‑71 (MQ=255)
                                      ccAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                <  1:2727966/71‑1 (MQ=255)
                                      ccAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                <  1:1946526/71‑1 (MQ=255)
                                       cAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCg                               <  1:164692/71‑1 (MQ=255)
                                        aaCGGTACTTCCACTGCGAACAAAATTGTTGGTATGTACTCTGCTCCag                                                    <  1:1047279/49‑2 (MQ=255)
                                         aCGGAACATCGACGTCGAACAAAATTGTGGGTATGt                                                                >  1:395519/1‑36 (MQ=255)
                                            gAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCg                           >  1:1014092/1‑70 (MQ=255)
                                             aaCATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGc                          >  1:1378356/1‑70 (MQ=255)
                                                aTCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTg                                   <  1:1098583/58‑1 (MQ=255)
                                                 tCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCt                                       >  1:3628666/1‑53 (MQ=255)
                                                 tCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGa                                  <  1:1796754/58‑1 (MQ=255)
                                                     cGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATc                                                   >  1:1086150/1‑37 (MQ=255)
                                                     cGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCaa                        >  1:1885878/1‑64 (MQ=255)
                                                     cGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAt                       >  1:3197765/1‑65 (MQ=255)
                                                     cGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAAt                       >  1:797214/1‑65 (MQ=255)
                                                     cGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCa                  >  1:3904913/1‑70 (MQ=255)
                                                      gTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCAt                 >  1:3034066/1‑70 (MQ=255)
                                                       tCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCa                                              <  1:2037503/40‑1 (MQ=255)
                                                       tCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATaa               <  1:2770815/71‑1 (MQ=255)
                                                       gcGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTaac                       >  1:2405916/2‑62 (MQ=255)
                                                        cGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCaca                <  1:3512941/69‑3 (MQ=255)
                                                        cGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCaca                <  1:273590/69‑3 (MQ=255)
                                                        cGAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCACaaa              >  1:1401832/1‑71 (MQ=255)
                                                        cGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCa                                              >  1:3220448/1‑39 (MQ=255)
                                                        cGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATaaa              >  1:3138174/1‑71 (MQ=255)
                                                         gAACAAAATTGTTGGTATGTACTCTGCTCCAGCAGGCAGCACCATTCTGATTGACCGTAACTGCCACAAAt             <  1:2307331/71‑1 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGt                                            >  1:2458660/1‑40 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTa                                           >  1:2325206/1‑41 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATc                                <  1:2889164/52‑1 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAAt             >  1:1743851/1‑71 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAAt             >  1:1477894/1‑71 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAAt             >  1:1496904/1‑71 (MQ=255)
                                                         gAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAAt             >  1:3339934/1‑71 (MQ=255)
                                                          aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAAt             <  1:3273680/70‑1 (MQ=255)
                                                          aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc            <  1:496206/71‑1 (MQ=255)
                                                          aaCAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATc            <  1:3982843/71‑1 (MQ=255)
                                                              aaaTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCg                               <  1:3631015/48‑1 (MQ=255)
                                                                aTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTgg       <  1:1411818/70‑1 (MQ=255)
                                                                aTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGgc      >  1:2122675/1‑71 (MQ=255)
                                                                aTTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGgc      <  1:2569125/71‑1 (MQ=255)
                                                                    tgGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCAt  >  1:1077841/1‑71 (MQ=255)
                                                                     |                                                                     
CGCGCGGACTTTTGGCGCGGAACAGAGTTATATCGTTACCAACGGAACATCGACGTCGAACAAAATTGTGGGTATGTACGCCGCGCCATCCGGCAGTACGCTGTTGATCGACCGCAATTGTCATAAATCGCTGGCGCAT  >  minE/189115‑189253

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: