Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 366,526 | 0 | T | G | 50.0% | 48.0 / 29.2 | 50 | T329P (ACC→CCC) | lacI | lactose‑inducible lac operon transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base T (5/20); new base G (25/0); total (30/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC > NC_000913/366379‑366675 | aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑ggggggg > 2:328146/1‑150 (MQ=255) gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTG‑GGATTTGGGGCCC‑gggggggt > 1:954399/1‑149 (MQ=255) gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGCGCCC‑gggggggt > 1:687568/1‑149 (MQ=255) gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGT‑GTATTGGGCGGCC‑GGGGGGGttttcttttt > 2:1196149/1‑150 (MQ=255) ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGTTTTGGGGCCCG‑GGGGGGtttttttttt > 1:75068/1‑150 (MQ=255) ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGCTTTGGGGGCCG‑GGGGGGtttttttttta > 1:1261256/1‑150 (MQ=255) tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTC‑GGCTTTGGCGGCC‑GGGGGGGttttctttttca > 1:333850/1‑147 (MQ=255) tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGC‑GTGTTGGGGGCCG‑GGGGGGGTTTTCTTTTCACCAgg > 1:437915/1‑150 (MQ=255) agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTg < 1:99953/150‑1 (MQ=255) gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGG‑GGTTTTGGGGCCA‑GGGGGGGTTTTTTTTTTACCCGTgag > 1:795688/1‑151 (MQ=255) tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTgaga < 2:1017862/151‑1 (MQ=255) gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGC‑CGCTTGGGGGCGACGGGTGG‑TTTTCTTTTTAccccagagaa > 2:994163/1‑143 (MQ=255) gCTAACTCCCATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGAGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTAACCCGTGAGACg < 2:1003694/150‑1 (MQ=255) gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTTTTGGGGGCCA‑GGGGGGGTTTTTTTTTCTcccgggggaggg > 1:772209/1‑141 (MQ=255) aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTT‑GGCTTGTGCGCCA‑CGGTGGGTTTTCTTTTTACCACTGCGGCgggg > 2:489003/1‑149 (MQ=255) tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCA‑GGGGGGTTTTTTCTTTTACCAGGGGGACGGGCAACAG‑CTGAt > 2:556512/1‑151 (MQ=255) tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGCGCCA‑GGGGGGTTTTTCTTTTCACCAGGGAGACGGGCAACAG‑CTGAt > 2:264559/1‑151 (MQ=255) gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCA‑GGGGGGGTTTTTTTTTTTCCAGTGAGAAGGGCAACAG‑CCTATTGcc > 2:543405/1‑150 (MQ=255) gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGC‑GTGTTGGGGGCCC‑GGGGGGGTTTTTTTTTCACCCGTGAGACGGGGAACAG‑GTGGTTggc > 2:740108/1‑148 (MQ=255) gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTC‑GTATTTGGCGCCA‑GGGGGGGTTTTTTTTTTACCAGGGAGGCGGGCCACAG‑CTGATTGCCCCTct > 1:805807/1‑149 (MQ=255) cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCAc < 2:638416/150‑1 (MQ=255) cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTG‑GTATTTGGGGCCA‑GGGGGGGTTTTTTTTTTCCCCGGGGGACGGGCAAAAG‑GTGGGTGccccttcccc > 1:742634/1‑143 (MQ=255) cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTG‑GTTTTTGGGGCCC‑GGGGGGGTTTTTTTTTTACCAGGGAGAAGGGGAACAG‑CTGGTTGGCCTTCACcccc > 1:1122833/1‑150 (MQ=255) cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTgag < 2:960349/151‑1 (MQ=255) cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGGGGGTTTTCTTTTTCCCAGTGGGAAGGGCAACAG‑CTGGTTGCCCCTTCCCGCCCGGccccggg > 1:666734/1‑147 (MQ=255) tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt < 2:864037/148‑1 (MQ=255) aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGGGCCC‑GGGGGGGTTTTTTTTTCACCAGGGAGACCGGCCACAG‑ATGAGTGCCCCTCACCGCCTGGCCCCGAGGGAGGTGCCGCaa > 2:802497/1‑151 (MQ=255) ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGT‑GGATTTGGGGCCA‑GGGGGGGTTTTCTTTTCACCACGGGGGAGGGCGAAAGACGGATTGCCC‑TCCCCGCCCTGCCCCGTGAGAGTTTCGGCCAGCg > 1:414624/1‑151 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTT‑GTATTGGGGGCCA‑GGGGGGGTTTTCTTTTCACCCGGGAGAAGGGCCAAAG‑CTGATTGGCCCTCACCGCCTTTCCCCGTGGGAGGTGGAGCACGGGGTCCCCCCTGGGTTg > 1:873520/1‑151 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGG‑GGATTGGGCGCCC‑GGGGGGTTTTTCTTTTTACCAGTGGGACGGGCAACAG‑CTGATTGCCCCTCCCCCCCTGGCCCTGAGGGAGTTGCAGCAAACGGGCCACCCCGGGTTg > 1:725620/1‑151 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGGTTTGGGGCCC‑GGGGGGGTTTTTTTTTTACCAGGGGGAAGGGCAAACG‑CTGATTGTCCTTCCCCCCCCGGGCCCGGGGGAGGTGCAGCCAACGGTCCACGCTGGTTTg > 1:278889/1‑151 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGGGCCA‑GGGGGGTTTTTTTTTTCACCAGGGAGACGGGGAACAG‑CTGATTTCCCCTCCCCCCCCGGCCCTGGGGGGGGTGGAGCAAGCGGGCCAccccgggttt > 1:735259/1‑142 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTGGGGGCCA‑GGGGGGGTTTTTTTTTTACCCGGGAGACGGGGAAAAG‑CTGAATGGCCTTCACCCCCCGGCCCTGCGGGAGGTGCAGCAAGCGGGCCCCGCTGGGTTg > 1:834801/1‑151 (MQ=255) ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGG‑GGATTGGGGGGCC‑GGGGGGGTTTTTTTTTCTCCACCGGGACGGGGAACAG‑CTGGTTGCCCTTCCCCCCCCTGCCCCCAGAGGGGTGCAGCAAGCCGCCCCCGCCGGGTTGcc > 2:893192/1‑150 (MQ=255) gCCAACGCGCGGGGAGAGGGGGTTGGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa < 1:781186/150‑1 (MQ=255) gCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCC‑GGGGGGGTTTTTTTTTTACCCGGGAGACGGGCAACAG‑CTGGTTGCCCTTCCCCGCCCGGCCCTGAGGGAGTTGCAGCAAACGGGCCACGCTGGTTTGCCCCACCAGGCgga > 1:189500/1‑149 (MQ=255) cgGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt < 2:1122407/151‑1 (MQ=255) gggAGAGGCGGTTTTT‑GTATTTGGCGCCC‑GGGGGGGTTTTTTTTTCACCCGTGAGAAGGGGAAAAG‑CTGATTGCCCCTTACCGCCTGGCCCCGAGAGAGTTGCAGCAAGCGGGTCCCGCTGGTTTGCCCCAGCAGGGGAAAATTCTGtttt > 2:865771/1‑150 (MQ=255) agGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg < 2:771893/150‑1 (MQ=255) gTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa > 2:903264/1‑151 (MQ=255) gTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa > 1:485810/1‑151 (MQ=255) ttGC‑GTATTGGGCGCAA‑GGGTGGTTTTTCTTTTCTCCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc < 2:950736/151‑1 (MQ=255) c‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg < 1:829218/150‑1 (MQ=255) tATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg > 1:406436/1‑151 (MQ=255) gcgcCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat < 1:828688/150‑1 (MQ=255) gcgcCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat < 2:518352/150‑1 (MQ=255) gcgcCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat < 2:42403/150‑1 (MQ=255) ccA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTTACCGCCTGGCACTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATAAAAc < 2:1164990/150‑1 (MQ=255) ccA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCCCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc < 2:687259/150‑1 (MQ=255) ccA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc < 1:1196643/150‑1 (MQ=255) ccA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa < 2:954050/151‑1 (MQ=255) ccA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa < 1:729202/151‑1 (MQ=255) ccA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa < 1:676868/151‑1 (MQ=255) ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt < 2:156845/148‑1 (MQ=255) tggtTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc < 2:660553/150‑1 (MQ=255) | AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC > NC_000913/366379‑366675 |
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: ‑ |