| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 523,335 | 0 | T | G | 58.8% | 2.7 / 12.2 | 17 | G25G (GGT→GGG) | rhsD | protein RhsD |
| Reads supporting (aligned to +/- strand): ref base T (2/5); new base G (10/0); total (12/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.39e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAAAATGTTAGTTTAATAAAATATTGAAAGTGACCTGTAATAACAGTTGTTGTTGATTGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCC‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATG > NC_000913/523086‑523585 | aaaaaTGTTAGTTTAATAAAATATTGAAAGTGACCTGTAATAACAGTTGTTGTTGATTGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGgtg < 1:60565/251‑1 (MQ=255) aataaAATATTGAAAGTGACCTGTAATAACAGTTGTTGTTGATTGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATAAGCGGACAACCAGAGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTATCCTGGGGTCGGCAGGGGTAAGAATTggtggg > 1:270185/1‑247 (MQ=255) tgaGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCCGCGGGGAGTTGTCGGTGTGTCcgggccggggggctt > 1:105500/1‑241 (MQ=255) aaGTTGATGTGAAAAATATATAAATACATTAGCTGGTCGTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACGTTATATGAGCGGAAAGCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACGTCGGGC‑ATCCCg < 2:109071/251‑1 (MQ=255) cTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGGGGCGGGCTCCGGGTGTTCGGGGGGGGTGACCTCGGGG‑AACCCCGGAAATCCCCTGCCggggggggagggggtgcc > 1:25555/1‑238 (MQ=255) attttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGTTCGGCCGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCGCGGGGTGTCCGGGGGGGGTGACTTCGGGG‑AACCCCGTAAATCCCCCTCTGGGGGCGGAGGTGGTGCCCGGCggggagggccc > 1:206529/1‑241 (MQ=255) ttgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGGTGACCTCGGGC‑AACCCGGTAAATCCGCTGCCGGGGGCGAAGGGGCCGCCCGGCGAGGCGGACCTTTCGCTgcccggg > 1:79824/1‑250 (MQ=255) tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGACCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGTGTGTCCGGGAGGGGTGACCTCGGGC‑AACCCCGGAAATACGCTGCTGGGGGCGAAGGGGCTGCCCGGCGAGAAGGGACTTGCGCTGACCGGCCCGccgc > 1:34517/1‑248 (MQ=255) cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAAATGGCGCGCCCACCGGCGTGGCGTGCTCCGGGTGTCCGGGCGGGGTGACCTCGGGC‑CACCCCGGAAATCCCCTGCTGGGGGCGAAGGGGCGGCCCGGCGGGGCGGGCCCTGGGCCGCCCCGCccgccgcc‑gg > 1:218921/1‑250 (MQ=255) aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGTGAGAATTGGCGCGCCCCCCCGCGTGGCGTGCTCGGGGGGTCCGGGGGGGGTGACTTTGGGCCACCCCGGAAAAACC‑CCGCCGGGGGGGGAGGGGCTGCCCCGCGAGACGGGCCTTGCGCTGCCCCGCCCCCTGCC‑gctcctttcctcccccc > 2:168110/1‑235 (MQ=255) gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGACTTCGGGC‑AACCCCGTAAATCCGCTGCTGGGGGCGGAGGTGCTGCCCGGCGGGACGGACCCTGCGCTGCCCGGCCCGCCGCC‑G‑GTCCTTTTCTCCCGCACCTAcagccgc > 1:228883/1‑248 (MQ=255) aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggc < 1:160563/120‑1 (MQ=255) cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGTGTGGCCGGGCGGGGTGACTTCGGGC‑CACCCGGTAAAACCGCTGCTGGGGGGGGAGGTGCTGCCCGGCGCGACGGACCCTGCGCCGCCCGGCCCCCCGCCCG‑TTAATT‑TCTCCCCCCACTACATCAGCCACCGGACGAAGACGCCCGCAcccggg > 1:163861/1‑247 (MQ=255) cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCC‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTgg < 2:251799/249‑1 (MQ=255) aTGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCCGGCGGGGTGACTTTGGGC‑AACCCGGGAAATCCGCTGCTGGGGGGGGAGGTGCTGCCCGGCGAGACGGACCCTTCGCTGCCCCGCCCGCTGCC‑GGTTC‑TTCTCTCCCGCCCCCACCGCAGCTACCCGACGAAGACCCCCGCCCCGGGGGGCGGTTTTGGCCCCCGCTGGGAAgcg > 1:279356/1‑251 (MQ=255) cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCC‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGa < 1:158942/251‑1 (MQ=255) ggcgggTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTTCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCC‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTATCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTAc < 2:151663/247‑1 (MQ=255) ggCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCC‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTAc > 1:7101/1‑251 (MQ=255) tgtAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGAGTGCCCGGGCGGGATGACTTAGGGC‑AACACGGTAAATCCGCTGCTGGGGGCGAAGGAGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCA‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTACCGGAAGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGCAAGCGCCTTCTGATATCCGCTTACCGCTACGtgatg > 1:294398/1‑251 (MQ=255) | AAAAATGTTAGTTTAATAAAATATTGAAAGTGACCTGTAATAACAGTTGTTGTTGATTGAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGC‑AACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCC‑G‑TTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATG > NC_000913/523086‑523585 |
| 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: ‑ |