Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,547,168 A→G 51.8% R107R (CGT→CGC rfaD ← ADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,547,1680AG51.8% 0.9 / 49.0 58R107R (CGT→CGCrfaDADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (11/17);  new base G (15/15);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGCC  >  minE/2547112‑2547233
                                                        |                                                                 
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGTCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:2491031/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:1952763/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:436099/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:1240843/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:2158555/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:3244287/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:3221200/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:3140906/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:2940758/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:1761505/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:165400/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:2424628/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:1720512/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:2832122/69‑1 (MQ=255)
 cGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCagt                                                      <  1:184404/69‑1 (MQ=255)
 cGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCagt                                                      <  1:2874552/69‑1 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCt                                          >  1:2390114/1‑66 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:2366116/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:1093009/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:2836277/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:2272925/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:2369348/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:2854941/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:3385141/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:845923/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:1261791/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:1451928/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:193739/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:1877340/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:1591206/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                       >  1:1758513/1‑67 (MQ=255)
                 gTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                       <  1:3080084/52‑1 (MQ=255)
                      tGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAg                                  <  1:1990076/68‑1 (MQ=255)
                      tGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAg                                  <  1:633679/68‑1 (MQ=255)
                      tGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAg                                  <  1:2836113/68‑1 (MQ=255)
                                  aagCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAgttgt                     <  1:2225749/69‑1 (MQ=255)
                                  aagCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAgttgt                     <  1:2366955/69‑1 (MQ=255)
                                   agCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGTTAgttgtt                    >  1:2856605/1‑69 (MQ=255)
                                         cAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCcat              <  1:2745271/69‑1 (MQ=255)
                                         cAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCcat              <  1:1827821/69‑1 (MQ=255)
                                         cAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCcat              <  1:1392791/69‑1 (MQ=255)
                                             aaCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCata          <  1:2717631/69‑1 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat         >  1:2593705/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat         >  1:3333919/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat         >  1:2382824/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat         >  1:2059504/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat         >  1:3206058/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCata          >  1:246756/1‑68 (MQ=255)
                                              aCGGGAATTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat         >  1:1498/1‑69 (MQ=255)
                                                     ttaaCGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:2458903/66‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGTGTATTGATAGTTGTTATCCATCATATACTTGcc  >  1:123664/1‑69 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  >  1:1102791/1‑69 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:821785/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:2249835/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  >  1:413654/1‑69 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:1404339/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:1471956/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:1675213/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc  <  1:2679364/69‑1 (MQ=255)
                                                        |                                                                 
TCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGCC  >  minE/2547112‑2547233

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: