| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,718,872 | G→C | 9.8% | D118H (GAC→CAC) | hmp → | fused nitric oxide dioxygenase and dihydropteridine reductase 2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,718,872 | 0 | G | C | 9.8% | 141.5 / 4.2 | 51 | D118H (GAC→CAC) | hmp | fused nitric oxide dioxygenase and dihydropteridine reductase 2 |
| Reads supporting (aligned to +/- strand): ref base G (18/27); new base C (2/3); total (21/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGAGGCGGAAATCTATAA > minE/1718805‑1718939 | caacaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAcg > 1:3300930/1‑71 (MQ=255)caacaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAcg > 1:2976003/1‑71 (MQ=255)caacaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGCAct > 1:3327883/1‑70 (MQ=255)caacaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGCAct > 1:949618/1‑70 (MQ=255) acaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGCActtc < 1:2293319/71‑4 (MQ=255) acaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGCActtc < 1:2147347/71‑4 (MQ=255) acaTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGCActtc < 1:4073941/71‑4 (MQ=255) aTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTg < 1:3457214/71‑1 (MQ=255) aTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTg < 1:265469/71‑1 (MQ=255) aTCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTg < 1:2595222/71‑1 (MQ=255) gtcgGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgggg > 1:897664/1‑71 (MQ=255) gtcgGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgggg > 1:3670856/1‑71 (MQ=255) tcgGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgggg > 1:1786796/1‑70 (MQ=255) tcgGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgggg > 1:974721/1‑70 (MQ=255) tcgGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGt < 1:2733259/71‑1 (MQ=255) tcgGTGAACACCTGTTGGCAACGCTGGACGAAATGATCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGt < 1:3090932/71‑1 (MQ=255) ggTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTa < 1:3347858/70‑1 (MQ=255) ggTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTa < 1:4054135/70‑1 (MQ=255) ggTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTa > 1:3771669/1‑70 (MQ=255) aaCACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGcc < 1:3155169/71‑1 (MQ=255) acacCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGcc < 1:1902523/70‑1 (MQ=255) cacCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTa < 1:1333009/71‑1 (MQ=255) cacCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTa < 1:2815218/71‑1 (MQ=255) cacCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTa < 1:1532103/71‑1 (MQ=255) acCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTAt > 1:1575364/1‑71 (MQ=255) tGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaaa > 1:1049577/1‑62 (MQ=255) tGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaaa > 1:2843829/1‑62 (MQ=255) tGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaaa > 1:2050514/1‑62 (MQ=255) tGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGcc < 1:3882725/65‑1 (MQ=255) tGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgt > 1:233395/1‑71 (MQ=255) tGTTGGCAACGCTGGACGAAATGTTAAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaaa > 1:1576834/1‑62 (MQ=255) ttGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgtg < 1:97988/70‑1 (MQ=255) ttGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgtg < 1:2479591/70‑1 (MQ=255) gCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACt < 1:3623645/71‑1 (MQ=255) gCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACt < 1:3711549/71‑1 (MQ=255) gCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACt < 1:609066/71‑1 (MQ=255) gCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACt < 1:4243772/71‑1 (MQ=255) aCGCTGGACGAAATGTTCAGCCCGGTGCAGGAAGTGCTGTACGCGTGGGGTAAAGCCTATGGTGTACTGGc > 1:108162/1‑71 (MQ=255) gCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTa < 1:2798882/71‑1 (MQ=255) aTGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTa > 1:3580839/1‑70 (MQ=255) tGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATc < 1:2169580/71‑1 (MQ=255) cAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATgg > 1:1232070/1‑45 (MQ=255) cAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATgg > 1:744031/1‑45 (MQ=255) gCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATcgc < 1:3337615/71‑1 (MQ=255) gCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATcgc < 1:93103/71‑1 (MQ=255) cccGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATcgc < 1:3307237/70‑1 (MQ=255) gCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTa < 1:716616/56‑1 (MQ=255) gCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTa < 1:815324/56‑1 (MQ=255) gAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGAGGCGGaaa > 1:2073728/1‑70 (MQ=255) tGCTGGTCGCTTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGAGGCGGAAATCtat < 1:3305070/71‑1 (MQ=255) cTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGAGGCGGAAATCTATaa > 1:2933297/1‑71 (MQ=255) gACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATg > 1:3951576/1‑37 (MQ=255) gACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTAttt > 1:1617839/1‑42 (MQ=255) gACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTa > 1:1135472/1‑43 (MQ=255) gACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTa > 1:1130678/1‑43 (MQ=255) | CAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGAGGCGGAAATCTATAA > minE/1718805‑1718939 |
| 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: ‑ |