Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,285,082 T→C 33.6% P152P (CCA→CCG argH ← argininosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,285,0820TC33.6% 67.9 / 49.8 92P152P (CCA→CCGargHargininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (24/37);  new base C (16/15);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGC  >  minE/2285028‑2285147
                                                      |                                                                 
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:62728/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:1097154/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:14965/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:1142129/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2308348/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:1980164/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2214221/69‑1 (MQ=255)
cAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2118923/69‑1 (MQ=255)
 aGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:2998934/68‑1 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1693620/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:97007/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1727930/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1949567/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1953782/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2035367/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2060026/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2269145/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2399520/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2498868/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2569526/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:2584852/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:3320688/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:1315801/1‑69 (MQ=255)
  ggCACCAGTGAGCGAACGTCACAGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGCATCACCGCGTCCt                                                   >  1:3283648/1‑69 (MQ=255)
  ggCACCAGTGAGAGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCCGGGATCACCGCGTCCt                                                   >  1:1277132/1‑69 (MQ=255)
     aCCAGTGCGCGAACGTCCCCGGCTGGGCGCGTTGCTGGTGAGTGTACCCTGGCATTACCGCGTc                                                     >  1:2058512/1‑64 (MQ=255)
     aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:1581728/1‑69 (MQ=255)
     aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:1116491/1‑69 (MQ=255)
     aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:2411294/1‑69 (MQ=255)
     aCCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGAt                                                >  1:1779627/1‑69 (MQ=255)
        aGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTc                                                     <  1:198879/61‑1 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:101529/1‑69 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:105507/1‑69 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:361989/1‑69 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:2816863/1‑69 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:2812569/1‑69 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:1140845/1‑69 (MQ=255)
           gcgcGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGtttt                                          >  1:1923777/1‑69 (MQ=255)
            cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:1570225/1‑69 (MQ=255)
            cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:3234703/1‑69 (MQ=255)
            cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:923778/1‑69 (MQ=255)
            cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:3321159/1‑69 (MQ=255)
            cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:1535560/1‑69 (MQ=255)
            cgcgAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:1140033/1‑69 (MQ=255)
            cgcgAACGTCACAGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTtg                                         >  1:2991061/1‑69 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTTTTTTGTGCg                                     <  1:2711000/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:50722/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1061800/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:3383346/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2740772/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1533822/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2381699/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1368591/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1605184/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2270039/69‑1 (MQ=255)
                aaCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:2208342/69‑1 (MQ=255)
                 aCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:1708219/68‑1 (MQ=255)
                 aCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:691397/68‑1 (MQ=255)
                 aCGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCg                                     <  1:951235/68‑1 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:3278141/69‑2 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:3230560/69‑2 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:330654/69‑2 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2044930/69‑2 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2073026/69‑2 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1615582/69‑2 (MQ=255)
                                     tGCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2285083/69‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCa               >  1:2875978/1‑69 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCTGTCTCGACCAGAGCGGATTgt                <  1:1679960/68‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGAtttt                <  1:1737729/68‑3 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:711793/68‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1413272/68‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:578272/68‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1469214/68‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:1355052/68‑2 (MQ=255)
                                      gCAGGTGAGTGTAACCCGGCATCACCGCGTCCTGATTGTTCTGCGCGGTCTCGACCAGAGCGGATTgt                <  1:2581416/68‑2 (MQ=255)
                                            gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc         >  1:3129575/1‑69 (MQ=255)
                                            gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc         >  1:1104946/1‑69 (MQ=255)
                                            gAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGcc         >  1:1221619/1‑69 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:42006/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1401683/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1318674/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1902326/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1150550/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:2449793/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:455884/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:3273357/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:2959996/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:33434/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:2981337/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:3299290/51‑1 (MQ=255)
                                                 tAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGcc                      <  1:1622076/51‑1 (MQ=255)
                                                   aCCTGTCATTACAGCGTCCTGATTGTTTTGTTCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGc  >  1:3431672/1‑69 (MQ=255)
                                                      |                                                                 
CAGGCACCAGTGCGCGAACGTCACCGGCTGGGCGCGTTGCAGGTGAGTGTAACCTGGCATTACCGCGTCCTGATTGTTTTGTGCGGTTTCCACCAGCGCCGATTGCAGCTGCCGGTTAGC  >  minE/2285028‑2285147

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: