Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,098,155 | 0 | A | C | 29.5% | 114.3 / 12.1 | 61 | L140F (TTA→TTC) | rhaS | DNA‑binding transcriptional activator RhaS |
Reads supporting (aligned to +/- strand): ref base A (29/14); new base C (0/18); total (29/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.86e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGAT > NC_000913/4098018‑4098293 | ctcGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTTGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTAcc > 1:856383/1‑140 (MQ=255) cTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGc < 2:163821/140‑1 (MQ=255) gAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCa < 1:279251/140‑1 (MQ=255) gAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGCCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCa < 2:735292/140‑1 (MQ=255) aaTCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAg < 1:207272/140‑1 (MQ=255) tCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTc < 2:985409/140‑1 (MQ=255) cAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcg < 1:307032/140‑1 (MQ=255) aGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcgc < 1:20877/140‑1 (MQ=255) gTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcgcg < 1:422353/140‑1 (MQ=255) gTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTcgcg < 1:489766/140‑1 (MQ=255) caAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGtt < 1:1012657/140‑1 (MQ=255) agagCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTa < 2:797463/140‑1 (MQ=255) agagCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTa < 1:1035588/140‑1 (MQ=255) gagCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTAt < 1:228026/140‑1 (MQ=255) gATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCATGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAAtt > 2:677581/1‑140 (MQ=255) ggCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTg > 1:380584/1‑140 (MQ=255) cAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGct > 2:949533/1‑140 (MQ=255) aTCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtg > 1:149380/1‑140 (MQ=255) cGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgctg < 1:100549/140‑1 (MQ=255) gTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTtgctgc > 1:844865/1‑140 (MQ=255) aCTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTaaa > 1:968635/1‑140 (MQ=255) tGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAg < 1:441613/140‑1 (MQ=255) gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa < 1:583692/140‑1 (MQ=255) gcgcGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCa < 1:780828/140‑1 (MQ=255) gTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGttt < 1:89692/140‑1 (MQ=255) ttAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTg < 2:576579/140‑1 (MQ=255) cacaGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGga > 2:291238/1‑140 (MQ=255) caGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga < 1:810800/140‑1 (MQ=255) gCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAAc > 1:62646/1‑140 (MQ=255) gTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCt < 2:149380/140‑1 (MQ=255) aTTGCAGCAGGTGCGACAGCTGGTTGCACACATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACGGCCAGTCGCGTGATCTTGTTTATGCAATTACTGCTCTCGCTGCGTACATGCATCTTGCAGGAAAGCCtcg > 2:868541/1‑138 (MQ=255) gcagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACa > 2:1032857/1‑140 (MQ=255) cagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAg < 1:150235/140‑1 (MQ=255) cagGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAg > 1:577824/1‑140 (MQ=255) ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAgcg > 1:51983/1‑140 (MQ=255) ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAgcg > 2:751842/1‑140 (MQ=255) ggTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAgcg > 2:564458/1‑140 (MQ=255) gCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCAt < 2:319604/140‑1 (MQ=255) gCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga < 1:824587/122‑1 (MQ=255) gCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGgaga > 2:824587/1‑122 (MQ=255) gCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCAt > 2:28029/1‑140 (MQ=255) gTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCaa < 2:968635/140‑1 (MQ=255) gTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCaa > 1:661112/1‑140 (MQ=255) gTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCaa > 1:308274/1‑140 (MQ=255) tGCACAGATGGCACAGCAGGAAGGGAAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt < 1:949533/140‑1 (MQ=255) tGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCACCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt < 2:551651/140‑1 (MQ=255) ggCACAGAGGAAACAGCAGGAAGGGAAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACt < 2:215387/139‑1 (MQ=255) gCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACtt > 2:658607/1‑140 (MQ=255) aTGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtct < 1:609911/128‑1 (MQ=255) aTGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGtct > 2:609911/1‑128 (MQ=255) tGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctg > 1:348285/1‑140 (MQ=255) tGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctg > 2:876859/1‑140 (MQ=255) gAACAGCAGGAAGGGGAAAAGGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggc < 2:618887/140‑1 (MQ=255) gAACAGCAGGAAGGGAAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggc < 2:856383/140‑1 (MQ=255) gcagGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTctggcctggc > 2:595502/1‑140 (MQ=255) cagGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggct > 1:391278/1‑140 (MQ=255) cagGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggct < 1:1032857/140‑1 (MQ=255) agGAAGGGGAAAATGATTTCCCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctg < 2:310105/140‑1 (MQ=255) agGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCctggctg > 1:72605/1‑140 (MQ=255) aTGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGCTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCAtttt > 1:460039/1‑140 (MQ=255) tGATTTACCCTCGACCGCCGGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTg > 2:650444/1‑140 (MQ=255) tttACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCg > 1:449102/1‑140 (MQ=255) tACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGAt > 2:240151/1‑140 (MQ=255) | CTCGCCGGGCTGAATCAGTTGCTGCCACAAGAGCTGGATGGGCAGTATCCGTCTCACTGGCGCGTTAACCACAGCGTATTGCAGCAGGTGCGACAGCTGGTTGCACAGATGGAACAGCAGGAAGGGGAAAATGATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTTGCCGAT > NC_000913/4098018‑4098293 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |