Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,056 | 0 | G | A | 43.2% | 39.2 / 14.0 | 37 | 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 (15/6); new base A (0/16); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.07e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT > NC_000913/1824910‑1825205 | aCATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaa < 1:862559/151‑1 (MQ=255) aCATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaa > 1:247567/1‑150 (MQ=255) ttttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga > 1:592483/1‑150 (MQ=255) tttACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaagaa > 2:777332/1‑150 (MQ=255) cATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgcc > 2:325955/1‑151 (MQ=255) tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg > 1:432835/1‑150 (MQ=255) aaGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg > 1:843893/1‑151 (MQ=255) tATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc > 2:652966/1‑151 (MQ=255) aTTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTcc > 1:476216/1‑150 (MQ=255) ttCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 2:80399/1‑151 (MQ=255) ccATACTCGCTCCCGGCCCGTTTTTTTCATTGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCAACTGGGGTGTTTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCg < 1:855471/151‑1 (MQ=255) cTCGCCCCCGGCCCGTTTTTTTCATTTTTATCGTTCAGTTACCCGTGGATGTTTATCGGCAACTCGGGTGTTTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTTTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCa < 1:624888/148‑1 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:656186/1‑151 (MQ=255) cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCa > 1:57992/1‑150 (MQ=255) tttttttCATTTTTATCTTTCGGGTAACCGTGGTTGTTTATCGCCAACTCGGGTGTTTTCC‑TTTAGTGGTGATCGGTTTTTTTTCACAACTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 2:720683/143‑1 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg > 2:11097/1‑150 (MQ=255) tttCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa > 1:340074/1‑149 (MQ=255) ttATCGTTCAGGTACCGGTGGTTGGTTACCGGCAACTGGGGTGATTTAC‑GTTAGTGGTGATCGGTTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc < 1:383133/149‑1 (MQ=255) tATCGTTCGGGTAACGGGGGATGTTTATCGGCACCTCGGGTGATTTAC‑GTTAGGGGTGATCGGGTTTTTTTTACAATTTTTTGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg < 2:675441/150‑1 (MQ=255) ttaaGGTAACCGTGGATGTTTATCGGCAACTGGGGTGTTTTAC‑TTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGt < 1:812212/143‑1 (MQ=255) ttCAGTTACCCTGGGGTGTTTATCGGCAACTGGGTTGATTTAG‑TTTAGTGGTGATCAGGTTTTTTTAACAACTTTTTGAAAAAAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTc < 2:101706/147‑1 (MQ=255) ttCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt < 2:274255/148‑1 (MQ=255) ccGTGGATGGTTATCGGCAACTGGGGTGTTTTCC‑GTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg < 1:160744/151‑1 (MQ=255) gTGGATGGTTATCGGCAACTGGGGTGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACACCTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt < 1:580677/151‑1 (MQ=255) gTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg > 1:547729/1‑149 (MQ=255) tATCGGAACTTCGGGGGTTTTAC‑TTTAGTGGTGACCGGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt < 2:609190/150‑1 (MQ=255) acagacaactgggTTATTTACGGTTGTTGGTGATCGGGTTTTTTTAACATCTGTTGGGAAAGAAAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt < 2:337915/139‑1 (MQ=255) aaCTGGGGTATTTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:174638/149‑1 (MQ=255) aaCTCGGGTGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc < 1:313578/147‑1 (MQ=255) aaCTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc > 2:630623/1‑150 (MQ=255) tggggTGATTTAC‑GTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc < 1:487056/145‑1 (MQ=255) gggTGATTTAC‑GTTAGTGGTGATCGGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt < 2:27532/148‑1 (MQ=255) gggTGATTTAC‑GTTAGGGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag < 1:204408/147‑1 (MQ=255) ggTGTTTTCC‑GTTAGGGGTGATCATGTTTTTTTCACACCTGTTTGAGAAAAAAAATGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 1:162472/151‑1 (MQ=255) c‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtttt > 1:871212/1‑150 (MQ=255) gtgttgTTCGGGTTTTTTTTACAACTGTTTGAAAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt < 2:325547/151‑1 (MQ=255) gtggtgATCAGTTTTTTTTTACACTTTTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc < 2:790945/148‑1 (MQ=255) ggtgAGCAGTTTTTTTTTACACCTGTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa < 2:886304/144‑1 (MQ=255) ggCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg > 1:663283/1‑150 (MQ=255) ttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:138521/1‑151 (MQ=255) tttACAACTTTTTGAAAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt < 2:277419/150‑1 (MQ=255) ttggaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc < 2:651339/141‑1 (MQ=255) agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAtt > 1:752616/1‑151 (MQ=255) agaaaaaaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAAt < 2:705018/150‑1 (MQ=255) | ACATTTTACGCATATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTAC‑GTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATT > NC_000913/1824910‑1825205 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |