| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 433,649 | 0 | T | G | 46.9% | 3.8 / 10.1 | 33 | G65G (GGT→GGG) | ribD | fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase |
| Reads supporting (aligned to +/- strand): ref base T (2/15); new base G (12/3); total (15/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.12e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCA > NC_000913/433421‑433872 | aggCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGCGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGACGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GGCAAAAGCAAAGTGGC‑GACCGcc > 2:641543/2‑238 (MQ=255) gcgcgcCTGGGGACCTAGGCCGTGCAGGACGAGTATCACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGAGGTATTCCCACGCAACCCACCCAGAATGTCGGGTGCGTCATTGTCAAAGGAGGGGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCAAAGTCCGCGCGTTGCGAATGGCGGGT‑AAAAAACCCAAAGGGGC‑GACCGCCTATGTCACACTCg < 2:436321/238‑1 (MQ=255) gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGGGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACGCGGGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCg < 1:653957/238‑1 (MQ=255) gcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGaa < 1:625315/238‑1 (MQ=255) gcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGGGCAACGAGGCCGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAAcc < 1:534693/238‑1 (MQ=255) ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGG‑GGAAAAAGCAAAGGGGG‑GACCGCCTATGTCACAACCGAAcccc > 1:436321/1‑238 (MQ=255) ccctggCAGGACGAGTTTTCCAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCGGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg < 2:273406/234‑1 (MQ=255) tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACACGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAAcc < 1:831292/220‑1 (MQ=255) tGCAGGACGAGTATTACAGGGCGCGGGCGCGAAGGCTGGCGCACCGAGGGCGTTTTACCACGCATCCCAGCCCGAATGTCGGTTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGTTTCCCACCAGCTAGCGGGAGAACCACATGCCTAAGTTCCCGCGTTGGGTATGGCGGGT‑GAAAAAGCCAAAGGTGG‑Ga < 2:93004/198‑1 (MQ=255) gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGG‑GGAAAAACCAAAAGGGG‑GACCGCCCATGTCACACTCCGACACCGAAGCcat > 2:831292/1‑220 (MQ=255) ttAAATGGCGCGGGCGCTAAAGCGGGCGCAACGAGGACGTTTTACCACGCAGCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGTTGGCGAAATGGTCGGTGAAGGTTACCACCAAGGGGCGGGGGAACCACATGCCGAAGTACAGGCGTTGCGAATGGCGGGT‑AAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCg < 1:127374/203‑1 (MQ=255) cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGCGCTCTGACGc > 1:720585/1‑238 (MQ=255) aGGACGTTTTCCCACCCATCCCACCCCGAATGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcgc < 1:706309/238‑1 (MQ=255) cccAACCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCAAAATGGTCGGGGAAGGTTCCCACCACCGTGCGGGTGAACCACAGGCCAAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCAAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa < 1:578729/236‑1 (MQ=255) aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGG‑GGAAAAACCAAAAGGGG‑GACCCCCTAAGTCGCACTCGCACCCCGTAGCCAACAAGGGTGGACCCCCCCCGGCCGTGACGCACCCATCCCCGCCGGCGGTACGCGCGGGGTTGCCTCCATGCACGATccctcccc > 2:578729/1‑232 (MQ=255) gtggggTGCGTCTTGGTCAAACTTGGCAAAGTTGTGGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACTGGCGGAAGTCCACGCTTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTCC‑AACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc < 2:720585/235‑1 (MQ=255) aaTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑Ga < 1:570693/92‑1 (MQ=255) gAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGG‑GGAAAAACCAAAAGGGC‑GACCGCCCATGt > 2:570693/1‑92 (MQ=255) gTTACCACAAACGTGCGGGTAAACCACATGCCGAGTTACGCGCGTGGCGAAGGGCGGGT‑GAAAAAGCCTAAGGTGC‑GACCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcg < 1:650229/163‑1 (MQ=255) ttACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGG‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCCCAGGTCGCT‑GGGCGTGG‑CCCTTACCGTCTGCAACAGGCTGGCATTg > 1:567736/1‑239 (MQ=255) caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGG‑GAAAAAACCAAAAGGGC‑GACCGCCTATGTCTCACTCGGAACCTGGAACCATCATGGTTGGACCCCACCGGGCtgcgcc > 1:698152/1‑128 (MQ=255) cGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa < 1:192912/239‑1 (MQ=255) gCGGGTGACCCACATGCCGAAGTACCGGGGTTGCGATGGGGGGGT‑GAAAGAGCCAAAGGGTG‑GGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg < 1:379008/238‑1 (MQ=255) tGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAGGGTGC‑GGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt < 1:84638/146‑1 (MQ=255) aaGTACACGCGTTGCGTATGGCGGGG‑GAAAAAACCAAAGGGGC‑GCCCGCCTATGTTACACTCGCACCCCGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACCCACCCACCGCCGCTGGCGGAGCGCGGGTCGGTGCCTCTATGCAAGATCCTAACCCGCCGGGCGCT‑GGGCGTGG‑ACCTTACCCGCTGCGACAGGGTTGGCTTGACGTCCGTCCCGGGCCGCTGATGGGTggaaga > 2:784257/1‑234 (MQ=255) aGTACACGCGTTGCGTATGGCGGGTGGAAAAACCCAAAGG‑GC‑GACCCCCTTTTTCCCACCCGAACCCCGTAGCCCTCATGGGCGTACGCCCCCCTGCTGTGACGCCCTCCTTCCCGCCGGCGTAGCGCGCGTGGGTTCCTCCATGCAAGATCCCAACCCCCAGGTCCCC‑GGGGGTGGCATTTTACGGTC‑GCAACAGGCCGGCATTTACGTCCGCCCCCGCCCGCTGGTGGGGGAAGACg > 2:395809/1‑239 (MQ=255) ccgggggTTGGGTAGGGCGGG‑‑GAAAAACCCCAAGGGCCGGGCCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg < 2:567736/232‑1 (MQ=255) cacGCGTTGCGTATGGCGGGG‑GAAAAAGCCAAAGGGGG‑GACCGCCCATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTAGCGCGCGGGGCTGCCCCCATGCAAGa > 2:84638/1‑146 (MQ=255) tGCGTATGGCGGGT‑GAAAAAAGCAAAAGGGC‑GACCCCCTATGCCACCCTTGAACCCTGTAGCCATCATCGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTAGCGCGCGTGGGCTCCTCGAAGCCAGACCCTAACCCGCAGCACTCT‑GGGCGGGG‑GCTTTTCCGTCTGCAACCGGCCCGCCTTGACGGCAGCCCAGGCCGGCTGACGAGTGCAGCCGAACAAATggaa > 2:470499/1‑236 (MQ=255) gCGTATGGCGGGT‑GAAAAAGCCAAAGGGGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCaa > 1:416903/1‑181 (MQ=255) cGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg < 1:702267/226‑1 (MQ=255) gggAAAAAGCCAAAAGGGG‑GAACCCCTCTTTGCCACTCCAAACCTGTAGCGATCATGGGCGGACGCCCACGGGGTGTGACGCCCTCATCACCGCTGGAGTAGCGCGCGGGGGTGCCTCGCTGGCAGATCCTAACCCGCAGGTCGCC‑GGGCGTGG‑GCTTTTCCGTCTGCAACACGCTTGGAATGACGTCCACCACCGCCCGAGGGTGAGGGAAGGCGGGCAATTGaa > 2:702267/3‑226 (MQ=255) tggaaaaACCAAAGGGGC‑GACCGCCTATGTCACACTCGAAACCTGTAGCCATCATGGTCGGACCCCCCCCTGCTGTGACGCACCCATCCCCGCTGGCGGAGCGCGCGTGGGTGGCCCGGTGCAAGATCCCAACCCCCACGGCGCTCGGGCG‑GG‑ACTTTACCGGCTGCCAAAGGCTGGGATTGGCGTCAGCCCCGGCCCGGTGATGAGTGGAAGCGCGCACATTAATAAAAGGTTTTTCa > 2:617716/4‑239 (MQ=255) | TGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGT‑GAAAAAGCCAAAGGTGC‑GACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCT‑GGGCGTGG‑ACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCA > NC_000913/433421‑433872 |
| 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: ‑ |