Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,056 | 0 | G | A | 46.8% | 48.7 / 26.2 | 47 | intergenic (‑119/+84) | ves/spy | cold‑ and stress‑inducible protein/periplasmic ATP‑independent protein refolding chaperone, stress‑induced |
Reads supporting (aligned to +/- strand): ref base G (20/5); new base A (0/22); total (20/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.41e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT > NC_000913/1824906‑1825205 | aCCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑gagaagag < 1:871328/151‑1 (MQ=255) gACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑gagaagagaa > 1:685124/1‑150 (MQ=255) cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGagaaaag < 1:719722/151‑1 (MQ=255) cATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGagaaa > 2:578916/1‑149 (MQ=255) ttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGagaaaagaaa < 2:703330/150‑1 (MQ=255) aCGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAAcc < 1:94178/151‑1 (MQ=255) cGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAAccg < 1:840264/151‑1 (MQ=255) atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACCGGGGTGATTTACGTTAGGGGTGATCAGTCTTTTTTTACAACTGTTT‑GAGAAGAAAAAAAAAAAccgccg < 1:449645/149‑1 (MQ=255) tCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATc > 2:953827/1‑151 (MQ=255) aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTg > 2:828643/1‑151 (MQ=255) gTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTc > 2:366233/1‑151 (MQ=255) aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc > 2:469755/1‑150 (MQ=255) aTACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGc > 1:241498/1‑150 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt > 2:336204/1‑151 (MQ=255) cTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt > 1:737178/1‑150 (MQ=255) cGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt > 2:502235/1‑151 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg > 2:137643/1‑151 (MQ=255) ccGGCCCGTTTTTTTATTTTTTATCGTTCAGGTACCGGT‑GGATGTTTATCGCCAACTGGGGTGATTTACGTTAGGGGTGATCGGGTTTTTTTTACAACTTTTT‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg < 1:756919/151‑1 (MQ=255) ggCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt > 1:299265/1‑151 (MQ=255) cccTTTTTTTTCATTTTTATCTTCCAGGTACCCGT‑GGATGGTTATCGCCAACTCGGTTAATTTACTTTAGGGGTAACCTGTCTTTTTTTACACCTGTTT‑AAAAAAAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 2:90639/151‑1 (MQ=255) cccGTTTTTTTCTTTTTTATCTTTCAGGTAACCGG‑GGATGTTTATCGGCAACTGGGTTTTTTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTGTTT‑AAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 1:635570/150‑1 (MQ=255) cGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg > 2:542296/1‑151 (MQ=255) cGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg > 1:973241/1‑151 (MQ=255) tttCATTGTTATCGTTCAGGTAACCGG‑GGTTGGTTATCGGCAACTGGGGTGATTTACGTTAGGGGTGATCGGTTTTTTTTTACACTTTTTG‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc < 2:829494/148‑1 (MQ=255) cttGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg > 2:602366/2‑151 (MQ=255) ttttATCTTTCAGGT‑ACCGTGGGGTGGTTTGCGGCCACTGGGGGGTTTTCCTTTAGGGGGGACCACTTTTTTTTAACAACTATTG‑GAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:575437/148‑1 (MQ=255) ttgggttAACCGT‑GGATGGTAGTCGGCACGTGGGGTATTTAACGTGGGTGTTTATGTGGCTTTTTTTCCACTTGTGA‑GAGAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa < 1:431032/144‑1 (MQ=255) ttcggttAACGGT‑GGTTGTTTATGGGCAACTCGGGTGATTTACTTAAGGGGTGATCAGGTTTTTTTTACACCTGTTG‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt < 2:971720/143‑1 (MQ=255) gttagccGG‑GGATGGTTATCGGCACCTGGGGTTTTTTACGTTAGTGGTGATCAGGTTTTTTTAACAACTTTTT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt < 1:782270/139‑1 (MQ=255) ggTAACGGG‑GGATGGTTACCGGCAACTGGGTTAATTTACGTGGGTGGTGACCGGGTTTTTTTTACACTTTTTT‑GAGAAAAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGt < 2:395511/149‑1 (MQ=255) ggTAACCGG‑GGAGGTTTACCGGAAACTGGGTTGTTTTACGTTAGTGGTAATCGGTTTTTTTTAACACTTGTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt < 2:976047/150‑1 (MQ=255) agccGT‑GGATGGTTATGGACAACTCGGGTAATTTACGTTAGTGGTGACCAGGCTTTTTTCAAAACTGTTT‑GAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt < 1:846636/145‑1 (MQ=255) aacgtt‑GGATGGTTATCGAAAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTT‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt < 2:790259/145‑1 (MQ=255) accgggggATGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCGGGTTTTTTTAACAATTTTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt < 2:897400/144‑1 (MQ=255) ccGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt > 1:907109/1‑149 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg > 2:441395/1‑149 (MQ=255) ggaacctggggTGTTTTACGTTGGTGGTGACCAGGTTTTTTTTACAATTGTTG‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc < 1:787722/142‑1 (MQ=255) ggCAACTGGGGTGTTTTCCGTTAGTGGTGATCAGTTTTTTTTTAAAACTGTTT‑GAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt < 1:9436/149‑1 (MQ=255) aaCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc > 1:229548/1‑150 (MQ=255) gggTGTTTTGCTTTGGTGGTAATCAGGTTTTTTTTACAACTTTTT‑GAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc < 1:756912/151‑1 (MQ=255) taaGTGGTGATCAGGTTTTTTTTACAACTGTTT‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt < 1:168311/146‑1 (MQ=255) ggtgATCAGGTTTTTTTTACACCTGTTT‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc < 1:183998/146‑1 (MQ=255) gtttttttttAAACTTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg < 1:786022/147‑1 (MQ=255) ttttttaACAGTTGTTT‑GAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc < 1:83172/149‑1 (MQ=255) ttttttACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt > 2:500027/1‑151 (MQ=255) ttACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGGTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCt > 1:825899/1‑149 (MQ=255) acacttGTTT‑AAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTg < 2:737157/143‑1 (MQ=255) aCTGTTG‑GAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa < 1:828663/148‑1 (MQ=255) tGTTT‑AAGAAAAGAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc < 1:155173/145‑1 (MQ=255) gTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt > 1:220129/1‑151 (MQ=255) ttt‑GAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg < 2:374470/148‑1 (MQ=255) ctt‑GAGACGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCCTTAGTTTTTTGGCTGAAAGTTTTGTCAGCAATTGCACGCGCTTTTCCTTTTGCCACTGGCCGGTTTGTCAAAc > 2:43853/2‑147 (MQ=255) agaaaaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa < 1:428859/151‑1 (MQ=255) tttgagaaGAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa < 2:753406/147‑1 (MQ=255) aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa < 2:658561/149‑1 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt > 2:643507/1‑150 (MQ=255) | ACCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGT‑GGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTT‑GAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT > NC_000913/1824906‑1825205 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |