Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,558,488 | A→G | 32.0% | E225E (GAA→GAG) | gpsA → | glycerol‑3‑phosphate dehydrogenase (NAD+) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,558,488 | 0 | A | G | 32.0% | 63.1 / 39.4 | 78 | E225E (GAA→GAG) | gpsA | glycerol‑3‑phosphate dehydrogenase (NAD+) |
Reads supporting (aligned to +/- strand): ref base A (38/15); new base G (13/12); total (51/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01 |
TCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGG > minE/2558426‑2558547 | tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:2318988/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:973152/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:583065/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:442720/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:386533/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:3304575/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:3155020/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:2286496/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:1882745/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:1815391/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:1723033/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:144655/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:1155232/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:1424741/1‑69 (MQ=255) tCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTcg > 1:1188448/1‑69 (MQ=255) aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc > 1:3070462/1‑69 (MQ=255) aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc > 1:1863774/1‑69 (MQ=255) aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTc > 1:1818865/1‑69 (MQ=255) aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTAGCGTc > 1:358289/1‑69 (MQ=255) aCGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGATGAAATGTCGCGTc > 1:1887143/1‑69 (MQ=255) aTCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGt > 1:2958844/1‑69 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGt > 1:1263267/1‑62 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:301407/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:2315955/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:2370741/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:2528285/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:2646968/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:3002511/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:855101/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:467302/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:789875/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:3045790/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:2538764/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:3243706/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:70263/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:1167557/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:2067837/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:1773304/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:1801138/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:1330267/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:1849929/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:1764783/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:1697462/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg < 1:2090468/69‑3 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:2233943/1‑67 (MQ=255) tCGGTTTTGGCGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAGATGTCGCGTCTTggcg > 1:1678912/1‑67 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:2945863/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:1456102/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:1323265/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:1105895/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:304039/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:678659/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:1897895/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:2858010/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa < 1:1716311/69‑1 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:1699790/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:2578558/1‑69 (MQ=255) tGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGa > 1:662648/1‑69 (MQ=255) gcgcGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGAc < 1:3304314/69‑1 (MQ=255) gcgcGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGAc < 1:1592503/69‑1 (MQ=255) aCGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGcc < 1:1420926/69‑1 (MQ=255) aCGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGcc < 1:2388427/69‑1 (MQ=255) aCGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGcc < 1:2478688/69‑1 (MQ=255) gcgcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCAcc > 1:2659863/1‑69 (MQ=255) gcgcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCAcc > 1:3400125/1‑69 (MQ=255) gcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt < 1:2608601/69‑1 (MQ=255) gcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt < 1:2339717/69‑1 (MQ=255) gcTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt < 1:881709/69‑1 (MQ=255) cTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCtt < 1:225322/68‑1 (MQ=255) cGTGGGCTGGCTTAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATg > 1:3370460/1‑69 (MQ=255) cGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCAt > 1:2648440/1‑68 (MQ=255) cGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATg > 1:1214385/1‑69 (MQ=255) cGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATg > 1:1427865/1‑69 (MQ=255) ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg < 1:1612056/68‑1 (MQ=255) ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg < 1:1587725/68‑1 (MQ=255) ctggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg < 1:2301751/68‑1 (MQ=255) tggctgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg < 1:1406831/67‑1 (MQ=255) tgAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCggg < 1:754775/63‑1 (MQ=255) | TCCGACGGTATCGGTTTTGGTGCGAATGCGCGTACGGCGCTGATCACCCGTGGGCTGGCTGAAATGTCGCGTCTTGGTGCGGCGCTGGGTGCCGACCCTGCCACCTTTATGGGCATGGCGGG > minE/2558426‑2558547 |
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: ‑ |