| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 875,327 | G→A | *783* (TAG→TAA) | yliE → | putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 875,327 | 0 | G | A | 97.7% | 126.3 / ‑6.3 | 43 | *783* (TAG→TAA) | yliE | putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (18/24); total (19/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGT > NC_000913/875196‑875455 | gacggCGTAGAGAGCGTAGTCGTTGCGCGCGAATTAATCAAATTACGCTGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCTCCTGGGATAAAAGTGTAAAACTAGTAAAAGTGGGATTTACGt < 2:563484/136‑1 (MQ=255) tAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa < 2:230033/139‑1 (MQ=255) tAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa < 1:1197171/139‑1 (MQ=255) tAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa < 1:364798/139‑1 (MQ=255) tAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTCGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa < 1:20526/139‑1 (MQ=255) agagagCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa > 2:1039953/1‑138 (MQ=255) agagagCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa < 1:1039953/138‑1 (MQ=255) cAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc < 1:1335659/84‑1 (MQ=255) cAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc > 2:1335659/1‑84 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc < 1:1266054/82‑1 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTc > 2:1266054/1‑82 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa > 2:1233302/1‑84 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCa < 1:1233302/84‑1 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg < 1:1254667/99‑1 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg > 2:1254667/1‑99 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg > 2:602700/1‑99 (MQ=255) ggggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCg < 1:602700/99‑1 (MQ=255) gggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagca < 1:638028/139‑1 (MQ=255) aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcag > 1:1118690/1‑133 (MQ=255) aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcag < 2:1118690/133‑1 (MQ=255) aGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt > 1:589969/1‑95 (MQ=255) aGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTAcagt < 2:589969/95‑1 (MQ=255) aaaCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCt < 1:540642/131‑1 (MQ=255) aaaCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCt > 2:540642/1‑131 (MQ=255) aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGtt < 2:1178197/79‑1 (MQ=255) aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGtt > 1:1178197/1‑79 (MQ=255) aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa > 2:686860/1‑139 (MQ=255) aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGGAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa > 2:830594/1‑139 (MQ=255) ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCaa < 1:28218/71‑1 (MQ=255) ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCaa > 2:28218/1‑71 (MQ=255) ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTg < 1:226857/124‑1 (MQ=255) ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTg > 2:226857/1‑124 (MQ=255) ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGt < 2:595494/51‑1 (MQ=255) ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGt > 1:595494/1‑51 (MQ=255) ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTACAGCAGTTGCCTGCGGGATCTACACTCAAATGGt > 1:1134244/1‑139 (MQ=255) ccTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAgg > 1:1297708/1‑139 (MQ=255) tGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa < 2:1297708/139‑1 (MQ=255) tGGGATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTACAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa < 2:1134244/139‑1 (MQ=255) ggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACt < 2:957684/64‑1 (MQ=255) ggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACt > 1:957684/1‑64 (MQ=255) ggATAAAAGTGGAAAATTAGTAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACg > 1:1180569/1‑139 (MQ=255) tAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGt < 2:1180569/139‑1 (MQ=255) tAAAAGAGTAATTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGt < 1:686860/139‑1 (MQ=255) | GAAGGCGTAGAGAGCGAAGACGTTGCGCGCGAATTAATCAAATTAGGATGTGTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGT > NC_000913/875196‑875455 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |