Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 188,635 A→G 31.1% K44K (AAA→AAG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE188,6350AG31.1% 81.0 / 45.6 93K44K (AAA→AAGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (27/37);  new base G (21/8);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 7.89e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

GGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGT  >  minE/188569‑188692
                                                                  |                                                         
ggAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAAtt                                                         >  1:3022899/1‑69 (MQ=255)
ggAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAAtt                                                         >  1:2828555/1‑69 (MQ=255)
ggAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAAtt                                                         >  1:3094352/1‑69 (MQ=255)
 gAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAAttt                                                        >  1:3369370/1‑69 (MQ=255)
 gAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAAttt                                                        >  1:1526875/1‑69 (MQ=255)
      ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCg                                                    >  1:552187/1‑68 (MQ=255)
      ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCg                                                    >  1:1811005/1‑68 (MQ=255)
      ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCg                                                    >  1:1582370/1‑68 (MQ=255)
      ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCg                                                    >  1:3347335/1‑68 (MQ=255)
      ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCg                                                    >  1:3216328/1‑68 (MQ=255)
      ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCg                                                    >  1:2793669/1‑68 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:1889507/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:2039703/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:1941827/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:3083029/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:1875729/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:2948363/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:2916404/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:1072745/68‑1 (MQ=255)
        cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAg                                                  <  1:1367761/68‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGct                                                <  1:2727434/69‑2 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:2712476/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:2433688/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:2452841/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:581145/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:1238479/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:2592052/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:519931/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:1908931/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:436094/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:1821772/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:2826549/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:3399653/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:106374/69‑1 (MQ=255)
         gcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTCATTTGCTGAAATTTATCGAGCa                                                <  1:3142176/69‑1 (MQ=255)
          cTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:3061612/68‑1 (MQ=255)
           tggtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:2242488/69‑1 (MQ=255)
           tggtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:2822449/69‑1 (MQ=255)
           tggtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:2560243/69‑1 (MQ=255)
           tggtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:1768143/69‑1 (MQ=255)
           tggtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:2220478/69‑1 (MQ=255)
           tggtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:1689060/69‑1 (MQ=255)
                  gcAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATa                                              <  1:1186026/62‑1 (MQ=255)
                           tCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCa                                                <  1:1808401/51‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:1542919/46‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:1698717/46‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:825973/46‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:738585/46‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:1135678/46‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:3349821/46‑1 (MQ=255)
                             aGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:352688/46‑1 (MQ=255)
                              gATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAGTTTATCGa                                                   <  1:3185214/45‑1 (MQ=255)
                               aTTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTg                          <  1:581533/69‑1 (MQ=255)
                                         cacaAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGa                  <  1:1822839/67‑1 (MQ=255)
                                         cacaAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGa                  <  1:239091/67‑1 (MQ=255)
                                         cacaAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGa                  <  1:315356/67‑1 (MQ=255)
                                               actgCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTg                   >  1:975134/4‑60 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:2834004/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1088439/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:816612/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1117978/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1371669/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:467596/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1628367/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:386395/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:384583/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1695598/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1717284/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:338060/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:1874082/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:3334123/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:2867365/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:2541918/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:2761891/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:3097318/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:2586794/1‑66 (MQ=255)
                                               aCAGCGTTGATTTGCTGAAGTTTATCGAACATAACCCACGAATTTGCGGCGTGATTTTTGACTggg             >  1:2669914/1‑66 (MQ=255)
                                                cAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATg          >  1:66560/1‑68 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTg                   >  1:2058182/1‑52 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTAc     >  1:3338215/1‑66 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:1009462/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:196775/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:971193/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:728656/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:727334/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:675628/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:1201389/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:1224180/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:123173/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:448399/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:2295135/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:211357/1‑69 (MQ=255)
                                                       gATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGt  >  1:2083755/1‑69 (MQ=255)
                                                                  |                                                         
GGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGT  >  minE/188569‑188692

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: