Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,475,665 | A→T | 34.4% | T244S (ACC→TCC) | nrdB → | ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,475,665 | 0 | A | T | 34.4% | 59.9 / 51.9 | 58 | T244S (ACC→TCC) | nrdB | ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like |
Reads supporting (aligned to +/- strand): ref base A (29/9); new base T (9/11); total (38/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01 |
AGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAT > minE/1475595‑1475729 | aGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCAGCGACGAAGCCCTGCACCTGa < 1:1696308/71‑1 (MQ=255) aaTTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCAcc > 1:4050549/1‑71 (MQ=255) tGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCa < 1:1933643/70‑1 (MQ=255) gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCatat > 1:2415058/1‑71 (MQ=255) gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCatat > 1:4019687/1‑71 (MQ=255) gAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCatat > 1:1856817/1‑71 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCTGCACCCAGCATATGCTGa < 1:3435938/71‑1 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata > 1:27364/1‑65 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCata > 1:2409515/1‑65 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGa < 1:2581722/71‑1 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGa < 1:2782624/71‑1 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGa < 1:2153528/71‑1 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGc > 1:2760631/1‑68 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGc > 1:2381027/1‑68 (MQ=255) cAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGc > 1:3496697/1‑68 (MQ=255) aaCGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGaa > 1:479337/1‑71 (MQ=255) aaCGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGaa > 1:544233/1‑71 (MQ=255) cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAAt > 1:3446092/1‑70 (MQ=255) cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATc > 1:2997529/1‑71 (MQ=255) cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATc > 1:2786424/1‑71 (MQ=255) cGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATc > 1:3876872/1‑71 (MQ=255) gCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATct > 1:3116808/1‑71 (MQ=255) ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAAt > 1:3847034/1‑68 (MQ=255) ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAAt > 1:854641/1‑68 (MQ=255) ccAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAAt > 1:903938/1‑68 (MQ=255) attattCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTgcg > 1:3695916/1‑71 (MQ=255) ttattCGCCTGATTGCCAGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctgc > 1:2123061/1‑66 (MQ=255) tCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAgcg > 1:136550/1‑71 (MQ=255) ccTGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAgcggcg < 1:2991391/71‑1 (MQ=255) tGATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg > 1:3193270/1‑71 (MQ=255) tGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTg < 1:697335/50‑1 (MQ=255) tGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTg < 1:3276048/50‑1 (MQ=255) gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg > 1:3034946/1‑70 (MQ=255) gATTGCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg > 1:2236627/1‑70 (MQ=255) aTTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg < 1:1052986/69‑1 (MQ=255) aTTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg < 1:1296140/69‑1 (MQ=255) aTTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGcgcg < 1:328083/69‑1 (MQ=255) gCCCGCGACGAAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCCGCTGCGCAGCGGCGCGGACg < 1:2263600/70‑1 (MQ=255) cccGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGa < 1:3145287/70‑1 (MQ=255) gacgaAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga < 1:785782/71‑1 (MQ=255) gacgaAGCCCTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga < 1:781035/71‑1 (MQ=255) gacgaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTga < 1:1982986/71‑1 (MQ=255) gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt > 1:3000900/1‑71 (MQ=255) gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt > 1:3113611/1‑71 (MQ=255) gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt > 1:3783330/1‑71 (MQ=255) gaAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGAt > 1:3561668/1‑71 (MQ=255) aGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATctgc > 1:1085509/1‑41 (MQ=255) aGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATga > 1:3807341/1‑70 (MQ=255) cTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa < 1:3885098/71‑1 (MQ=255) cTGCACCTGTCCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa < 1:1212125/71‑1 (MQ=255) cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa < 1:1929052/71‑1 (MQ=255) cTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGa < 1:3312226/71‑1 (MQ=255) gCACCTGACCGGCACCCAGCATATGCTGAATCTTCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa > 1:2761962/1‑71 (MQ=255) gCACATGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGAGCCTGAGATGGCGGaaa > 1:1708579/1‑71 (MQ=255) cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGaaa > 1:670780/1‑70 (MQ=255) cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt > 1:3538549/1‑71 (MQ=255) cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt > 1:3282385/1‑71 (MQ=255) cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt > 1:897421/1‑71 (MQ=255) cACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAt > 1:3325665/1‑71 (MQ=255) | AGAACGCGAATTGATGGAAGGCAACGCCAAAATTATTCGCCTGATTGCCCGCGACGAAGCCCTGCACCTGACCGGCACCCAGCATATGCTGAATCTGCTGCGCAGCGGCGCGGACGATCCTGAGATGGCGGAAAT > minE/1475595‑1475729 |
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: ‑ |