Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 366,526 | 0 | T | G | 56.5% | 14.3 / 24.0 | 23 | T329P (ACC→CCC) | lacI | lactose‑inducible lac operon transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base T (4/6); new base G (13/0); total (17/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.08e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-03 |
AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCA‑AGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG > NC_000913/366379‑366674 | aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCCggggggg > 2:477365/1‑150 (MQ=255) ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCGCCCGGGGGGGTTTTTTTTTCa > 1:83396/1‑151 (MQ=255) gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCCGGGGGGTTTTTCTTTTTACCAGTg > 1:280587/1‑151 (MQ=255) cTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc < 1:442951/151‑1 (MQ=255) caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGTTTGGGGGGCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAAA‑Ag > 1:16878/1‑149 (MQ=255) caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAGTGGGCGCCCGGGGGGGTTTTCTTTTTACCAGGGAGAggggggaaagg > 1:213428/1‑143 (MQ=255) cgcttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTTTTTTTACCAGTGAGAAGGGCAAA‑AGCTGAt > 2:279782/4‑151 (MQ=255) aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTg > 2:236789/1‑151 (MQ=255) cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCAc < 2:247882/150‑1 (MQ=255) cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGTGGGTTTTCTTTTTACCAGGGAGACGGGCAAC‑AGCTTATTTGCCCTCCCCGc > 2:179743/1‑149 (MQ=255) cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAAC‑AGCTGATTGCCCTTCACCGcc > 1:351094/1‑150 (MQ=255) cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGTGGTTTTTCTTTTCACCAGGGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTgg > 2:357644/1‑151 (MQ=255) cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTTTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGGCCTTCACCGCCTgg > 1:541092/1‑151 (MQ=255) gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGTTTTTCTTTTTACCCGGGAGAGGGGCGAC‑AGCTGATTTCCCTTCCCCGCCTGGc > 1:91318/1‑150 (MQ=255) gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGATTTGGGGGCAGGGGGGGTTTTTTTTTTAACCGGGGGACGGGGAACAAGCGGTTTGCCCtccccccccggc > 1:314124/1‑139 (MQ=255) gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCCCCTGGCCCTga > 2:500579/1‑151 (MQ=255) cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTgag < 2:83398/151‑1 (MQ=255) cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTgag < 1:326850/151‑1 (MQ=255) cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGGGGTTTTTTTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca > 2:426504/1‑151 (MQ=255) gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTATTGGGCGCCAGGGGGGTATTTCTTTTCCACAGGGAGGAGGGCAAC‑AGCTGATTGCCCTTCACCCCCTGGCCCTGAGGGAAGTGCAGCA‑AGCggg > 1:54179/1‑149 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGGGAGAGTTGCAGCA‑AGCGGTCCACGCTGGTTTg > 2:330399/1‑151 (MQ=255) gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCAGGGGGGGTTTTTTTTTCCCCCGGGGGACGGGCAAC‑AGCTGATTGGCCTCCACCGCCCTGCCCCGGGAGAAGTGCAGCACAGCGTCCCACCCGGGttg > 1:112237/1‑150 (MQ=255) cGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCA‑AGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg < 2:16837/150‑1 (MQ=255) gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCA‑AGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat < 2:532848/150‑1 (MQ=255) gtggtTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCA‑AGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAAATGAg > 1:226140/1‑150 (MQ=255) | AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAC‑AGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCA‑AGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG > NC_000913/366379‑366674 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |