Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,356,201 | 0 | T | G | 51.4% | 40.2 / 57.5 | 36 | T304P (ACC→CCC) | sapA | antimicrobial peptide transport ABC transporter periplasmic binding protein |
Reads supporting (aligned to +/- strand): ref base T (0/17); new base G (18/0); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.20e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-06 |
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGG > NC_000913/1356054‑1356344 | cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGtt < 1:497335/151‑1 (MQ=255) aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCa < 1:499589/150‑1 (MQ=255) cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATa < 1:450363/151‑1 (MQ=255) ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGAGTTAAA‑TGCCAGATAg < 2:631697/151‑1 (MQ=255) cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAgg < 2:169419/151‑1 (MQ=255) cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCg < 2:670852/151‑1 (MQ=255) gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAA‑TGCCAGATAGGCg > 2:639777/1‑150 (MQ=255) aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAA‑TGCCAGCTaag > 1:529716/1‑145 (MQ=255) aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGAc < 1:32741/151‑1 (MQ=255) gagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCa < 2:875417/151‑1 (MQ=255) gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTa < 1:784506/151‑1 (MQ=255) gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAAc < 2:136285/151‑1 (MQ=255) cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCaa < 2:79372/151‑1 (MQ=255) aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGTGGGGGTTAAA‑T‑CCAGAAAAGGGGCGTTCATTCCCGGACGCAGCGGTAAACGCAAgcgc > 1:529252/1‑148 (MQ=255) cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtc < 1:22937/151‑1 (MQ=255) cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAA‑AGCCAGAAAGGGGACGTTTATTCCAAGACGCAGCGTTTAAAACAACAGCggg > 1:311168/1‑149 (MQ=255) taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAA‑TGCCAGATAGGCGGCGTTTATCCCCCGACCCCGCGGTTAAAGCAAACGCGGGGCGTCAc > 1:304370/1‑151 (MQ=255) gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGGTTAACGCAAGCGCGGGGGGTCACGc > 2:743356/1‑151 (MQ=255) ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑aaa < 1:399686/151‑1 (MQ=255) tCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAA‑‑GCCAGAAAGGGGGCGTTTCTCCCCGGACGCCGAGGTTAACACAAACGGGGGGCGGCACGC‑AAAAAAGATa > 1:824691/1‑150 (MQ=255) cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAA‑AGCCCGATAGGGGAAGTTCATCCCAGGACGCAGCGCTTAAAGCAAGCGCGGGGCGGCACGC‑AAAATGGAt > 2:496504/1‑148 (MQ=255) gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGC‑AAAATGGATAgctgg > 2:2431/1‑149 (MQ=255) ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAA‑ATGCAAAAAGGCGACGGTTCTCCCCGGACGCCACGGTAAAAGCAAACGCGGGGCGGCAAGCAAAAATGGATGGcgggcggg > 1:547040/1‑151 (MQ=255) aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATTGCCA‑AGACGGGGCGGTCATCCCCGGAACCAGCGGTAAAAGCAAGCGCGGGGCGTGCCGC‑AAAAAAGGTAGCCGGCCGGGTGCaggcccggcc > 1:962754/1‑150 (MQ=255) aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAA‑TGCCAGATAGGGGACGTTCATCCCCGGACGCAGCGGTTAACGCAAGCGCGGGGCGTCACGC‑AAAAAGGGTAGCTGGCTGGGAGCagggcaggcc > 2:889799/1‑150 (MQ=255) gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAA‑TGCCAGATAAGGGGGGTTCATCCCAGGACGCAGCGGTAAAAGCAAGCGGGGGGCGGCCCGC‑AAAAAGGGTAGCTGGGTGGCAGcaggccagg > 1:19158/1‑145 (MQ=255) gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGa < 1:949699/150‑1 (MQ=255) cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAA‑TGCCAAAAAGGCGCCGTTCATCCCAGGACGCCGCGGTTAACGCAAAGGCGGGGGGTGCCGC‑AAAATGGGAAGCTGGGCGGCAGCAGGCCAGGGCaaaa > 2:381318/1‑149 (MQ=255) aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAA‑AGCCCGGTAGGGGACGTTTATCCCCGGGCGCAGCGTTTAACGCAAACGCGGGGGGTCACGC‑CAAAAGGATAGCTGGCTGGCAGCAAGGCAGGCCAGAACggc > 2:721671/1‑145 (MQ=255) ggACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTTAA‑TGCCAGATAGGGGGCGTTTATCCCCGGACGCCGCGGTAAAAGCAAACGCGGGGCGGCACGC‑AAAAAGGGTAGGTGGGTGGGAGCCGGGCCGGCCCGAACGGCGCCATcccc > 1:910020/1‑151 (MQ=255) tattTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGtt < 2:135398/150‑1 (MQ=255) tattTAGCGGCGG‑TTTGGCGGGGGTAAA‑TGCCAGATAAGCGGCGTTCATCCCCCGACGCAGCGGTTAAAGCAAACGCGGGGCGGCACGC‑AAAAAGGGTAGCTGGCTGGGAGGCGGCCAGGCCCGAACGGCGCATTCCCCCGGCAgggggtt > 1:227308/1‑146 (MQ=255) tAGCGGCGG‑TTTTGTGGGGTTAAA‑TGCCAGATAAGCGACGTTTATCCCAGGAAGCAGCGTTAAACGCAAGCGCGGGGCGTCACGC‑AAAAAAGGAAGCTGGCTGGCAGCAGGGCCGGGCAGAAAGGCGCATTCCCCGGTCAGGAGGTTCGa > 1:806905/1‑150 (MQ=255) cggcgg‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa < 2:466347/151‑1 (MQ=255) cggcgg‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa < 2:633504/151‑1 (MQ=255) tGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg < 1:834681/150‑1 (MQ=255) | CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAA‑TGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGC‑AAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGG > NC_000913/1356054‑1356344 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |