| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 17,000 | 0 | C | T | 45.5% | 37.2 / 20.6 | 44 | intergenic (‑40/‑489) | mokC/nhaA | regulatory protein for HokC, overlaps CDS of hokC/sodium‑proton antiporter |
| Reads supporting (aligned to +/- strand): ref base C (18/6); new base T (0/20); total (18/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.57e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAT > NC_000913/16853‑17146 | ggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCttt < 1:187062/151‑1 (MQ=255)ggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCttt > 1:378684/1‑151 (MQ=255) gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtc < 2:392042/150‑1 (MQ=255) tACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTctct > 1:342280/1‑151 (MQ=255) gACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATg > 2:252427/1‑151 (MQ=255) aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGc > 1:230056/1‑151 (MQ=255) tcGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGagac > 1:543208/1‑151 (MQ=255) ccTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGagac > 1:515864/1‑148 (MQ=255) ttATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGagacagac > 2:431819/1‑150 (MQ=255) gctgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCt > 1:92060/1‑151 (MQ=255) ctgctTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTc > 2:348428/1‑151 (MQ=255) ctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACccgcc < 1:186663/149‑1 (MQ=255) tcCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCtattat > 2:78120/1‑151 (MQ=255) ttGCCCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGAATATGGGAGGCCTCGGGTTGATGAAAAAATATCCCTCGGGGTTTTTCTCTATCTCCCTTTCAGCAAATGCC‑TGAGACAGACAGCCTCAAGCCCCCGCCGCTATTatat < 2:543177/150‑1 (MQ=255) tgcccTTACGGTAAGAAAGAGGCACTCTACATTTGTTCAGAATATAGGGGGCCTGGGTTTGATGGAAAAATACCCCTGGGGCTTTTTCTCTACCTCCCTTTCAGCAAATCCC‑TGAGACAGACAGCCTCAAGCCCCCGCCGCTATTatat < 2:109382/146‑1 (MQ=255) ggTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCttt > 2:286063/1‑150 (MQ=255) ggTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCttt > 1:291632/1‑150 (MQ=255) gggcccctctACATGTTTTCAGCATATGGGAGGCCTCGGGTTAATGGTAAAATATCCCTCGGGGTTTTTCTTTATCTCCCGTTCAGCTAATCCC‑TGAGACAGACACCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCAtttt < 2:147475/145‑1 (MQ=255) gctctatACGTGTGTTCATCATATGGGGGGCCTGGGTTTGGTGGAAAAATATCACTCGGGTTTTTTCTTTATCCGCCGTTCA‑CCAATGCCATGAGAAAAACAGCCTCAACCAGCCGCCGCTATTATATCGCTCTCTTTAACCCAtttt < 1:545926/142‑1 (MQ=255) ctctACTTGTTCTAAACATATGGGAGGCCGGGGTTGGATGGAAAAATACCCCCCGGGCCTTCTTTCTACCTCCCGTTCAGCTAATCCC‑TGATAGAGACACCCCCGAGCGCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGttt < 2:531490/149‑1 (MQ=255) aCATGTGTTCAGCATATAGGGGGCCTGGGGTTGATGGAAAAATATCCCTGGGGTTTTTTCTTTATCTGCCTTTCAGCTAATGCC‑TGAGACAGACAGCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCg < 1:515457/150‑1 (MQ=255) aaaatatAGGGGCCCGGGGGTTGATGATAAAATACCCCTCGGGGTTTTTCTCTATCTGCCGTCCATCTAATGCT‑TGAGACAACCACCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTa < 1:545928/143‑1 (MQ=255) aGCATAAGGGGGGCCTGGGGTTGATTAAAAAATTTCTTTGGGGGCTTTTCTCTCCCCGCCGTTCCGCAAATCCT‑TGAGAAAGACACCCTCATTCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATc < 1:326446/149‑1 (MQ=255) atatAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATc > 1:455171/1‑146 (MQ=255) tgggaggCCTGGGTTTGATGTAAAAATACCCCCCGGGTTTTTTCTCTATCTCCCGTTCAGCTAATCCC‑TGAGACAGACACCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATc < 2:455163/141‑1 (MQ=255) ggggcccgggggTTGATGGAAAAATATCACTGGGGTTTTTTCTCTATCTGCCTTACAGTTAATCCC‑GGAAAAAGACAGCCACAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAg < 1:369826/138‑1 (MQ=255) ggggcccTGGGTTTGATGAAAAAATTCCCCGCGGGGTTTTTCTCTCTCTCCCGTTCACCTATTCCG‑TGAAACACACACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGaa < 2:213372/141‑1 (MQ=255) ggCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc > 2:46893/1‑151 (MQ=255) ggtttGATGAAAAAATACCCCGCGGGTTTTTTCTTTATCGCCCTTAAAGTAAATGCC‑TGAGACAGACCCCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc < 2:140066/142‑1 (MQ=255) gggTTGATGGAAAAATATCCCTCGGGTTTTTTCTCTATCTCCCTTTCAGCAAATGCC‑TGAGAAAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttt < 2:72859/151‑1 (MQ=255) gTTGATGGTAAAATATCACTCGGGACTTTTATCTATCTGCCGTTCAGCTACTGCC‑TGAGACAGTCAGCCTCAAGCAACCGCCGCTATTATATCGCTCTCCTTCACCCATTTTGTTTTATCGAGTCTAATTCTGAAGACGCCTCGTAt > 2:437256/1‑147 (MQ=255) aTGGAAAAGTATCTCGCGGGTTTTTTCTCTTTCGGCCGTTCACCTACTGTC‑GGATAAAGCCAGCCTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAtttt < 2:560542/146‑1 (MQ=255) tgaaaaaaTATCCCGGGGGTTTTTTCTCTATCTGCCTTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtgt < 2:414040/144‑1 (MQ=255) ggaaaaaTACCCCTCGGGGTTTTTCTCTTTCTCCCGTTCAGCTAATCCC‑TGAGACAGCCAGCCAAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTtgt < 1:181172/144‑1 (MQ=255) aaaaTCCCACGCGGGGTTTTTCTCTCTCTCCCTTTCATCAAATCCC‑TGAAACACACCGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttttt < 2:32667/142‑1 (MQ=255) aaaaTACCCCTGGGGGTTTTTCTCTATCTGCCGTTCAGCTAATGCC‑GGAGACAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTgg < 1:363047/146‑1 (MQ=255) ataccccggggggTTTTTTTCTACCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAttt < 1:443568/142‑1 (MQ=255) tatCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAtttt > 2:313518/1‑151 (MQ=255) atCACTCGGGTCTTTTCTCTACCTCCCTTTCAGCTAATGCC‑TGAGACAGACACCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAtttt < 2:251010/150‑1 (MQ=255) atCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAttttt > 1:148577/1‑151 (MQ=255) cccTCGGGTTTTTTCTCTATCTGCCTTTCATTTAATCCC‑TGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGTGTAAttattt < 1:83393/148‑1 (MQ=255) ggggTTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATg < 1:348435/149‑1 (MQ=255) ggggTTTTTCTCTACCTGCCTTTCAGCTAATGCT‑TGAGAAAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt < 2:90346/151‑1 (MQ=255) ggggCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt > 2:159602/1‑151 (MQ=255) | GGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCC‑TGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAT > NC_000913/16853‑17146 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |