Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 29.2% | 87.5 / 14.3 | 49 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (32/2); new base A (0/14); total (32/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.49e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCT > NC_000913/1824922‑1825195 | tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaag > 1:26942/1‑140 (MQ=255) cAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaaa > 2:1716436/1‑140 (MQ=255) tttcccATCCTCCCCCCCGGCCCTTTTTTTTCATTGTTATCGTTCAGGTAACGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTGAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAAccgc < 2:1101545/136‑1 (MQ=255) ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg > 1:285763/1‑140 (MQ=255) tcccATCCTCGCCCCCGGCCCGTTTTTTTTATTGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCACCTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAAccgccg < 2:26942/138‑1 (MQ=255) tcaaTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTTTGTTAGCGGAGATCTGGCCTTTTTTACAACTGATTGAGAAGTGAAAAGAAAACTGCCGa > 2:1365726/4‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 2:1226652/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 2:1011686/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:710966/1‑140 (MQ=255) cATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATc > 2:1054790/1‑140 (MQ=255) aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc > 1:642949/1‑140 (MQ=255) tACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCt > 1:637585/1‑140 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCa > 1:1534131/1‑140 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:1541594/1‑140 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:1304172/1‑140 (MQ=255) cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg > 2:732983/1‑140 (MQ=255) ggCCTGTTCTTTTTATTTTTTTGGTGCAGGTAACGGTGGATGTTTATGGGCAAGTGGGGTAATTTACGTTGGTGGTGATCAGGTTTTTTTTACAATTGTTGGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:832837/140‑1 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:316252/1‑140 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGGGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAGAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 2:1569565/1‑140 (MQ=255) tCTTTTCATTTTTATCGTTCGGGAACCGGTGGATGTTTTTGGCCAATTGGGGTAATTTACGTTGGGGGTGATCAGTTTTTTTTTACACCTTTTTAAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCa < 2:1501031/140‑1 (MQ=255) tCTTTTCATTTTTATCGTTCAGGTAACGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa < 2:714198/140‑1 (MQ=255) tttttttATCGTTCAGGTAACGGTGGTTGTTTATCGGCACCTGGGGTTTTTTACGTTGGTGGTGACGGGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 2:386758/135‑1 (MQ=255) tGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg > 2:1236963/1‑140 (MQ=255) ttctttCAGGTACCGGTGGTGGTTTATCGGCACCTGGGGTGTTTTACTTTGGGGGTAATCGGGTTTTTTTTACATCTGTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa < 2:1760954/136‑1 (MQ=255) tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAACAAAACCGACGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg < 2:1452937/140‑1 (MQ=255) gaGGTAACGGGGGTTGTTTAGCGGCAATTGGGGTGTTTTACGTTGGGGGTGATCGGGTTTTTTTTAAATGTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:1642793/139‑1 (MQ=255) ggaaccggTGGGTGTTTAGCGGAAACGGGGGTGTTTTACGTTAGTGGTGATGGGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg < 2:595340/133‑1 (MQ=255) aCCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt > 1:221281/1‑140 (MQ=255) cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa > 2:225199/1‑140 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt > 1:116820/1‑140 (MQ=255) ggTTATCGGCAACACGGGTGATTTACGTTAGAGGAGATCAGGCTTTTTTTACAAATGTTTGAGAAGAGAAAAGAAAACCGCCGATCCAGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGatttt > 2:1326644/1‑140 (MQ=255) cctggggTGTTTTACGTTGGTGGTGATGGGTTTTTTTTAACACTTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa < 2:393332/136‑1 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGAAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGa > 1:217296/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 1:1463739/1‑140 (MQ=255) ggTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACACTTGTTTGAAAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt < 1:1156830/140‑1 (MQ=255) gTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc > 1:262086/1‑140 (MQ=255) gttttACTTTGGTGGTGATCAGTTTTTTTTTACAATTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:1541594/138‑1 (MQ=255) ttAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 1:1054790/140‑1 (MQ=255) taggtggtgATCGGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:803219/137‑1 (MQ=255) tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 1:1733747/1‑140 (MQ=255) tggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 2:35733/1‑140 (MQ=255) aTCGGGTTTTTTTAACAATTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa < 2:1698807/140‑1 (MQ=255) aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa > 1:62413/1‑140 (MQ=255) ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt > 1:1646004/1‑140 (MQ=255) ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt > 1:1727863/1‑140 (MQ=255) gCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGGTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:789740/1‑140 (MQ=255) tttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTg > 1:1287598/1‑140 (MQ=255) ttttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGc > 1:1698807/1‑140 (MQ=255) tttACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc < 2:1463739/140‑1 (MQ=255) aCAACTGTTGGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg < 1:1716436/140‑1 (MQ=255) gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCCGTTGCAGGCCTTTTACCTTTTGCCGCTGGACGttt > 1:258464/1‑139 (MQ=255) tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt > 2:220168/1‑140 (MQ=255) tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGt > 1:271363/1‑140 (MQ=255) gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACg > 2:824278/1‑140 (MQ=255) aaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg < 2:217296/137‑1 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCGGACGCt > 1:1208129/1‑140 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTTGAGGCATTTTACCTTTTTCCGCTGGACGTTCTTTCAGACGCt > 1:364693/1‑140 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCT > NC_000913/1824922‑1825195 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |