Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,547,177 G→A 49.1% C104C (TGC→TGT 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,1770GA49.1% ‑1.0 / 52.7 61C104C (TGC→TGTrfaDADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (11/20);  new base A (15/15);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGCCGTCCCA  >  minE/2547112‑2547239
                                                                 |                                                              
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGTCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:2491031/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:3244287/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:2158555/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:3221200/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:1952763/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:1761505/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:1720512/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:165400/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:2424628/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:3140906/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:2940758/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:436099/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:2832122/69‑1 (MQ=255)
tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                                                             <  1:1240843/69‑1 (MQ=255)
 cGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCagt                                                            <  1:2874552/69‑1 (MQ=255)
 cGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCagt                                                            <  1:184404/69‑1 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCt                                                >  1:2390114/1‑66 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:845923/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:1093009/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:2369348/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:2366116/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:2272925/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:2836277/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:2854941/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:193739/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:1877340/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:1758513/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:1591206/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:3385141/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:1451928/1‑67 (MQ=255)
                cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt                                             >  1:1261791/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:1392791/69‑1 (MQ=255)
                                         cAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCcat                    <  1:1827821/69‑1 (MQ=255)
                                             aaCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCata                <  1:2717631/69‑1 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat               >  1:3333919/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat               >  1:2382824/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat               >  1:3206058/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat               >  1:2593705/1‑69 (MQ=255)
                                              aCGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCatat               >  1:2059504/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:2679364/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        <  1:2249835/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        <  1:1675213/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        <  1:1471956/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        >  1:413654/1‑69 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        <  1:1404339/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        <  1:821785/69‑1 (MQ=255)
                                                     ttCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGcc        >  1:1102791/1‑69 (MQ=255)
                                                           cTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGCCGTCCCa  <  1:307755/69‑1 (MQ=255)
                                                           cTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGCCGTCCCa  <  1:775308/69‑1 (MQ=255)
                                                                 |                                                              
TCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCTTGGAGTATTGATAGTTGTTATCCATCATATACTTGCCGTCCCA  >  minE/2547112‑2547239

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: