Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 190,078 A→T 17.0% V525V (GTA→GTT ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE190,0780AT17.0% 114.9 / 17.4 47V525V (GTA→GTTldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (26/13);  new base T (6/2);  total (32/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

GCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAA  >  minE/190008‑190142
                                                                      |                                                                
gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtagta                                                                  <  1:2844800/71‑1 (MQ=255)
    tatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTagag                                                               >  1:3342698/1‑70 (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)
                      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:1680610/71‑1 (MQ=255)
                                   gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               <  1:1466335/71‑1 (MQ=255)
                                   gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               <  1:2393807/71‑1 (MQ=255)
                                   gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt                               >  1:2480545/1‑71 (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:1837572/1‑36 (MQ=255)
                                            ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc                                                         >  1:3688156/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)
                                                           ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa        >  1:3703444/4‑70 (MQ=255)
                                                           ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa        >  1:1480316/4‑70 (MQ=255)
                                                            ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt                         >  1:3199325/1‑52 (MQ=255)
                                                            ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa           >  1:1260951/1‑66 (MQ=255)
                                                               tCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCa   <  1:3279586/71‑1 (MQ=255)
                                                                cGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa  >  1:4054129/1‑71 (MQ=255)
                                                                  tagtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt                           <  1:3837901/44‑1 (MQ=255)
                                                                   agtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt               >  1:4042529/1‑55 (MQ=255)
                                                                      |                                                                
GCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAA  >  minE/190008‑190142

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: