Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 433,649 | T→G | 29.9% | 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 | 29.9% | 22.4 / 3.1 | 32 | G65G (GGT→GGG) | ribD | fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase |
Reads supporting (aligned to +/- strand): ref base T (6/16); new base G (6/4); total (12/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01 |
TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCT > NC_000913/433419‑433870 | tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGAACg > 1:370558/1‑238 (MQ=255) aaGAGAACGCGCGCATGAAGGACTAAGCCGTGCGGGACGGGTATTACATGGCGGGGGCGCTAAAGCTGGCGCCAAGAGGACGTTTTACCACGAATCCCAACCCGAGTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAGGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGt < 2:70904/238‑1 (MQ=255) ggaggaCAAAGCCGTGCAGGACGAGTATGCCAGGGCGCGGGCTGTAAAGCGGGCGCAACGAGGCCGTTTTCCCACGCTACCCACCCCGAATGTCGGGGGCGTCATTGTGAAAGATGGCGAAATTGTCGGTGAAGTTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGt < 2:350628/238‑1 (MQ=255) cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGGGGACGCCCCAGTCACCCCCGCACCCTTTAGCcat > 2:287053/1‑238 (MQ=255) cggtgaGGACGAGTATTGCAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTGGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc < 2:166748/234‑1 (MQ=255) gaGGACGAGTTTTACATGGCGGGGGCGCTAAAGCGGGCACAACGAGGACGTTTTACCACGCACCCCAACCCGAATGCCGGGTGCGTCATTGTCAAAGATGCCAAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTa < 2:255401/237‑1 (MQ=255) gACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGACCGCCTATGGTACAATCGAACCCTGTAGCCATCATGGTCGTACGc > 1:337784/1‑238 (MQ=255) cgggtcttACAGGGCGCGGGCCGTAAAGCGGGCGAACAGAGGGCGTTTCACCCCGCATCCCACCCGGAATGGCGGGGGCGTCATTGTCAAAGATGGCAAAATGGTCGGTGAAGGTTTCCACAAACGTGGGGGTGACCCACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa < 2:391009/232‑1 (MQ=255) gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGTGCGACCGCCGATGTCACACTCGAACCCCGTAGCCATCATGGGCGTGAGCCACCgg > 1:211973/1‑238 (MQ=255) tATTACATGGGGCGGGCGTTAAAGCGGGCGCAACGAGGACGTTTTACCACGCAGCCCAACCCGAATGTCGGGGGCGCAATTGACAAAGAGGGCAAAATTGTGGGGGAAGGTTACCCCCAACGTGCGGGTGAACCACAGGCCGAGGTACACGGGTGGCGTATGGCGGGGGAAAAACCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg < 2:262926/239‑1 (MQ=255) aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGCGGAAAAAACCAAAGGGGCGGCCGCCTATGTTCCACTCGGACCCTGTAGGCACCAACGGCGCACGCCACCGTGCtgt > 1:68769/1‑239 (MQ=255) cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCCCACTCGAACCCTGTAGCCATCCTTGTCGTACGCCCCCGTGCTGTGACGc > 2:178563/1‑238 (MQ=255) ccTAAAGCTGGCGAAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTAACCCACATGCCGAAGTACACGCGTGGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATc < 2:383641/237‑1 (MQ=255) tACCAAGCATCCTACCCCGAATGTCGGGTGCGTCATTGTCAAGGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt < 2:333928/239‑1 (MQ=255) ccACGCACCCCAACCCGAATGGCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg < 2:245655/239‑1 (MQ=255) cATCCCAACCCGAAGGTGGGGGGCGGCATTGTCAAAGAGGGCGAAATTGGCGTGGAAGGTTCCCCCCAGCGTGCGGGTGAACCCAATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa < 2:337784/239‑1 (MQ=255) aacccagcgggAATGTGGGGTGCGTAATGGTCAAAGATGGCGAAGTTGTGGGTTAAGTTTACCACAAAGGTGCGGGTGAACCACATGCGGAAGTACACGCGTTGCGTATGGCGGGTGAAAAGGCCAATGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa < 2:68769/228‑1 (MQ=255) tCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCACAAGGGCGCCCCCCCCCGGCACCCTCTGACCCTCTACCCCTCATGGGCGTCCGCCCCCGCCCTGTGACCCCCTCCGCGCCCCCGCCGTGGAGCGCGCGGTTACCTcccctgacg > 2:254920/1‑231 (MQ=255) gtggggTGCGGAATTGTAAAAGGTGGCAAAATTGGCGGTGAGGTTTACCGCCACCGTGCGGGTGAACCACATGCGGAAGTCCCCGCGTGGCGTAGGGGGGGTGAAAAAGCCAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc < 2:132066/235‑1 (MQ=255) cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCCCCTATGTCACACCCGAAACCTGTAGCCCTCATGGGCGTACGCCACCGTGCTGTGACGCACTCATTCCCGCCCGCGTCACGCGCGTGGTTGCCTCGATGCCAGAACCTAACCCCCGGGTGCCTGGGCGTGGaacttaacct > 2:2056/1‑229 (MQ=255) aTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCAAATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCa < 1:287053/239‑1 (MQ=255) gTGAAGGTTACACCCAACGTGCGGGGGACCGACATGCGGAAGTACACGGGTTGCGTATGGCGGGGGAAAAAGCAAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc < 1:34361/238‑1 (MQ=255) aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCCAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCCTCCTGGTCGTACGCCACCGTGCTGTGACGCACTCCTCCCCGCCGGCGTAGCGCGCGGGGTTGCCCCGATGCAAGATCCCAAACCGCCGGGCGCTGGGCGGGGGATTTCCCGTCCGCACCCGGCCGGCCCTGACCTCCGc > 2:292427/1‑238 (MQ=255) ggggggcgaaaCCCCAAGCCGAAGTACGCGCGTGGCGTATGGCGGGGAAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa < 1:178563/229‑1 (MQ=255) tGCGGGGGAACCACAGGCGGAAGTACACGCGTTGCGTATGGCGGGTGAAAAGGCCAAAGTGGCACCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa < 1:16329/237‑1 (MQ=255) gggggTGAAACACATGCCGAAGTACACGCGGTGCGTAGGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTCCGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg < 2:471326/236‑1 (MQ=255) aGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGCCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCCCGATGCAAGAGCCTAACCCCCGCGTCCCTGGGCGTGGCCTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCCGATGGTGAGTGAAGcca > 2:50828/1‑238 (MQ=255) acggggtGCGTAGGGGGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa < 2:370558/233‑1 (MQ=255) ggctttGCGTATGGCGGGTAAAAATCCAAAAGGTTCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt < 2:165296/235‑1 (MQ=255) gggttaggggggggaaaaaaGCCAAAGGGGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAt < 2:342379/224‑1 (MQ=255) gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGGACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGGGGACCTTACCCGCTGCAAAAGGCTGGCATTGACCTCAGCCCCGGCCCGAGTGTGAGTGAAGCCGAGCAATTTAATAAAAGCGTtct > 2:322576/1‑238 (MQ=255) gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtct > 2:460751/1‑238 (MQ=255) | TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCT > NC_000913/433419‑433870 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |