| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,685,303 | T→C | 18.2% | Y42H (TAT→CAT) | sseA → | 3‑mercaptopyruvate sulfurtransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,685,303 | 0 | T | C | 18.2% | 97.7 / 16.5 | 44 | Y42H (TAT→CAT) | sseA | 3‑mercaptopyruvate sulfurtransferase |
| Reads supporting (aligned to +/- strand): ref base T (11/25); new base C (5/3); total (16/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.17e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
CCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCC > minE/1685233‑1685370 | cccGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGt < 1:1886231/71‑1 (MQ=255) ccGGAAATTCAGATTATCGATGCCGGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTa < 1:3131137/71‑1 (MQ=255) cGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGCAt < 1:1962992/71‑1 (MQ=255) ggAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATc > 1:2150502/1‑71 (MQ=255) cAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATg < 1:261106/70‑1 (MQ=255) aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATgg < 1:529065/70‑1 (MQ=255) aGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATgg < 1:2809994/70‑1 (MQ=255) tCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACatat < 1:2548652/70‑1 (MQ=255) cccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGgc > 1:1338797/1‑71 (MQ=255) ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt < 1:1158583/49‑1 (MQ=255) ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt < 1:459778/49‑1 (MQ=255) ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt < 1:247374/49‑1 (MQ=255) ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt < 1:2402950/49‑1 (MQ=255) ccGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTAt < 1:2083697/49‑1 (MQ=255) cATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCt < 1:3140603/48‑1 (MQ=255) cGTCGCCTGGACAGGAGGAGCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGttt < 1:1708135/71‑1 (MQ=255) tCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGttttt < 1:913630/71‑1 (MQ=255) gCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGtttttt < 1:2201177/70‑1 (MQ=255) ccTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACata > 1:2257445/1‑49 (MQ=255) ccTGGACAGGAGGATCGTAACGTTGCTCAGGAGCATCTGAATGGACATATTcc < 1:1874599/53‑1 (MQ=255) ggACAGGAGGATCGTAACGTTGCTCAGGAGCATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGatat > 1:263942/1‑71 (MQ=255) aCAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATc > 1:1242493/1‑70 (MQ=255) aCAGGAGGATCGTAACGTTGCTCAGGAGCATCTGAATGGACat < 1:3139981/43‑1 (MQ=255) aggaggATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCg > 1:2570978/1‑69 (MQ=255) aggaggATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCg > 1:588821/1‑69 (MQ=255) aggaggATCGTAACGTTGCTCAGGAGCATCTGAATGGACatat > 1:1127061/1‑43 (MQ=255) aggaggATCGTAACGTTGCTCAGGAGCATCTGAATGGACatat > 1:2501366/1‑43 (MQ=255) gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc < 1:1803145/71‑1 (MQ=255) gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc < 1:1798885/71‑1 (MQ=255) gaggaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAgc > 1:2529086/1‑71 (MQ=255) gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTc < 1:965393/40‑1 (MQ=255) gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTc < 1:2389655/40‑1 (MQ=255) gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTc < 1:3705291/40‑1 (MQ=255) gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCt < 1:3237934/71‑1 (MQ=255) gaTCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCt < 1:775032/71‑1 (MQ=255) gTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCt > 1:2997785/1‑71 (MQ=255) aCGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGAt < 1:2518296/71‑1 (MQ=255) aCGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGAt > 1:577740/1‑71 (MQ=255) ttGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCg < 1:2829671/54‑1 (MQ=255) cTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCacac < 1:3389062/70‑1 (MQ=255) cTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCacac < 1:3497530/70‑1 (MQ=255) aGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAg > 1:527517/1‑51 (MQ=255) aGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcc > 1:1205152/1‑71 (MQ=255) gAGCATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcccc > 1:3585313/1‑71 (MQ=255) gAGCATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTcccc > 1:2827037/1‑71 (MQ=255) tATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATc < 1:3158695/42‑1 (MQ=255) tATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCtt < 1:567408/55‑2 (MQ=255) | CCCGGAAATTCAGATTATCGATGCCCGCATGGCGTCGCCTGGACAGGAGGATCGTAACGTTGCTCAGGAGTATCTGAATGGACATATTCCCGGCGCAGTGTTTTTTGATATCGAAGCGCTTTCTGATCACACTTCCCC > minE/1685233‑1685370 |
| 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: ‑ |