Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,485,630 A→T 100% I70I (ATT→ATA menF ← isochorismate synthase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,485,6300AT92.5% 119.5 / 1.8 40I70I (ATT→ATAmenFisochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base T (0/37);  total (0/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTA  >  minE/1485625‑1485692
     |                                                              
cGGGTTATCGCACCCAGGACGTCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1155116/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:2081467/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:965008/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:2088394/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:2139934/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:2157436/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:2277956/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:248765/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:342062/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:392262/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:485011/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:548629/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:611720/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:61487/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:692303/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:745211/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:754271/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:881671/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:882197/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1029931/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1943513/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1183968/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1230626/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:127899/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:130347/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1330974/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1337466/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1432514/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:145239/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1537427/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1559485/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1574674/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1575638/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1588538/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1735032/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:1916310/67‑1 (MQ=255)
cGGGTTATCGCACCCAGGACGGCAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGt   <  1:2038173/67‑1 (MQ=255)
 gggTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTa  <  1:509115/67‑1 (MQ=255)
 gggTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTa  <  1:252444/67‑1 (MQ=255)
 gggTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTa  <  1:157717/67‑1 (MQ=255)
     |                                                              
CGGGTAATCGCGCCCAGGACGACAGCTTCTTCATCACCATTACGTTGTTGCCAGTAGAATTGCGGGTA  >  minE/1485625‑1485692

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: