Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,056 | 0 | G | A | 44.4% | 67.0 / 26.5 | 54 | 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 (25/5); new base A (0/24); total (25/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.09e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT > NC_000913/1824907‑1825205 | ccGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTgagaagaga < 1:507446/151‑1 (MQ=255) gACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa < 2:914569/151‑1 (MQ=255) gACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTGTTTTACAACTGTTTgagaagagaa > 1:155900/1‑150 (MQ=255) cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc > 1:40797/1‑151 (MQ=255) cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc < 2:269361/151‑1 (MQ=255) tCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTG‑TTGAGAAGAGAAAAGAAAACCGCCGATcc > 1:482668/1‑151 (MQ=255) aaattttcccATCCCCGCCCCCGCCCTTTTTTTTTCTTTTTTATCGTTCAGGAACCGGGGGATGTTTATCGGCACCTGGGGTGATTTACGTTA‑GTGGTGATCGGGTTTTTTTTACAACTTTTTGAGTAGAAAAAAAAAATCCGCCGATCCt < 1:444995/143‑1 (MQ=255) ttCCATACTCGCTCCCGGCCCGCTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:408257/1‑151 (MQ=255) ccATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCg > 2:904746/1‑151 (MQ=255) agccgcGCCCCAGGCCCGTTCTTTTTTTTTTTCTCTTGAAGGTAACCGTGAATGTTTGT‑GGGAACGCGGTGTATTTTAGTTAGGTGGTAATAAGTTTTTTTGAACACTTGTGTGAGAGGAAAAAAAAAAACCGCCGATCCTGTCCAc < 1:836754/142‑1 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc > 2:733134/1‑151 (MQ=255) cGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt > 2:119470/1‑151 (MQ=255) cGCCCCCGGCCCTTTTTTTTATTTTTTATCGTTCAGGTACCCGTGGTTGTTTATCGGCACCTGGGGTGTTTTACGTTA‑GTGGTAATCAGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa < 2:817774/149‑1 (MQ=255) ccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg > 1:132319/1‑150 (MQ=255) cTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc > 1:539712/1‑150 (MQ=255) ttttCATTTTTATCGTTCAGGTACCCGTGGATGTTTATCGGCAACTGGGGTGATTTACTTTA‑GTGGGGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc < 2:738275/149‑1 (MQ=255) tttCATTTTTATCGTCCAGGTACCCGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTA‑GTGGGGATCAGGTTTTTTTTACACCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc < 2:722930/148‑1 (MQ=255) tttCATTGTTATCGTTCAGGTACCCGTGGATGTTTACCGGCACCTGGGGTGATTTACGTTG‑GTGGTGACCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa < 2:519641/150‑1 (MQ=255) cATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc > 2:988810/1‑151 (MQ=255) ttGTTATCGTTCAGGTAACGGTGGATGGTTTTCGGCAACTGGGGTGATTTACTTTA‑GTGGTAATCGGTCTTTTTTTACAATTGTTGGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc < 1:309929/149‑1 (MQ=255) ttGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCACCTGGGGTGATTTACGTTA‑GTGGTGATCAGGTTTTTTTTACAACTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAcc < 2:160726/149‑1 (MQ=255) tATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt > 2:813322/1‑151 (MQ=255) aTCTTTCAGGTACCGGTGGATGGTTATCGGCAACTGGGGTTATTTACGTTA‑GGGGGGATCAGGTTTTTTTTACATCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt < 2:241539/150‑1 (MQ=255) aTCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc > 2:612568/1‑151 (MQ=255) aGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAggt > 1:688332/1‑149 (MQ=255) ggTACCGGTGGATGTTTATCGGCACCTGGGGTGTTTTACGTTA‑GTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTaa < 2:333213/147‑1 (MQ=255) accggtggtgggTTATCGGCAACTGGGGTGTTTTACGTTA‑GTGGTGATCGGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt < 2:637985/141‑1 (MQ=255) cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg > 1:733122/1‑150 (MQ=255) tGGATGTTTATCGACAACTGGGGTTATTTACTTTA‑GTGGTGACCAGTTTTTTTTTACAACTTTTAGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt < 1:87576/150‑1 (MQ=255) tGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg > 2:179783/1‑151 (MQ=255) ggttgtttATCGGCAACTGGGGTAATTTACGTTA‑GTGGTGATCAGGTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg < 1:693452/144‑1 (MQ=255) ggATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTg > 1:42448/1‑150 (MQ=255) ggATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg > 2:900742/1‑150 (MQ=255) ggTTATCGGCAACTGGGGTGATTTACTTTA‑GTGGTAATCGGTTTTTTTTAACAACTTTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg < 2:749368/150‑1 (MQ=255) ggTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaa > 2:591377/1‑148 (MQ=255) ttttCGCCAACGCGGGTATTTGGCGTTT‑GTGGTAATCGGTCTTTTTTTACAATTTTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt < 1:749125/136‑1 (MQ=255) ttATCGGAAACGGGGGTGATTTACGTTA‑GTGGTGATCAGGTTTTTTTAACACTTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAg < 1:269951/148‑1 (MQ=255) tagcGGAAACTGGGGTGATTTATGTTA‑GTGGTGATCGGTTTTTTTTTACACCTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGaaa < 2:198268/143‑1 (MQ=255) aCTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc > 2:657833/1‑149 (MQ=255) gggTGTTTTACGTTG‑GTGGTGATCAGGTTTTTTTTACAACTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag < 2:991596/147‑1 (MQ=255) tttgggTTA‑GGGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc < 1:97070/144‑1 (MQ=255) tttAGGTTA‑GTGGTGATCGGGTTTTTTTTACATCTGTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc < 2:889561/149‑1 (MQ=255) tttACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt > 2:463337/1‑151 (MQ=255) ttACTTTA‑GTGGTGATCAGGTTTTTTTTCCACCTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt < 2:105316/150‑1 (MQ=255) tACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt > 1:765782/1‑151 (MQ=255) ttA‑GTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa < 2:81771/149‑1 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 2:218797/1‑151 (MQ=255) tgATCAGGTTTTTTTTACATCTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGc < 2:640342/150‑1 (MQ=255) tgATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTcc > 1:875602/1‑151 (MQ=255) aTCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGc > 2:761687/1‑151 (MQ=255) aGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg > 1:384298/1‑151 (MQ=255) tttttttACATCTTTTGGAGAAGAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGAc < 2:227993/149‑1 (MQ=255) ttttttACACCTGTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGt < 2:579972/150‑1 (MQ=255) ttttttACAACTG‑TGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 2:482770/150‑1 (MQ=255) ttttcccgaCTGTTTGAAAAAAAAAAGAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt < 2:408346/128‑1 (MQ=255) ttACATCTTTTGGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTg < 1:354417/150‑1 (MQ=255) ttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt > 2:63673/1‑150 (MQ=255) aagaaaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCaa < 2:688264/150‑1 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt > 1:992488/1‑150 (MQ=255) | CCGACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTA‑GTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT > NC_000913/1824907‑1825205 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |