Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 366,519 | G→A | A331V (GCG→GTG) | lacI ← | DNA‑binding transcriptional repressor LacI |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 366,519 | 0 | G | A | 93.8% | 37.4 / ‑7.2 | 16 | A331V (GCG→GTG) | lacI | DNA‑binding transcriptional repressor LacI |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (2/13); total (3/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.88e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCA > NC_000913/366412‑366759 | caatcgcgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGG‑GTTTTTTTTTTACCAGAGAGACGGGCAACCGATGATTGCCCCTC‑ACCCCCTCGCCCTGCGGGAGGTGTGGCAAGGCGG‑CCCCGCTGGGTTGGCCCCGCCAGG‑GAAAATCCAGTTTGGTTGTGGTTAACgggggg > 2:325992/8‑251 (MQ=255) cgcgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑Tttttcttttc < 1:200506/121‑1 (MQ=255) cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc < 1:157271/251‑1 (MQ=255) gCCCGCCTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGACCAGGGGGGGTTTTTTTTTTC‑CCAGTGAGACGGGCAAAAGCTGATTGCCCTTCCACC‑CCCGGCCCCTAGGGAGGTGCAGCAAG‑CGGTCCACGCTGG‑GTTGCCCCAGCAGGCGGAAATCCTGGTTTATGGTGGTTTACGGCggggtgttacataatgtgtattct > 1:120365/1‑230 (MQ=255) cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGTTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg < 1:167617/224‑1 (MQ=255) ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAGCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGt < 1:357201/251‑1 (MQ=255) tacaggcGGGAAACCTGTCGTGCCAGCTGCATTAATGAACCGGCCAACGCGGGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGCCGGGCAACAGCTGATTGCCCTTC‑TTCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCTAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTAtcgt < 1:361901/244‑1 (MQ=255) ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGGGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCGCCAGCAGGCGAAAATCCTGTTTTATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTAtcgtc < 1:381993/250‑1 (MQ=255) gggAAACCGTTCGGGCCCGCTGCATTAATGAATCGGCCAACGCGCGGGGAGGGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATccc < 2:225312/251‑1 (MQ=255) ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGTCTTGAGATAACATGAGCTGTCTTCGGTATCGTCGTATCCCa < 1:52818/251‑1 (MQ=255) aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCca < 1:247804/160‑1 (MQ=255) aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGcc < 1:161686/251‑1 (MQ=255) gcgcgGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGt < 1:103237/251‑1 (MQ=255) cGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCa < 1:213100/251‑1 (MQ=255) tGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg > 1:324221/1‑135 (MQ=255) tGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg < 2:324221/135‑1 (MQ=255) | ATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCA > NC_000913/366412‑366759 |
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: ‑ |