| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,466,824 | C→T | pseudogene (1433/2518 nt) | ydbA → | putative outer membrane protein, N‑terminal fragment |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,466,824 | 0 | C | T | 81.2% | 19.4 / 2.7 | 17 | pseudogene (1433/2518 nt) | ydbA | putative outer membrane protein, N‑terminal fragment |
| Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (7/6); total (9/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01 | |||||||||||
AAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCA > NC_000913/1466593‑1467064 | gtgctcctccgatcgtCACCAACGGTGGCACCGGCACGCAGAATAACGGTGATGACGCCACGGCAAACAACAACGGCAAACCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGTCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGTCGACAGCGCAACGGTGGATAACAAGGGCACCATGATCGTCACCGATCCGGAGTc < 1:241531/235‑1 (MQ=17) ccATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATcc > 1:89748/1‑249 (MQ=255) gCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACCACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGATCGTCACCGATCCGGAGTCCATCGGTATCCAGAt > 1:129841/1‑238 (MQ=255) naGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGATCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc > 1:94932/2‑242 (MQ=255) cGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGGACGGAAATTGCTGGCAATAACGGGAAGGTGATTTAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATGGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGATCGTCACCGATc < 1:293800/188‑1 (MQ=17) ccaccGTTGACGGAAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGATCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt > 1:88636/1‑251 (MQ=255) caccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGTTCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAACGGcggcaccgg > 1:319688/1‑246 (MQ=255) ggcGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGTTCGNCGANCGGGNTGNCAGCGACGTCGGGCANGGTATTGATATCTCCGGCGACAGCGTAACGGTGGATAACAAGGGNNCNNNGNTCGTCAcc < 2:227616/148‑1 (MQ=17) gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGATCGTCAcc > 1:227616/1‑150 (MQ=255) aCGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCCTGACCGTCACCGATCNGGANNNNNTcnnnnnnnnnn > 2:202688/1‑144 (MQ=17) aCGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGcggttgttaacaan < 1:202688/183‑14 (MQ=255) gATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGCGACAGAGGAAAGGTGGATAACAATGGCACCATGATCGTCACCGCTACGGAGTCGATCGGTATACAGCTTGACGGAGACCAGGACATCGTCAATAACGAAGGCGAGAGCACTATCCCCAATGGCGGCACCGGCAATACGATCAACGGAAAAGACGCCACCGCGACTAACAGTGGAAAAAACACTGTTGATGGAAAAAATTTCACggg > 1:222949/1‑251 (MQ=18) ggggcggTCACGGTATGGATACCACCGGCGAAAGCGCAACGGGGAAGAAAAAGGGCACCATGATCGACACCGATCCGGAGTCAATCGGTACCTAGATTGACGTCGCCCAGGCCATCGTCAATAACGAAGGCGAGAGAAGTATCACAAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCCCCGCGAATAACAGTGGAAAAACCCCAGTGGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGc < 1:144160/249‑1 (MQ=25) gtatggagatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCGTGAGCGTCGCCGATCCGGAGGCCATCGTTATCCAGATTGACGTGGACCAGGCCATCGTCAATAACGAAGGCgagan < 1:33410/116‑2 (MQ=17) gcncAACGNTGGANAACAAGGGCACCATGNTCGTCACCGAGCCGGAGTCNNTNGGTATNCAGAt < 2:129841/64‑1 (MQ=11) cgcAACGGTGGATAACAAGGGCACCATGATCGTCACCGATCCGGAGTCCATCGGTCTCCAGATCGACGGCGACCCGGCCCTCGACAATAACGAAGGCGAGAGCACTATCACCAATGGCGACACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTg > 1:52931/1‑179 (MQ=25) ggCACCATGATCGTCACCGATCCGGAGTCCATCGGTATCCAGATGGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa < 1:284101/251‑1 (MQ=255) | AAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCA > NC_000913/1466593‑1467064 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |