| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,384,808 | C→T | K305K (AAG→AAA) | yfbK ← | IPR002035/DUF3520 domain‑containing protein YfbK |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,384,808 | 0 | C | T | 85.7% | 30.2 / ‑5.5 | 14 | K305K (AAG→AAA) | yfbK | IPR002035/DUF3520 domain‑containing protein YfbK |
| Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (8/4); total (8/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACCAGA > NC_000913/2384578‑2385053 | gcTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCACCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGCTTTTGGATCGTCAATGACAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGc > 1:218407/1‑251 (MQ=255) gagagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCAGTTc > 1:43188/1‑251 (MQ=255) aTCAATGGAGCTGTAGTTGCCGTTACCAACACCGGCAATTCGCACCAACATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTTACCTTTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCgg < 2:51450/251‑1 (MQ=255) aCCAACATCGGCAATGCGCACCATCATTGCCGCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGTTCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGAGGCCGCCCTAAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCATTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCgg < 1:256230/251‑1 (MQ=255) cccGCTGTTTTTTGACCATTGATTCAATCGCTTTTGGATCGTCACTGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATACGCCAGTTCCA‑GCCCGACACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGCTAGAAGGCAATGCAGTACGGGAGTCGCAAGCg > 2:280047/1‑251 (MQ=255) ctgtttgttgcacttttatccAATCGATTTTGGTTCTTCAAACCCAACTTTAAAGTCACC‑TCCGTGGGGTAATAAATTGCGTTGGATTCCGCCCTTA‑TTAACCCCTTCCGCGGCCTGCGGATAAGCCTGTTCCT‑GCCCGTCCCCGCCATTGGTTCTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGTAGATAGAAGGCAATGCAATACGGGAGTCTCCAGCGTa < 2:123497/229‑1 (MQ=255) ttttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGTGTCGCCAGCGTAGGTCACg > 1:49358/1‑251 (MQ=255) tttGACGATTGTTTAAATCGATTTTGGATTGTCAATGCCAAAGTTAAAGTCACAGTCAGTGGT‑TAATAAAATGAGATTGATTCCGCACTTAATAAACACTTTG‑GTCGCATGCTGATAAGACAGTTCCT‑GCGAGGAAACTCAATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGAGGCTTTAATTTCCGCCTTATGACTACAGGGGATAGAAGGAAATGATATACggggaacgtctagggaggtcacg > 2:153268/1‑230 (MQ=255) aCGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCaata > 1:118481/1‑250 (MQ=255) gTCAGTGGC‑TAAAAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACt < 1:2708/251‑1 (MQ=255) cccaacaggccccctttccattcttGATGCCGCCTTTAATAC‑CCCCTTCCGCCGCCTGCTGATAACCCAGT‑CCACGCCCGGCGCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGTCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGaa < 2:252364/225‑1 (MQ=255) aaTGCGATTGCTGCCTCCCTTAATAACCCATTTC‑GGCGCCTGTTGATAATCCATTTCCA‑GCCGGGCTCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGTCCAATAATTTCAACGAAGACTGGATAAGTGGCAAAc < 2:263478/251‑1 (MQ=255) tGATGCCGCCCTTAATAAACCCTTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGaa > 1:252364/1‑224 (MQ=255) aCCCTTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACCAGa > 1:139325/1‑250 (MQ=255) | GCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGC‑TAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTC‑GTCGCCTGCTGATAAGCCAGTTCCA‑GCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACCAGA > NC_000913/2384578‑2385053 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |