| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,947,850 | 0 | C | T | 51.5% | 25.6 / 50.0 | 33 | E40K (GAA→AAA) | hdfR | flhDC operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base C (4/12); new base T (17/0); total (21/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.13e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-06 | |||||||||||
TCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATG > NC_000913/3947707‑3947998 | tCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAGGGTCACACCCAGTTGATTTTcccctt > 2:615977/1‑145 (MQ=255) aTGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGa < 2:172938/151‑1 (MQ=255) tGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAAGTTCACACCCAGATGATTTTTCAGTTGTCTGAtt > 2:54695/1‑147 (MQ=255) cgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGCGGGGG‑GAAAAGGTTC‑CACCCAGTAGGTTTTTCAGTTGTTTGATTTGaaa > 1:621402/1‑150 (MQ=255) cgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAg < 1:721292/151‑1 (MQ=255) aCGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTCACACCCCGTTGATTTTTCAGTTGTCTGATTTGAAAGCt > 2:318113/1‑141 (MQ=255) aCGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GGAAAGGGTTACAACCCGTTGATTTTTCCGTTGTCTGATTTGAAACGGCActgccgg > 2:266770/1‑147 (MQ=255) ggcCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTc < 2:81682/151‑1 (MQ=255) tCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTTACACCCAGTTGATTTTCCCGGTGTCTGGTTTGGAAGCTCACTGCTGTCTGGGGCCGATAGAGCGAtt > 2:239733/1‑150 (MQ=255) cATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGAtt < 1:302668/149‑1 (MQ=255) gTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GGAAAGGTTTACACCCAGCTGGTTTTTCTGTTGTCTGGTTCGAAAACTTAATGCTGATTGGGGCCAATAGAGGGCTTCattccc > 2:66701/1‑144 (MQ=255) tttCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GGAAAAGGTCACACCCAGTTGGTTTTTCAGTTGTCTGATTTGAAAGCTCACTGCTGACTGGGGCAGATAGAGCGATTTAGCCGccc > 1:197686/1‑148 (MQ=255) ttCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTCACACCCAGTTGATTTTTCAGTTGTTTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTTAGCCGccc > 2:188739/1‑147 (MQ=255) tCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAGGGTCACACCCCGTTGATTTTTCAGTTGTTTGTTTTGAAAGCTCACTGCTGGCTGGGTCAGATAGAGCGATTCAGCCGccccac > 1:611925/1‑146 (MQ=255) ggCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAAGGTTACACCCCGTTGATTTTTTAGTTGTTTGATTCGAAAACTCACTCCTGAATGGGGCAGATAGAGCGATTCAGCCGCCCGACCcaa > 2:650327/1‑146 (MQ=255) agtTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAGGTTCACACCCAGTTGATTTTTCCGTTGTCTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTCAGCCGCCCGACCAAAGTGACGCGtt > 2:460026/1‑150 (MQ=255) cACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTCACACCCCGGTGATTTTTCCGTTGTCTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTCAACCGCCCGACCAAAGAGACGCGTTCGGCTaa > 2:692867/1‑150 (MQ=255) aCCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACt < 2:382990/151‑1 (MQ=255) cGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTc < 1:255289/151‑1 (MQ=255) gCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTCTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGa < 1:160321/151‑1 (MQ=255) ttAAACGGATATTGTTTCTGTGGC‑GGGTCGAAAAGG‑TCACACCCAGGTGGTTTTTCTGCTTTGTGTTTTTAAAACCAAACGCTGTCCGGGTCTGCTAAAACGATTCTACCGCCCGACCAAAGAGATGCGGTCCGCTCACTTCCAGGaa > 2:719433/1‑148 (MQ=255) ttAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTCACACCCAGGTGATTTTTCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGGCAGGTAGAGGGATTCAGCCGCCCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGaa > 2:349545/1‑148 (MQ=255) ttAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAAGGTTACACCCAGTTGATTTTTCAGTTGTCTGATTCGAAAACTCACTGCTGACTGGGGCAGATAGAGCGATTCAGCCGCCCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGaa > 2:28418/1‑148 (MQ=255) tAAACGGATATTGTTTCTGTGGC‑GGGG‑GTAAAGGTTCACACCCACTTGATTTTTCAGTTGTCTGATTTGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGAGTTAGCCCCCCGACCAAAGAGACGCGTTCGGCTAACTTTCAGGaaa > 1:333323/1‑148 (MQ=255) tAAACGGATATTGTTTCTGTGGC‑GGGG‑GAAAAAGATCACACCCAGATGATTTTTCAGTTTTCTGATTTGAAAACTCACTTCTGACTGGGGCAGATAGAGGGATTCAGGCCCCCGACCCAAGAGAAGCGGTTGGGTAAATTCCAGGGAAggt > 1:664459/1‑149 (MQ=255) aaaCGGATATTGTTTCTGTGGC‑GGGT‑GAAAAAGGTCACACCCCGTTTATTTTTCAGTTTTCTGATTTGAAAACTCACTGCTGACTGGGGCAGATAAAGCGATTCAGCCGCTCGACCAAAGTGACGCGGTCGGCTAACTTCCACGaaa > 1:703811/1‑147 (MQ=255) tCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGt < 2:723802/150‑1 (MQ=255) aaGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTcc > 1:706800/1‑151 (MQ=255) aGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCg < 2:326380/151‑1 (MQ=255) aGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCg < 2:589685/151‑1 (MQ=255) tCACACCCAGTGGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCaa < 1:609542/150‑1 (MQ=255) acCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAatta < 1:433237/151‑1 (MQ=255) ttttCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGtt > 2:624135/1‑150 (MQ=255) tCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATg > 2:235865/1‑150 (MQ=255) | TCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTTAAACGGATATTGTTTCTGTGGC‑GGGT‑GAAAAGGTTCACACCCAGTTGATTTTCCAGTTGTCTGATTCGAAAGCTCACTGCTGACTGGGTCAGATAGAGCGATTCAGCCGCTCGACCAAAGTGACGCGTTCGGCTAACTTCCAGGAAAGTTTTTAACAATTCCGTATCCACAGTGCTCTCCGCAAAATTATTTGTCGTTATG > NC_000913/3947707‑3947998 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |