Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,547,132 G→A 71.5% T119T (ACC→ACT 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,1320GA71.5% 39.1 / 15.5 42T119T (ACC→ACTrfaDADP‑L‑glycero‑D‑mannoheptose‑6‑epimerase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (3/9);  new base A (15/15);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

ACCGTAGACGTTCAACGGTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCT  >  minE/2547066‑2547196
                                                                  |                                                                
aCCGTAGACGTTCAACGGTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAgg                                                                 >  1:144749/1‑68 (MQ=255)
aCCGTAGACGTTCAACGGTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAgg                                                                 >  1:3002045/1‑68 (MQ=255)
aCCGTAGACGTTCAACGGTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGt                                                                >  1:1999029/1‑69 (MQ=255)
                ggTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCg                                                <  1:1563662/69‑1 (MQ=255)
                ggTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCg                                                <  1:302697/69‑1 (MQ=255)
                                gggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                <  1:1656967/69‑1 (MQ=255)
                                gggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                <  1:2080875/69‑1 (MQ=255)
                                gggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                >  1:1816744/1‑69 (MQ=255)
                                gggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                >  1:1051303/1‑69 (MQ=255)
                                gggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                <  1:1709050/69‑1 (MQ=255)
                                gggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                <  1:2718336/69‑1 (MQ=255)
                                 ggATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGAttt                                <  1:858717/68‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGTCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:2491031/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:3221200/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:2424628/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:2832122/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:2158555/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:2940758/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:1952763/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:3140906/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:1761505/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:3244287/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:1720512/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:165400/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:436099/69‑1 (MQ=255)
                                              tCGGAGGTGCGTCCACCGTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:1240843/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:2272925/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:3385141/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:845923/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:2854941/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:2836277/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:2369348/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:2366116/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:1451928/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:1261791/1‑67 (MQ=255)
                                                              cgTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACAGTAGTGCAGCAGCTctt  >  1:1093009/1‑67 (MQ=255)
                                                               gTAAGTGGCTGCGGAAGAGGCATACAGGAACGGGATTTCGCGCTCCAGACag                  <  1:3080084/52‑1 (MQ=255)
                                                                  |                                                                
ACCGTAGACGTTCAACGGTTTTTCGTACTCGCGGGATTCAATAAAGTCGGAGGTGCGTCCGCCGTAGGTGGCTGCGGAAGAAGCGTACAGGAACGGGATTTCACGCTCCAGGCAGTAGTGCAGCAGCTCCT  >  minE/2547066‑2547196

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: