Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 24.5% | 111.2 / 15.7 | 53 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (38/2); new base A (0/13); total (38/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.25e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824930‑1825196 | tATTCCATACTCCCCCCCGGCCCGTTTTTTTCATTTTTATCGTTCAGGTAACGGTGGATGGTTATCGGCAACGCGGGTGTTTTAGGTTAGTGGTGATCAGGTTTTTTTA‑ACAACTGTTTGAAAAGAAAAAAAAAAAccgc < 2:882162/140‑1 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:1516643/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:98000/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:1934441/1‑140 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg > 1:1404973/1‑140 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg > 2:416542/1‑140 (MQ=255) cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGt > 2:1990278/1‑140 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTc > 1:1050089/1‑140 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:1037742/1‑140 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:1951294/1‑140 (MQ=255) cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGGGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg > 2:1035099/1‑140 (MQ=255) tcttttctttGTTATGTTTAAGTTACCGGTGGATGTTTAGCGGCACGTGGGGTGTTTTACTTTGGTGGTGATCAGGTTTTTTTT‑ACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa < 2:1153382/140‑1 (MQ=255) cTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCaa > 1:500241/1‑140 (MQ=255) ttttCTTTGTTTTCGTTAGGGAAGCGGTGGTTGGTTACCGGCACCTGGGGTGTTTTCTGTGAGGGTTGACCGGTTTTTTTTT‑ACAATTTTTTGAAAAAAAAAAAGAAATCCGCCGATCCTGTCCACCGCATTACTGCAAg < 2:500241/140‑1 (MQ=255) tCATTTTTCTCGTTCAGGAAACGGTGGTTGGTTTTCGACAACGCGGGTGTTTTACGTTAGTGGGGAGCGGGTTTTTTTT‑ACAAGTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 2:1044266/140‑1 (MQ=255) tCATTTTTATCGTTCAGGTAACGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGACCAGGTTTTTTTT‑ACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 2:1934441/140‑1 (MQ=255) tGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg > 1:2007961/1‑140 (MQ=255) ttATTGTTCAGGTAACGGGGGTTGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGGTTTTTTTT‑ACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc < 2:98000/140‑1 (MQ=255) cgggTACCGGTGGATGGTTATCGGCAACTGGGTTGATTTACGTTAGTGGTGATCGGTTTTTTTTT‑ACAATTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:2134990/138‑1 (MQ=255) ggTACCCGGGGATGGTTATGGGCAACTGGGGTGTTTTCGGTTAGTGGTGACCAGGTTTTTTTT‑AAAATTTTTAGGGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg < 2:1273873/140‑1 (MQ=255) ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg > 1:936939/1‑140 (MQ=255) aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt > 2:613324/1‑140 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt > 2:1038849/1‑140 (MQ=255) aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt > 1:1351423/1‑140 (MQ=255) tttATCGGCACCTGGGTTGTTTTAGTTTAGTGGTGA‑CCGGCTTTTTTTAACACTTGTTGAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt < 2:2207618/139‑1 (MQ=255) gTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt > 2:493216/1‑140 (MQ=255) tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg > 1:1681803/1‑140 (MQ=255) ggCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGa > 2:419124/1‑140 (MQ=255) ggCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGa > 1:2167155/1‑140 (MQ=255) cAACTCGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTT‑ACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa < 1:416542/140‑1 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTAAAATTTTTTGCTGAAAgg > 2:5571/1‑139 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 2:2002007/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 2:485674/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 1:224807/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTAAGTTTTTTTTCTGAAAAat > 2:1875385/1‑140 (MQ=255) cTCGGGCGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGGTTTTTGGTCGAAAGat > 1:154870/1‑140 (MQ=255) gttttACGTTGGTGGTGATCAGGTTTTTTTT‑ACAATTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 1:1990278/138‑1 (MQ=255) gttttACGTTAGTGGTGATCAGGTTTTTTTT‑ACAATTTTTTAAAAAAAAAAAAAAAAACCGCCGATTCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:629017/138‑1 (MQ=255) ttggtggtgATCGGTTTTTTTTT‑ACATTTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:1516643/137‑1 (MQ=255) gtggtgATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGCAAAATTATTTAACAGTTGCAGGc > 2:2225869/1‑140 (MQ=255) gtggtgATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc > 1:1859810/1‑140 (MQ=255) tggtgATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 1:1550520/1‑140 (MQ=255) aTCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa > 1:765586/1‑140 (MQ=255) cAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc > 2:493969/1‑140 (MQ=255) aGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCt > 1:1478151/1‑140 (MQ=255) ggCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt > 2:1944849/1‑140 (MQ=255) gCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 2:1896132/1‑140 (MQ=255) gCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:659045/1‑140 (MQ=255) ttttttt‑AAAATTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg < 2:936939/140‑1 (MQ=255) tttttt‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGc > 2:225874/1‑140 (MQ=255) ttacttGTTGGAGAAGAGAAAAAAAAGCCGCCGATCCTGTCCGCCGCATTGCTGCAAGGTAGTGGACAAGATCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg < 2:154870/135‑1 (MQ=255) aaaattGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg < 2:224807/135‑1 (MQ=255) tGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 1:360818/140‑1 (MQ=255) gaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg > 1:548242/1‑140 (MQ=255) aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa < 2:765586/140‑1 (MQ=255) aaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa < 2:1478151/135‑1 (MQ=255) gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg > 2:2120841/1‑140 (MQ=255) aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt > 2:1989886/1‑140 (MQ=255) | TATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824930‑1825196 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |