| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 433,649 | T→G | 45.1% | 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.1% | 5.0 / 3.8 | 27 | G65G (GGT→GGG) | ribD | fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase |
| Reads supporting (aligned to +/- strand): ref base T (6/8); new base G (7/6); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.06e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01 | |||||||||||
CGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATA > NC_000913/433433‑433860 | cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGGG‑ACCGCCTa > 2:440933/1‑228 (MQ=255) gcgtGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGAGGGCGCAAAAGCTGGCCCA‑ACGAGGACGTTTTACCACGCATCCCACCCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGAAAGTT‑TACCACCAACGTGGGGGTGAA‑CCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTca < 2:512035/225‑1 (MQ=255) tGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTAAACGCGTTGCGGATGGCGGGTGGAAACGCCAAAGGTGCG‑CCGCCCCATGTCacac > 2:287682/1‑228 (MQ=255) ggggggggCGC‑AAAGGTGGCGGACAAGAGGTCGTTTTACCACGCTCCCCAACCC‑TAATGGTGGGGGCGTTATTGTCAAAGAGGGCAAATTTGGCGGTGAAGGT‑TACCACCACGGGGCGGGGGAA‑CCACATGCCGAAGTCCCCGCGTTGGGTAGGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcc < 1:418277/223‑1 (MQ=255) cgcgGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGACAAAGGGGCG‑ACCCCCCATGTTCCACCCGCACCCCTTCGCcatccaggggggaaccacccg > 1:648396/1‑210 (MQ=255) cgGGCGCTAAAGCTGGCGCA‑ACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAGCAAAGAGGCG‑ACCGCCCTAGGCACCCCCCAACCCGTGAGCCACAATGGTGCTACGCCACCGTGc > 2:89013/1‑228 (MQ=255) acccacgccaccAACCCGGAATGGCGGGTGCGCAGTTGAAAAAGATGGGAAAATGGTGGGT‑AAGGTGTACCACAAACGGGGGGGAAAC‑CCCAAGGCCGAAGCACGCGGGTTGCGATGGGCGGGGAAAAAAGACAAGGGGGGG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcg < 1:378841/218‑1 (MQ=255) ccACCCC‑GAATGCCGGGTGCTTCATTGTAAAAGGTGGCGAAATTGGCGGTGAAGGT‑TCCCACCAACGTGGGGGTGAA‑CCACAGGCCGAAGTACCCGCGTTGCTTATGGCGGGTAAAAAGGCAAAGGGTGCG‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑t < 2:19008/228‑1 (MQ=255) acatccc‑GGATGTCGGGGGTGTCATTGTCAAAGAGGGCGAAGTGTTCGGGGAGGGT‑CACC‑CCAATGGGGGGGGCAACCCACAGGCGGATATACGCGCGTGGCTTATGGGGGGTAAAAAAGCCAAAGGTGCG‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑t < 1:129885/224‑1 (MQ=255) ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGCGCTGTGACGCACCCCACGCCCCTGGCGTAGCGCGCGCGGTTGCC‑TCGATGCAAGATcc > 2:270211/1‑226 (MQ=255) cAAAGTGGGCGAAATTGGAGGGGAGGGT‑TCCCACCAACGGGGGGGTGAA‑CCACATGCCGAAGTCCCGGCGTTGCGTATGGCGGTGGAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGc < 2:195505/228‑1 (MQ=255) cGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCG‑ACCGCCCATGTCACACTCGAACCCTGTAGCCACCATGGTCGTACGCCACCCTGCTGTGACGCACTCATCGCCCCTTGCGTAGCGCGCGTGGTTGCC‑TCGATGCACGATCCTAACCCGCAGTGCCCTGGGC‑ggg > 2:137617/1‑225 (MQ=255) cGAAATTGTCGGTGAAGGT‑TACCACCAACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCG‑ACCGCCCATGTCACACTCGAACCCCGTAGCCATCATTGTGGTACGCCCCCCTCCGGTTACGCACTCCTCCCCCCCGGCGCAGGCCCCGGGGTGGCC‑CCGCTGGACAATACCAACCCCCAggcgcccgggcgggg > 1:195505/1‑215 (MQ=255) cgggTGACGGT‑TCCCACCAACGTGAGGGTGAA‑CCCCATGCCGAAGTACCCGCGTTGCGTATGGGGGTTAAAAAAGCCAAAGGTGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTa < 2:416687/224‑1 (MQ=255) caccaACGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCG‑ACCGCCTATGTCCCCCTCGAACCCTGTAACCATCATGGTCGGACGCCACCCTGCTGGGACCCACCCCTCGCCGCCTGCGTGACCCGCGTGCTGGCC‑CCGCTGGCAAATCCCACCCCCCAGCGGCCCGGGC‑GCGGCGTTTACCGTCCCCCACAGGCt > 1:58482/1‑228 (MQ=255) aCGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCGCCCATGTCCCACTCGCACCCCGTAGCCATCATGGTCCGACCCCACCGTGCTTTGCAGCACTCAACGCCCCCGGGGTCACCCCCGTCGTGGCC‑TCGCTGCAACGTCCCAACACCCAAGAGCCTGGGC‑CGGGGCTTTACCTTCtcgcacaaggcctga > 1:436590/1‑213 (MQ=255) aCGTGCGGGTGAA‑CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCA‑ACCGCCTATGTCACAATCGAACCCCGTAACCATCAAGGTCGGACGCCACCGTGCCGTTACGCCCTCATCACCGCTGGCGGAGGGCCCCGGGTTGCCTTCGTGGCAAGTACCAAACCCCCAGG‑GGCGGGGCGGGGGACTTACCCGTCTCCAAAAGGCGGGc > 1:491339/1‑227 (MQ=255) cccACAGGCCGAAGTCCCGGCTGGGCGTGGGGCGGTTGAAAAACCGAAGGGGGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCGTGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa < 1:554376/227‑1 (MQ=255) tGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGGC‑GCCGCCTATGTCACCCTCGAACCCCGTAGCCCACATGGTAGTACGCCCCCGTGCTCTGACGCCCTCATCCCCGCACGCGCATCGCGCGTCGGTGCC‑TCGGCGCAAACTCCTAACCCCCCGGGCGCCGGGC‑GGGGGCTTTTCCGCCTGCAACAGGCTCGCCTTTGCGCAAGCCACGGccg > 1:379341/1‑227 (MQ=255) acgttcaccgggtGCGTGGGGGGGGTAAAAAGGCCAAGTGTGCG‑TCCGCCTATGTCACACTCGAACCCTGTAGGCATCATGGTCGTACGCCACCGTGCTGTGCCGCACTAATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg < 2:94699/215‑1 (MQ=255) aaGTACACGTGTTGGGTTGGGCGGGTAAAAAAGCAAAAGGGGCG‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg < 2:436590/227‑1 (MQ=255) ttGCGTATGGCGGGTGAAAAAGCCAAAGGTGCG‑ACCCCCTATGTCACACTCGAACCCTGTAGCCATCAAGGTCGTACGCCACCGTGCTGTGCCGCACTCCTCCCCCCTGGCGTAGCACGCGTGGGTTCC‑CCGCTGCAAGCTCCCAACCCCCCGGTGGCCGGGC‑GTGGCGTTTTCCGTCTGCCACAGGCTGGCATCGCCGGCAGCCACgggcgcagggggagggaaaca > 1:10932/1‑209 (MQ=255) ggggaagggggttaaaaagGCCAAAGTGGCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg < 2:326606/213‑1 (MQ=255) gttGGGGGGTGAAAAAGCCCAAGGGTCGCACCGCCTATGTCACACGCGAACCCTGTAGCCAGCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg < 2:210069/225‑1 (MQ=255) aTGGGGGGTGAAAAGCCCAAGGGGGCA‑GCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa < 2:58482/228‑1 (MQ=255) tGGCGGGTGAAAATGCAAAGTGTCCG‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAACGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCC‑TCGATGCAAGATCCTAACCCGCAGGTCGCTGGGC‑GTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt < 2:379341/228‑1 (MQ=255) gTGAAAAAGCCAAAGGTGCG‑ACCCCCTATGTCACACTCGAACCCTGTAGCCATCAAGGTCGTACGCCACCGTGCTGTGACGCACCCATCGCCGCTGGCGGCTCGCCCGGGGTTGCC‑TCGATGCCAGATCCTAACCCGCAGGGCGCCGGGC‑GTGGTCTTTACCGCCCGCACCAGGCCGGCGTTGACGTCAGGCACGGCCCGACGATGAGGGAAACCGAGGAAATGaaaa > 1: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 < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |