| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,158,393 | 0 | G | A | 29.4% | 56.4 / 10.2 | 34 | E248E (GAG→GAA) | yqhG | DUF3828 domain‑containing protein YqhG |
| Reads supporting (aligned to +/- strand): ref base G (21/3); new base A (0/10); total (21/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.18e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAAAAAGAACAGCCGGAAGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTAGCCGGAACAATAAATAAAAGGGAACACTATATGAAAACGATTTTCAC > NC_000913/3158144‑3158643 | gCAAAAAGAACAGCCGGAAGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGgag < 2:155967/250‑1 (MQ=255) aaaaaGAACAGCCGGACGCCCTTGTAGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCAGTGGACGTGGGTAGTTTACGACTCTTACAGCCAGGAGGTT‑AA‑TGCCTTCTATAAAACCACCTTCCAGACGGCCAATAATCCAGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCAGGAGTTGATGAAATTCGCGTCC‑TGGAGaa > 2:354917/1‑250 (MQ=255) aaaGAACAGCCGGAAGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTGCCGACTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAAGCCGAGCTGTTTTGGCGAAGAGGCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAAc < 1:298697/251‑1 (MQ=255) aaCAGCCGGACGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAAccgcc > 1:182005/1‑251 (MQ=255) cccGAACTTTTTAAACATATTGCTGTTTATAGCGGGCGGTGGAGGTGGGTGGTTCCCGCCTCTTACCGCCAGG‑GGTT‑AAATGCCTTAAAAAAACCCCCCTTCAAGACGCCAAATAACCCCGATGAAGATATGCAAATAGAACGGCAATTTATGTACGACAACCCGATCTGTTTGGGGGAAGGGTCGCTATTTTCACGTGTGGATGAAATTCGAGGCC‑TGGAAAAAACCGCCGATTCCGCCCGCATTCATGt < 1:180227/250‑1 (MQ=255) gggggccggggctgtgggtgcTCTCCGACTCATCCCCCCAGGCGGTA‑AA‑TGCCTTATATAAACCCACCTTCAAGACGACCAATAACCCCTATGAATATATCCAAATAAACCGGCAATTTATTTCCTACAATCCGATCTGTTTTGGCGAAGACTCGCTATTTCCACGCTTGGAGAAAATTCGACTCC‑TGGAAAAACCCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGaa < 2:455901/238‑1 (MQ=255) ggCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAag > 2:205723/1‑249 (MQ=255) gCCGTGGATGGGGTTGGTTGCCGACTCTCACCGCCAGGCGGTT‑AA‑GGCCTTCTAAAAAACCACTTTAAAGACGGCAAATAACCCCGATGAAGATATGAAAAAAAAACGGCATTTTATTTACGACAAGCAGATCTGTTTGGGCGAAGAGTCGCTTTTTTCCCGCTTTGATGAAATTCGAGTCA‑TGGAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGc < 2:124361/251‑1 (MQ=255) cGTGGGGGGTTCCCGACTCTTACC‑CCAGGGGGTTTAA‑TGCCTTCTATAAACCCCCCTTAAAGACGGCCAATAATCCCGATGAAGATATGAAAAAAGAACGGCAATTTATTTGCGACAACCCGTTCTGTTTGGGCGAAGAGTCGCTTTTTTCACGCGTTGATAAAATTCGATTCC‑GGGAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTg < 1:427472/251‑1 (MQ=255) gtggtTTCCGCCTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCCAGACGGCCAATAATCCCGATGCAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGAACAATGGCAACAAAGAAGAGCAAGAACTGGtttt > 2:185305/1‑251 (MQ=255) tggtTTCCGACTCTAACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCAGCATTCATGTTCGTTTTACGATGACCAATGGCAACAACGAAGAGCACGAACTGGtttt > 2:147719/1‑250 (MQ=255) ctTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGACCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGa > 1:504758/1‑250 (MQ=255) ggTT‑AA‑GGCCTTCAAAAAACCCACCTTCAAGCCGGCCAATAACCCGGATGAAAATATGCAAAAAAAACGGAAATTTATTTACGACAATCCGATTTGTTTGGGCGAAGAGTCGCTTTTTCCACGCGTGGATGAAATTGGAGTCC‑GGGAAAAAGCCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGaaa < 1:205723/251‑1 (MQ=255) taAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATcc > 1:69834/1‑251 (MQ=255) taAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGGGGGAAATCGCTGATTTTATcc > 1:14569/1‑251 (MQ=255) caccTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg > 2:350054/1‑250 (MQ=255) cTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAAAGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCAGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATAGCTGATTTTATCCGCCCGAAc > 1:329559/1‑250 (MQ=255) gCCAATAACCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCt > 1:437486/1‑251 (MQ=255) ccAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCt > 2:89007/1‑250 (MQ=255) aGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAgg > 2:260008/1‑251 (MQ=255) atacAAATAGACCGGCAATTTATTTACGACAACCCGATCTGTTGGGGCGAAGAGTCCCTTTTTTCACGCGTTGATGAATTTCGGGTCC‑GGGAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCaa < 1:182666/247‑1 (MQ=255) aaaaaTAGAACGGAAATTTATTTACGAAAATCCGATCTGTTTGGGCGAAGAGTCGTTTTTTTCCCGCGTTGATGAAATTCGGGTCC‑GGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCaaaaa < 1:169576/249‑1 (MQ=255) aTAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTgc > 1:328976/1‑251 (MQ=255) aTAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGAAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTgc > 1:61212/1‑251 (MQ=255) aaGAACGACAATTTATTTACGACAATCCGTTCTGTTTTGGCGAAGAGTCCCTTTTTCCACGGGTTGATAAAATTCG‑GTCCCTGGAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTgcc < 2:153000/250‑1 (MQ=255) caacccGATCTGTTTGGGAGAACAGCCTCTTTTTTCACGCGTTCAGAAAATTCGAGCCC‑TGAAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGACTAAAGCAATGAGCTg < 1:147719/247‑1 (MQ=255) tGTTTTGGCGATGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACCGCGGCAGCCTACTTAAGCAGATTGAGGAAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATa > 1:489309/1‑250 (MQ=255) gagcgggcaaaGAGTCGCTTTTTTCGCGCTTTGAAAAATTTCGAGCCC‑GGGAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATGGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAAc < 1:273246/242‑1 (MQ=255) ttttGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAAc > 2:160743/1‑250 (MQ=255) ggagAAGAGTCGCTATTCTCACGCTTTGAGAAAATTCGAGTCC‑TGAAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTa < 1:260008/248‑1 (MQ=255) ggCGACGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGCTTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAAc > 2:273246/1‑246 (MQ=255) aaGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGTAAACAGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCACCAACGACGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCACTGAGCTGAATTAAAt > 1:552585/1‑239 (MQ=255) gtatgaAATTCGAGTCT‑TGAAAAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTGTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTAGCCGGAACAATAAATAAAAGGGAac < 2:148324/248‑1 (MQ=255) gTCC‑TGGAGCAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTAGCCGGAACAATAAATAAAAGGGAACACTATATGAAAACg > 1:456247/1‑251 (MQ=255) cc‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTAGCCGGAACAATAAATAAAAGGGAACACTATATGAAAACGAt > 1:412930/1‑251 (MQ=255) gaAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTAGCCGGAACAATAAATAAAAGGGAACACTATATGAAAACGATTTTCAc > 1:571594/1‑251 (MQ=255) | GCAAAAAGAACAGCCGGAAGCCTTTGTTGCCGAACTCTTTGAACATATTGCTGATTATAGCTGGCCGTGGACGTGGGTGGTTTCCGACTCTTACCGCCAGGCGGTT‑AA‑TGCCTTCTATAAAACCACCTTCAAGACGGCCAATAATCCCGATGAAGATATGCAAATAGAACGGCAATTTATTTACGACAATCCGATCTGTTTTGGCGAAGAGTCGCTATTTTCACGCGTTGATGAAATTCGAGTCC‑TGGAGAAAACCGCCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGCCTACTTAAGCAGATTGAGGCAAAAACTGCCGCCAGATTAAAGCAATGAGCTGAATTAAATAACAATTAGCCGGAACAATAAATAAAAGGGAACACTATATGAAAACGATTTTCAC > NC_000913/3158144‑3158643 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |