Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,056 | 0 | G | A | 35.5% | 45.0 / 11.8 | 31 | intergenic (‑119/+84) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (18/2); new base A (0/11); total (18/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.21e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTC > NC_000913/1824807‑1825301 | gcgATGGATGCAAGGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGTTG‑ATCAGGCTTTTTTT‑ACAACTGTTTgag‑aagaga < 1:288974/251‑1 (MQ=255) gCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAgaaaagaaaa < 1:547352/250‑1 (MQ=255) cAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑cc < 1:498589/251‑1 (MQ=255) gCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATcc > 2:509864/1‑251 (MQ=255) gCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTc > 1:355599/1‑251 (MQ=255) gCCGGTGGGAACGTGCAAATTTAGCGCGTTTCACCGGCAGCATATCGCCACAGATTCACCGACATTTTACGCATATAAAAGTATTCTATACTCGCTCCCGGCCCGTACTTTTCACTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACCA‑GGGTGATTTACGTTAGAGGTG‑ATCAGGCTTTTTTT‑ACAATTGTTTGAG‑ATGTGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCACGGtc > 2:43491/1‑250 (MQ=255) gtgaaaaTTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACTTTTTACGCATATAAAATTATTCCATACCCGCCCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTACCCGTGGATGGTTATCGGCAACTG‑GGGTTTTTTACGTTAGTGGTG‑ATCGGGTTTTTTTT‑ACACCTGTTGGAA‑AAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc < 1:54260/247‑1 (MQ=255) aatttccgggctttacccgggCATCATTTCCCCACAGTTCCCCCGATTTTTTCCGCATATAAAAGTTTTCCATCCTCGCTCCGGGCCCGTTTTTTTTATTGTTATCGTTCAGGTACCGGTGGTTGGTTATGGGAAACTG‑GGGTTTTTTACTTTAGGGGTG‑ACCGGGTTTTTTTT‑ACAGCTGTGTGAG‑AAGAGAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg < 2:337406/232‑1 (MQ=255) aTTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcggg > 2:272247/1‑250 (MQ=255) gTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCAGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcggacctacgatt > 2:288832/1‑241 (MQ=255) aCCGGCAGCATTTAGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAACA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGtttttttg > 2:429484/1‑249 (MQ=255) ccGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACAGGCGGTCTTACATTTTTTGGc > 1:494377/1‑251 (MQ=255) ccGGCAGCATTCCCCCACAGTTCCACCGACTTTTGACGAATAACAAAGTTTTCCATCCCCGCTCCCGGCCTGTTTTTTTAATTTTTTTCTTTAGGGAAACGGTGGAGGGTTACGGGCACCGG‑GGGTGTTTTATGTTGGTGGTG‑ACCGGGTTTTTTTA‑ACAATTGTTTGAA‑AAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc < 2:87317/251‑1 (MQ=255) ccccacatattccccatcATTTTCCGCATACCAAAGTATTCCATCCTCCCCCCCGGCCCGTTCTTTTAATTGTTATCTTTCAGTTACCCGGGGATGTTTATCGGCA‑CTCCGGGTGTTTTACTTTAGTGGTG‑ATCAGGCTTTTTTT‑AAACCTTTTGGAA‑AAGAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattat < 1:361073/233‑1 (MQ=255) aCCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCTGGAAGATTATTCAGCAGTTGc > 1:209002/1‑249 (MQ=255) ttcacgataataaagttatccaatccTCGCTCCCGTCCGTTTTTTTTTATTGTTTTCGTTCAGTTACCGGTGGATGGTGCCCGACAACTG‑GGGTAATTTACGTTGGGGGTG‑ATCGGGTTTTTTTA‑ACACCT‑TTTGGG‑AAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt < 2:150764/225‑1 (MQ=255) tCCATTCTCGCTCCCGGCCTGTTTTTTTCTTTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTG‑GGGTGTTTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACACCTGTTTGAA‑AAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc < 1:478114/248‑1 (MQ=255) ccGCTCCCGGCCTGTTTTTTTCATTGTTATCGTTCAGGTACCCGTGGATGTTTATCGCCAACTG‑GGGTGATTTCCGTTAGTGGTG‑ATCAGTTTTTTTTT‑ACAACTGTTTAAG‑AAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAg < 1:195337/250‑1 (MQ=255) cGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGGTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc > 1:47351/1‑248 (MQ=255) tttttatttGTTATTGTTCAGGTACCCGGGGTGGGTTATGGGCAAGTG‑GGGTGATTTCCGTTAG‑GGTGGAT‑AGGTTTTTTTTTACACCTGTTGGAA‑AAAAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt < 2:549175/243‑1 (MQ=255) cTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTGCCGCTGGGAGGTCTGTCAGACGCTTCTCAAAAtt > 2:297925/1‑250 (MQ=255) gtacgagtaccggTGGATGTTTACCGCCAACTG‑GGGTTATTTACGTTAGTGGTA‑AACAGGTTTTTTTT‑AAATCTGTT‑GAGGAGGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCttt < 1:431117/239‑1 (MQ=255) gTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGc > 1:31400/1‑251 (MQ=255) gtttatggtcacctggggTGATTTACGTTCGGGGTG‑ATCAGGTTTTTTTA‑ACAACTGTTTGAG‑AAAAAAAAAGAAAC‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAgg < 1:509864/236‑1 (MQ=255) ttATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑AAAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTGCCCCTGGACGTTCTTTCATACGCTTCTCAAAATTAGCCTTATATTGCTTTTTCTGTTCCTGCGTCAGGAt > 1:399494/1‑251 (MQ=255) cAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGGTTTTTGGCTGAAAGATTATTTAGGAGTTGCAGGCATTTTAACTTTTGCCGCGGGAAGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGGtgttgtg > 2:178022/1‑249 (MQ=255) aaCTC‑GGGAGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGCAAGATTATTCAGCAGTTGCAGGCATTTTTCCTTTTGCCCCTGGACGTTTTGTCATACGCTTCTCAAAAATAGCATTTAATTTCTTTTTCTGTTTCGGCGTCAGGATGTTGTaa > 2:120395/1‑250 (MQ=255) gATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt > 2:202107/1‑110 (MQ=255) g‑ATCAGGTTTTTTTT‑ACAACTGTTTGAA‑AAGAAAAAAAAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCg < 1:295568/151‑1 (MQ=255) g‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGGTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCg > 2:295568/1‑151 (MQ=255) gCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTTTTCTTGGTTTCCATGTGCGCCAGCAt > 2:154813/1‑251 (MQ=255) t‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCtt > 2:99039/1‑251 (MQ=255) aaCTGTTTGAA‑AAGAAAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgc < 1:297925/251‑1 (MQ=255) ngag‑aagagaaAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTTTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTTTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttt > 2:250618/2‑250 (MQ=255) g‑aagagaaAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGAAAGGTAGTGGACAAGACCGGCGGTCCTTAGTTTTTTGGGTGAAAGATTATTCAGCAGCTGTAGCCATTTTAACTTTTGCCGCTGGTCCTTCTGTCAGACGCTTCTCAAAATTTGCATTACATTGCTTTTTCTGTTCCGGCGTCCGGATGTTGTACATTTTGTTCTGGGGTTCCATGTGCGCCCGCAAGTTAGATTTTGGCCGttcttc > 1:169921/1‑251 (MQ=255) | GCGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTC‑GGGTGATTTACGTTAGTGGTG‑ATCAGGCTTTTTTT‑ACAACTGTTTGAG‑AAGAGAAAAGAAAA‑CCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTC > NC_000913/1824807‑1825301 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |