| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 366,526 | 0 | T | G | 54.8% | 25.2 / 44.4 | 31 | T329P (ACC→CCC) | lacI | lactose‑inducible lac operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base T (4/10); new base G (17/0); total (21/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.26e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-06 | |||||||||||
AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAA > NC_000913/366378‑366668 | aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGTGTATTGGGCGCCC‑gggggg > 2:178815/1‑151 (MQ=255) aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑Gggtggt < 2:390953/151‑1 (MQ=255) cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTTGGCGCCC‑GGGGGGttt > 2:83047/1‑150 (MQ=255) tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCC‑GGGGGGgtttt > 2:555537/1‑151 (MQ=255) tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGCGTATTGGGGGCCC‑GGGGGGgtttt > 1:477802/1‑151 (MQ=255) gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGTATTTGGGGCCC‑GGGGGGGttttttttt > 2:376633/1‑150 (MQ=255) aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGTGTATTTGGCGCCA‑GGGGGGGTTTTCTTTTTAcc > 2:263097/1‑150 (MQ=255) aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTATTGGGGGCCA‑GGGGGGGTTTTTTTTTTAcc > 2:248672/1‑149 (MQ=255) gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGGATTGGGGGCCA‑GGGGGGGTTTTCTTTTTACCAGt > 2:401215/1‑150 (MQ=255) aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCa < 1:44877/151‑1 (MQ=255) caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑CTTTCGCTTTGTGGTCCC‑GGGTGGGTTTTTTTTTTACCAACGAGAGAGGCGACa > 1:113463/1‑148 (MQ=255) aTTAATTGCGTTGCGCTCACTGCCCGCTTGCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGc < 1:652554/149‑1 (MQ=255) aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt < 2:240993/151‑1 (MQ=255) aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg < 1:434054/151‑1 (MQ=255) tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGGTTGGGCGGCA‑GGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGGTTGc > 2:605220/1‑150 (MQ=255) cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑ATTGCGGAATGGGGGCCA‑GGGGGGGTTTTTTTTTTTCCCGAGAGGCGGGGCACAGCTGAtttgcc > 1:694704/1‑146 (MQ=255) ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGGTTTGGGCGCCC‑GGGGGGGTTTTTTTTTTCCCCGGGAGACGGGGAACAGCTGATTGcccct > 2:674340/1‑148 (MQ=255) ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTTCGTATTTGGGGCCC‑GGGGGGGTTTTTTTTTTTCCCGTGAGAAGGGGAACAGCTGATTGgcc > 1:346993/1‑148 (MQ=255) gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGGGGGTTTTTTTTTTACCAGGGAGGCGGGCAACCGCTGATTGCCCTTCAcc > 1:153190/1‑150 (MQ=255) cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCCGCTTTGGG‑GCCC‑GGGGGGGTTTTTTTTTTACCCGAGGGAAGGGCGACAGCCGGTTGGCCCTCTcccc > 1:156369/1‑149 (MQ=255) cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt < 1:376597/151‑1 (MQ=255) ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt < 1:170100/151‑1 (MQ=255) gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca < 1:172845/151‑1 (MQ=255) ccTGTCGCGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTCCGTCTTTGGTGCCA‑GCGGGGGTTTTCTTTTCACCCATGAGACCGGGAACAGCCGGATGCCCCTCACCCCCCAGCCCCGAGAGAGGTggagca > 2:554774/1‑147 (MQ=255) cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGGATTTGGCGGCAGGGGGGGTTTTTTTTTTCA‑CAGTGAGAAGGGGAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGTTGCAGCAAACGGGCCACGCTGGgtt > 1:109777/1‑149 (MQ=255) ggCCAACGCGCGGGGAGAGGCGG‑TTTGCGGATTGGGCGCCA‑GGGGGGTTTTTCTTTTCACCCGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCTGGCCCCGAGAGAGTTGCAGCAAACGGTCCACGCTGGGTTGGCCCCGCAGGCg > 2:420019/1‑150 (MQ=255) aaCGCGCGGGGAGAGGCGG‑TTTGTGGATTGGGGGCCA‑GGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCTCGCCCCGGGAGAGTTGGAGCAAGAGGGCCCCGCCGGGTTGCCCCCGCAGGCGgaaa > 1:585075/1‑150 (MQ=255) gggAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg < 1:274781/151‑1 (MQ=255) gg‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt < 2:360015/150‑1 (MQ=255) gg‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt > 2:235494/1‑151 (MQ=255) gCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg > 1:575474/1‑151 (MQ=255) gTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg < 1:514403/151‑1 (MQ=255) gcCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATaa > 1:412186/1‑150 (MQ=255) | AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCCA‑GGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAA > NC_000913/366378‑366668 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |