Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,069 G→C 12.2% G522G (GGG→GGC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0690GC12.2% 131.7 / 9.3 49G522G (GGG→GGCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (27/16);  new base C (4/2);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAC  >  minE/190001‑190134
                                                                    |                                                                 
gAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGa                                                                  >  1:550116/1‑70 (MQ=255)
  gCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                               <  1:893249/71‑1 (MQ=255)
     gggCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag                                                            >  1:1280975/1‑71 (MQ=255)
       gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtagta                                                          <  1:2844800/71‑1 (MQ=255)
       gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg                                                               <  1:1913072/66‑1 (MQ=255)
           tatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTagag                                                       >  1:3342698/1‑70 (MQ=255)
             tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta                                                             >  1:995718/1‑62 (MQ=255)
             tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta                                                             >  1:2987016/1‑62 (MQ=255)
             tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                                                     >  1:927767/1‑70 (MQ=255)
                      aggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaaa                                                   >  1:2771031/1‑63 (MQ=255)
                        gggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa                                                    <  1:2156612/60‑1 (MQ=255)
                           aTCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa                                      >  1:3132968/1‑71 (MQ=255)
                             cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta                                                             >  1:3784601/1‑46 (MQ=255)
                             cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                 >  1:3721888/1‑58 (MQ=255)
                             cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc                                    >  1:1751251/1‑71 (MQ=255)
                             cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc                                    >  1:397263/1‑71 (MQ=255)
                             cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                     >  1:3751403/1‑70 (MQ=255)
                              ccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAc                                                  >  1:1664714/1‑56 (MQ=255)
                                ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc                                 <  1:3767392/71‑1 (MQ=255)
                                 gcggagCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgct                                >  1:2416242/1‑71 (MQ=255)
                                      gcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtc                           <  1:2875556/71‑1 (MQ=255)
                                         ggTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                     >  1:3534410/1‑58 (MQ=255)
                                          gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga                                                        <  1:628541/38‑1 (MQ=255)
                                          gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc                         <  1:4049635/69‑1 (MQ=255)
                                          gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                       >  1:2480545/1‑71 (MQ=255)
                                          gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                       <  1:1466335/71‑1 (MQ=255)
                                          gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                       <  1:1680610/71‑1 (MQ=255)
                                          gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                       <  1:2393807/71‑1 (MQ=255)
                                            aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa                                                     >  1:1306369/1‑39 (MQ=255)
                                            aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                       >  1:1335423/1‑69 (MQ=255)
                                             gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                    <  1:738377/71‑1 (MQ=255)
                                                 aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc                >  1:989448/1‑70 (MQ=255)
                                                  ataccTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt                 <  1:1156892/66‑1 (MQ=255)
                                                   ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                 >  1:3688156/1‑36 (MQ=255)
                                                   ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                 >  1:1837572/1‑36 (MQ=255)
                                                    accTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcctg                         <  1:2034573/58‑2 (MQ=19)
                                                     ccTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCa            >  1:2442605/1‑71 (MQ=255)
                                                      cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt           >  1:3512026/1‑71 (MQ=255)
                                                      cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt           >  1:3564507/1‑71 (MQ=255)
                                                       tCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc          >  1:3450749/1‑71 (MQ=255)
                                                         gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTta                                         >  1:2721051/1‑38 (MQ=255)
                                                         gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc          <  1:1918597/69‑1 (MQ=255)
                                                         gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa        >  1:1369840/1‑71 (MQ=255)
                                                            gaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg                    <  1:178726/56‑1 (MQ=255)
                                                            gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct                          >  1:2411829/1‑50 (MQ=16)
                                                            gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct                          >  1:1916936/1‑50 (MQ=16)
                                                               cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc                                     >  1:1728748/1‑36 (MQ=255)
                                                               cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaaa  <  1:2821963/71‑2 (MQ=255)
                                                                gTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt                   <  1:2450675/53‑1 (MQ=255)
                                                                  gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTa                     >  1:244781/1‑49 (MQ=255)
                                                                  gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt       >  1:3267638/1‑63 (MQ=255)
                                                                   ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                 >  1:3199325/1‑52 (MQ=255)
                                                                   ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa   >  1:1260951/1‑66 (MQ=255)
                                                                    |                                                                 
GAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAC  >  minE/190001‑190134

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: