| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 433,649 | T→G | 45.0% | G65G (GGT→GGG) | ribD → | fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 433,649 | 0 | T | G | 45.0% | 0.5 / 2.4 | 21 | G65G (GGT→GGG) | ribD | fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase |
| Reads supporting (aligned to +/- strand): ref base T (4/7); new base G (5/5); total (9/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01 | |||||||||||
CGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATA > NC_000913/433433‑433860 | cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGGG‑ACCGCCTa > 1:440933/1‑228 (MQ=255) gcgtGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGAGGGCGCAAAAGCTGGCCCA‑ACGAGGACGTTTTACCACGCATCCCACCCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGAAAGTT‑TACCACCAACGTGGGGGTGAA‑CCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTca < 1:512035/225‑1 (MQ=255) tGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTAAACGCGTTGCGGATGGCGGGTGGAAACGCCAAAGGTGCG‑CCGCCCCATGTCacac > 1:287682/1‑228 (MQ=255) ggggggggCGC‑AAAGGTGGCGGACAAGAGGTCGTTTTACCACGCTCCCCAACCC‑TAATGGTGGGGGCGTTATTGTCAAAGAGGGCAAATTTGGCGGTGAAGGT‑TACCACCACGGGGCGGGGGAA‑CCACATGCCGAAGTCCCCGCGTTGGGTAGGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcc < 2:418277/223‑1 (MQ=255) cgGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAGCAAAGAGGCG‑ACCGCCCTAGGCACCCCCCAACCCGTGAGCCACAATGGTGCTACGCCACCGTGc > 1:89013/1‑228 (MQ=255) acccacgccaccAACCCGGAATGGCGGGTGCGCAGTTGAAAAAGATGGGAAAATGGTGGGT‑AAGGTGTACCACAAACGGGGGGGAAAC‑CCCAAGGCCGAAGCACGCGGGTTGCGATGGGCGGGGAAAAAAGACAAGGGGGGG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcg < 2:378841/218‑1 (MQ=255) ccACCCC‑GAATGCCGGGTGCTTCATTGTAAAAGGTGGCGAAATTGGCGGTGAAGGT‑TCCCACCAACGTGGGGGTGAA‑CCACAGGCCGAAGTACCCGCGTTGCTTATGGCGGGTAAAAAGGCAAAGGGTGCG‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑t < 1:19008/228‑1 (MQ=255) acatccc‑GGATGTCGGGGGTGTCATTGTCAAAGAGGGCGAAGTGTTCGGGGAGGGT‑CACC‑CCAATGGGGGGGGCAACCCACAGGCGGATATACGCGCGTGGCTTATGGGGGGTAAAAAAGCCAAAGGTGCG‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑t < 2:129885/224‑1 (MQ=255) ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGCGCTGTGACGCACCCCACGCCCCTGGCGTAGCGCGCGCGGTTGCC‑TCGATGCAAGATcc > 1:270211/1‑226 (MQ=255) cAAAGTGGGCGAAATTGGAGGGGAGGGT‑TCCCACCAACGGGGGGGTGAA‑CCACATGCCGAAGTCCCGGCGTTGCGTATGGCGGTGGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGc < 1:195505/228‑1 (MQ=255) cGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCG‑ACCGCCCATGTCACACTCGAACCCTGTAGCCACCATGGTCGTACGCCACCCTGCTGTGACGCACTCATCGCCCCTTGCGTAGCGCGCGTGGTTGCC‑TCGATGCACGATCCTAACCCGCAGTGCCCTGGGC‑ggg > 1:137617/1‑225 (MQ=255) cgggTGACGGT‑TCCCACCAACGTGAGGGTGAA‑CCCCATGCCGAAGTACCCGCGTTGCGTATGGGGGTTAAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTa < 1:416687/224‑1 (MQ=255) caccaACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCCCCCTCGAACCCTGTAACCATCATGGTCGGACGCCACCCTGCTGGGACCCACCCCTCGCCGCCTGCGTGACCCGCGTGCTGGCC‑CCGCTGGCAAATCCCACCCCCCAGCGGCCCGGGC‑GCGGCGTTTACCGTCCCCCACAGGCt > 2:58482/1‑228 (MQ=255) aCGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCA‑ACCGCCTATGTCACAATCGAACCCCGTAACCATCAAGGTCGGACGCCACCGTGCCGTTACGCCCTCATCACCGCTGGCGGAGGGCCCCGGGTTGCCTTCGTGGCAAGTACCAAACCCCCAGG‑GGCGGGGCGGGGGACTTACCCGTCTCCAAAAGGCGGGc > 2:491339/1‑227 (MQ=255) cccACAGGCCGAAGTCCCGGCTGGGCGTGGGGCGGTTGAAAAACCGAAGGGGGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCGTGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa < 2:554376/227‑1 (MQ=255) tGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGGC‑GCCGCCTATGTCACCCTCGAACCCCGTAGCCCACATGGTAGTACGCCCCCGTGCTCTGACGCCCTCATCCCCGCACGCGCATCGCGCGTCGGTGCC‑TCGGCGCAAACTCCTAACCCCCCGGGCGCCGGGC‑GGGGGCTTTTCCGCCTGCAACAGGCTCGCCTTTGCGCAAGCCACGGccg > 2:379341/1‑227 (MQ=255) aaGTACACGTGTTGGGTTGGGCGGGTAAAAAAGCAAAAGGGGCG‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg < 1:436590/227‑1 (MQ=255) gttGGGGGGTGAAAAAGCCCAAGGGTCGCACCGCCTATGTCACACGCGAACCCTGTAGCCAGCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg < 1:210069/225‑1 (MQ=255) aTGGGGGGTGAAAAGCCCAAGGGGGCA‑GCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa < 1:58482/228‑1 (MQ=255) tGGCGGGTGAAAATGCAAAGTGTCCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAACGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt < 1:379341/228‑1 (MQ=255) gTGAAAAAGCCAAAGGTGCG‑ACCCCCTATGTCACACTCGAACCCTGTAGCCATCAAGGTCGTACGCCACCGTGCTGTGACGCACCCATCGCCGCTGGCGGCTCGCCCGGGGTTGCC‑TCGATGCCAGATCCTAACCCGCAGGGCGCCGGGC‑GTGGTCTTTACCGCCCGCACCAGGCCGGCGTTGACGTCAGGCACGGCCCGACGATGAGGGAAACCGAGGAAATGaaaa > 2:335094/1‑226 (MQ=255) | CGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATA > NC_000913/433433‑433860 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |