Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,583,233 T→A 44.8% D54E (GAT→GAA gadW → DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,583,2330TA44.8% 21.5 / 94.6 87D54E (GAT→GAAgadWDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base A (9/30);  total (35/52)
Fisher's exact test for biased strand distribution p-value = 4.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AATTTGGCAGACGATTTTGCGTTTTGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCAATTTAGGGTTA  >  minE/2583177‑2583293
                                                        |                                                            
aaTTTGGCAGTCGATTTTGCGTTTTGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAg                                                    <  1:883725/67‑1 (MQ=255)
           cGATTTTGCGTTTTGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCAttact                                        >  1:4406/1‑68 (MQ=255)
                       ttGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:2752479/63‑1 (MQ=255)
                       ttGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:777230/63‑1 (MQ=255)
                       ttGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:2348350/63‑1 (MQ=255)
                       ttGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:2261980/63‑1 (MQ=255)
                       ttGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:1417860/63‑1 (MQ=255)
                        tGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:475061/62‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:895557/68‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:2950894/68‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:3170881/68‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:1003081/68‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:1647067/68‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:1200908/68‑1 (MQ=255)
                              ccGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:1023857/68‑1 (MQ=255)
                               cGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCt                     <  1:2946895/67‑1 (MQ=255)
                                   tACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACaaaa                                 <  1:1912707/51‑1 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTAc                   >  1:522933/1‑63 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:3123233/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:2904200/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             <  1:1024354/69‑1 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:3007590/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             <  1:2519817/69‑1 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:2477735/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:2366186/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1072578/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:2281343/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:3252510/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:2255862/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             <  1:42414/69‑1 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1997433/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1855374/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:60473/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1378476/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:138697/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:489739/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:721474/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1852881/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1131471/1‑69 (MQ=255)
                                     cgcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAAAATTCGTCAGCTACCACGCa             >  1:475726/1‑69 (MQ=255)
                                      gcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCaa            >  1:584988/1‑69 (MQ=255)
                                      gcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCaa            >  1:866624/1‑69 (MQ=255)
                                      gcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCaa            >  1:1638381/1‑69 (MQ=255)
                                      gcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:1797541/1‑68 (MQ=255)
                                      gcgACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCa             >  1:2757274/1‑68 (MQ=255)
                                             cTGGATATCGAAGAGCTGAGAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1413155/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTatga     <  1:1699301/69‑4 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:3368738/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:3155873/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:66838/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:710543/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:3224348/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:3017990/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:552435/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2424771/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:753898/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2277193/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1224763/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1221059/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1518334/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1639138/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1821443/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1216548/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1914824/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1127931/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1125905/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2108150/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2310360/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2593781/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1068111/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1067834/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1020961/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2897359/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:1011061/69‑2 (MQ=255)
                                             cTGGATATCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAgga     <  1:2945550/69‑2 (MQ=255)
                                                  tatCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCAATTTAGGGTTa  <  1:2910683/67‑1 (MQ=255)
                                                  tatCGATGAGCTAACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCAATTTAGGGTTa  <  1:2943612/67‑1 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:188801/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:1171036/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:393511/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:3413571/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:2609820/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:3145222/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:2774271/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGTTGCCACGCAATTTAGGAttg  >  1:2980762/1‑66 (MQ=255)
                                                  tatCGAAGAGCTGACAGTTTGCCATTACTTACAAAATATTCATCAGTTGCCACGCAATTTAGGAttg  >  1:141004/1‑66 (MQ=255)
                                                   atCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCAATTTAGGGTTa  <  1:1261081/66‑1 (MQ=255)
                                                        |                                                            
AATTTGGCAGACGATTTTGCGTTTTGTTCACCGGATACGCGACGACTGGATATCGATGAGCTGACAGTTTGCCATTACTTACAAAATATTCGTCAGCTACCACGCAATTTAGGGTTA  >  minE/2583177‑2583293

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