Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 188,890 A→G 36.7% T129T (ACA→ACG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE188,8900AG36.7% 49.6 / 54.2 87T129T (ACA→ACGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (12/43);  new base G (19/13);  total (31/56)
Fisher's exact test for biased strand distribution p-value = 9.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCA  >  minE/188833‑188958
                                                         |                                                                    
ggAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCa                                                            >  1:927220/1‑68 (MQ=255)
ggAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCa                                                            >  1:818041/1‑68 (MQ=255)
ggAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCa                                                            >  1:1909437/1‑68 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTAcgccgc                                                                  >  1:1546676/1‑56 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:1104483/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:2395121/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:2781501/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:2250722/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:2192234/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:2144776/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:3182545/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:1942835/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:1882989/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:1876913/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:1850552/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:560633/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:964179/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:844118/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:1291254/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:809748/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:796339/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:139033/1‑69 (MQ=255)
      tatCGCCATTCGTATGCGTCAGTACACCAACGAATATCTTGATAACATTACGCCGCCGTTCACGAAAGc                                                     >  1:883199/1‑69 (MQ=255)
        tCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:1320726/69‑1 (MQ=255)
        tCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:222117/69‑1 (MQ=255)
        tCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:2015466/69‑1 (MQ=255)
        tCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:1224822/69‑1 (MQ=255)
        tCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:2994005/69‑1 (MQ=255)
         cGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:1851386/68‑1 (MQ=255)
         cGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:1996966/68‑1 (MQ=255)
         cGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCt                                                   <  1:3020155/68‑1 (MQ=255)
                cGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAg                                                      >  1:41577/1‑58 (MQ=255)
                cGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGttt                                              <  1:1395998/66‑1 (MQ=255)
                                  gacggATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1584742/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGTGGAAg                         <  1:3089724/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGaaa                         <  1:3211101/69‑2 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGaa                          >  1:587225/1‑68 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:45281/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:468919/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2985674/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         >  1:2922191/1‑69 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2785703/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1038586/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:995010/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2678766/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2655830/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1702530/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1163743/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2555389/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2427789/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1218688/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2586348/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1239883/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2258141/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2085798/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:1291266/69‑1 (MQ=255)
                                  gacgaATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAg                         <  1:2349146/69‑1 (MQ=255)
                                        tatCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACAcc                   >  1:336399/1‑69 (MQ=255)
                                        tatCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACAcc                   >  1:2768402/1‑69 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:3402768/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:3044050/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:452497/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:3141855/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:2828846/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:2269397/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:1386023/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:2786609/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:2781479/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:2724244/69‑1 (MQ=255)
                                          tCTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:2591969/69‑1 (MQ=255)
                                          tCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 >  1:1659151/1‑69 (MQ=255)
                                          tCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 >  1:2035047/1‑69 (MQ=255)
                                           cTTGATAACATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCtt                 <  1:3213092/68‑1 (MQ=255)
                                           cTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCttt                >  1:1991481/1‑69 (MQ=255)
                                           cTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCttt                >  1:2079379/1‑69 (MQ=255)
                                                tatcATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:2767627/63‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:3348254/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:1146536/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:3377374/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:352868/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:380399/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:1941338/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:3029377/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:3005703/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:921946/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:2703323/66‑1 (MQ=255)
                                                tAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTg              <  1:2355182/66‑1 (MQ=255)
                                                   cATTACGCCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACTCCtt                 <  1:1196550/60‑1 (MQ=255)
                                                       acacCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCgggg    >  1:3359001/1‑69 (MQ=255)
                                                       acacCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCgggg    <  1:1511835/69‑1 (MQ=255)
                                                       acacCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCgggg    <  1:2115683/69‑1 (MQ=255)
                                                       acacCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCgggg    >  1:3305491/1‑69 (MQ=255)
                                                       acacCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCggg     >  1:1189641/1‑68 (MQ=255)
                                                       acacCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTTACGCCgggg    >  1:2625317/1‑68 (MQ=255)
                                                       acacCGCCGTTCACCAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCgggg    <  1:2815653/69‑1 (MQ=255)
                                                         acCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCa  <  1:1649416/69‑1 (MQ=255)
                                                         |                                                                    
GGAAGATATCGCCATTCGTATGCGTCAGTACACCGACGAATATCTTGATAACATTACACCGCCGTTCACGAAAGCCTTGTTTACCTACGTCAAAGAGCGGAAGTACACCTTTTGTACGCCGGGGCA  >  minE/188833‑188958

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: