Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 33.3% | 62.8 / 15.1 | 39 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (24/2); new base A (0/13); total (24/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACG > NC_000913/1824922‑1825193 | tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGagaaaag > 1:1068161/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATTGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGa > 1:657965/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGa > 1:99309/1‑140 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTg > 1:1206294/1‑140 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTc > 1:1282751/1‑140 (MQ=255) cccccGGCCCTTTTTTTTTTTTGTTATCGTTCAGGTAACGGTGGTTTTTTATCGGCAACTGGGGTGATTTACGTTGGTGGTGATCAGTTTTTTTTTA‑AACTTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCAc < 2:657965/138‑1 (MQ=255) cccccGGCCCGTTTTTTTCATTGTTATCGTTCAGGTACCGGTGGTTGTTTACCGGCAACTGGGGTTTTTTAGGTTAGTGGTGATCAGTTTTTTTTTACAATTTTTT‑AAAAAAAAAAAAAAAAACCGCCGATCCTGTCCAc < 2:1206294/138‑1 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 2:443796/1‑140 (MQ=255) ggCCCTTTTTTTTCTTTTTTATCGTTCAGGTACCGGTGGATGTTTATCGGCACCTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTTTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:1282751/140‑1 (MQ=255) ggCCCTTTTTTTTCTTTGTTTTCGTTCAGGTACCGGTGGATGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTTTTG‑GAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:1053572/140‑1 (MQ=255) ggCCCGTTCTTTTCTTTGTTATCGTTCAGGTACCGGTGGATGTTTATCGGCAACTGGGGTGATTTACGTTAGTGGTGATCGGGTTTTTTTTACAATTGTTT‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:99309/140‑1 (MQ=255) gccccTTTTTTTTTTTTGTTTTTGTCGAGGTAACCGTGGATGTTTCTGGGAAACTGGGTTTATTTACGTTGGTGGTAATCAGGTTTTTTTTACAATTGTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:1271871/138‑1 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:274555/1‑140 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:1271871/1‑140 (MQ=255) ttttttCTTTTTTATCGTTCAGGTACCGGTGGTTGGTTATGGGCAACTGGGGTGTTTTATGTTGGGGGTGTTCAGGTTTTTTTTACAATTGTTT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa < 2:274555/138‑1 (MQ=255) ttttttCTTTGTTATGGTCCAGGTACCCGTGGGTGTTTAGCGGCACCTGGGTTGTTTTGGGTTAGTGGGGATCCGGCTTTTTTAACAGGTGTTT‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa < 2:1406792/138‑1 (MQ=255) cAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg > 1:537075/1‑140 (MQ=255) ggTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg > 1:1156480/1‑140 (MQ=255) cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCttt > 2:843690/1‑139 (MQ=255) aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt > 1:1480455/1‑140 (MQ=255) aTGGTTATCGGAAACTGGGGTGTTTTGCTTTGGTGGTGATCGGTTTTTTTTTACAACTGTTT‑AAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt < 2:1480455/140‑1 (MQ=255) tttATCGGCAACTGGGGTGTTTTACGTTGGTGGTGATCAGGTTTTTTTTACAATTTTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt < 2:330850/139‑1 (MQ=255) gTTATCGGAACCTGGGGTGATTTCGTTTAGTGGTGATCAGTTTTTTTTTACAATTGTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt < 2:587071/140‑1 (MQ=255) ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTg > 1:376145/1‑140 (MQ=255) ttCGGCAACTGGGGTGATTTACGTTGGGGGTGATCAGGTTTTTTTTAAAACTGTTG‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc < 1:1290540/139‑1 (MQ=255) tCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt > 2:37875/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 1:63836/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 2:1416632/1‑140 (MQ=255) ggggttttttAGGTTAGGGGTGAGGGGTTTTTTTTTAAAATTGTTG‑GAGAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta < 2:883849/133‑1 (MQ=255) ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGGAAGATTTTTcagcag > 1:883849/1‑140 (MQ=255) ttAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa > 1:673922/1‑140 (MQ=255) tAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg < 1:37875/140‑1 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc > 2:1642274/1‑140 (MQ=255) tggtgATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCa > 2:1097409/1‑140 (MQ=255) gATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCAtttt > 1:1023264/1‑140 (MQ=255) gATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt > 1:905285/1‑140 (MQ=255) aTCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa > 2:161995/1‑140 (MQ=255) gCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:1254783/1‑140 (MQ=255) cAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGa > 1:179757/1‑140 (MQ=255) tcTGTTT‑GAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 2:1023264/139‑1 (MQ=255) ttGTTA‑GAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGt < 2:179757/139‑1 (MQ=255) tGTTG‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 2:673922/140‑1 (MQ=255) aaaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa < 2:1156480/135‑1 (MQ=255) aaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg < 2:376145/137‑1 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg > 2:1279872/1‑126 (MQ=255) gaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 1:1279872/126‑1 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACG > NC_000913/1824922‑1825193 |
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: ‑ |