Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,363,339 A→G H59H (CAT→CAC yfaU ← predicted 2,4‑dihydroxyhept‑2‑ene‑1,7‑dioic acid aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,363,3390AG83.5% 149.0 / 16.9 79H59H (CAT→CACyfaUpredicted 2,4‑dihydroxyhept‑2‑ene‑1,7‑dioic acid aldolase
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base G (0/66);  total (13/66)
Fisher's exact test for biased strand distribution p-value = 3.79e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCA  >  W3110S.gb/2363306‑2363385
                                 |                                              
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:432262/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:522124/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:518038/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:487090/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:48502/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:454877/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:447510/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:438121/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:435499/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:534110/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:430519/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:420657/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:400999/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:391120/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:342322/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:341384/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:244606/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:209341/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:867510/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:996475/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:981397/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:972190/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:924969/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:913262/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:889635/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:8855/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:86934/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1025814/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:862137/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:841045/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:712959/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:700396/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:661046/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:601891/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:579295/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:551021/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1159889/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:128615/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1261670/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1261622/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1243396/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:124026/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1184099/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1181128/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1172624/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:188378/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1143273/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1135669/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1133425/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1122488/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1112068/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:104513/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1030098/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1293474/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1318900/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1320181/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1335474/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1289543/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:195238/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:188095/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:182271/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:14925/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1453162/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1445388/40‑1 (MQ=255)
gtggctggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1445996/40‑1 (MQ=255)
gtggcgggcATAGGGCGCTACCGCCTGTAGCTGGTGATaaa                                         <  1:1007662/40‑1 (MQ=38)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGATTGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:923329/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGgtgt                          >  1:328984/1‑37 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCt              >  1:1147069/1‑49 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:814569/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:1377857/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:654881/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:370062/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:162437/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:605013/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:1199028/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:600431/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:102868/1‑61 (MQ=255)
                   aCCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCa  >  1:1204758/1‑61 (MQ=255)
                                 |                                              
TTGGCTGGCATAGGGCGCTACCGCCTGTAGCTGATGATAAAGATCCTGAATGGTGTTTGGCGCGTGCTCCCCGTCAATCA  >  W3110S.gb/2363306‑2363385

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: