Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 291,942 | A→G | K148R (AAG→AGG) | yagJ → | CP4‑6 prophage; protein YagJ |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 291,942 | 0 | A | G | 95.2% | 129.6 / ‑2.2 | 42 | K148R (AAG→AGG) | yagJ | CP4‑6 prophage; protein YagJ |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (18/22); total (19/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01 |
CTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAA > NC_000913/291807‑292072 | ctctGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACCGGGAATCACTTATGATGAAGCCGAACAGTACAGGAtt > 1:350013/1‑140 (MQ=255) tGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGtt > 2:706377/1‑140 (MQ=255) ttGAACCATTGGGAAACCAAGGATTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTAAGTTCAAACTGAACAGTACAGGTTTGTTCg > 1:335796/1‑140 (MQ=255) ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCg < 1:170035/140‑1 (MQ=255) ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCg < 1:706377/140‑1 (MQ=255) ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCg < 1:611126/140‑1 (MQ=255) ccATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCCTATGTTGTAGCCGATCAGTACAGGATTGTTCGTCCCg > 1:398373/1‑140 (MQ=255) aTTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGAt < 1:563632/140‑1 (MQ=255) ttGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGAtt > 2:540894/1‑140 (MQ=255) ttGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGAtt > 2:549581/1‑140 (MQ=255) gggAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTc > 1:189743/1‑140 (MQ=255) gggAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTc > 1:760475/1‑140 (MQ=255) ccAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCGTATGTTGAGGCTGAACCGTACAGGAGTGTTGGTCCCGATTTCCTGTTc > 2:212329/1‑140 (MQ=255) gTTTTTGAGATCGAAACCAAACGTTAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGc > 2:648117/1‑140 (MQ=255) cAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCaaaa < 2:484872/140‑1 (MQ=255) gTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGtt < 2:751895/140‑1 (MQ=255) gTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGtt < 2:561677/140‑1 (MQ=255) tGAAGGATCTGCTTTCTTGTATTGTAATCCAAACTCCACAGGGCAGTCATCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTtg < 2:398373/140‑1 (MQ=255) ggATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGa < 2:189743/140‑1 (MQ=255) gATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGAc < 2:650337/140‑1 (MQ=255) tGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTgg < 2:762424/140‑1 (MQ=255) tGCTTTCTGGTATCGTAATCCAAAGTACACAGGTCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTgg < 2:350013/140‑1 (MQ=255) ttCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGa < 1:692602/140‑1 (MQ=255) tATCGTAATCCACAGTAAACAGTGCAGTCGTCACTGGGAGTCGCTTATATTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCcaca < 2:335796/140‑1 (MQ=255) tCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATa > 1:163309/1‑140 (MQ=255) cacaGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAGAATGGTTATGGGCCGGGTAGAGCCACATAGCCTACAt > 1:235725/1‑140 (MQ=255) acaGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATc < 1:549581/140‑1 (MQ=255) aGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAGATGATTGCGGACCTGGGAGATCCACATAGCCTACATCTg > 2:289538/1‑140 (MQ=255) tACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGc > 1:670399/1‑140 (MQ=255) cacaGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTg < 1:109026/140‑1 (MQ=255) acaGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGa > 1:500990/1‑140 (MQ=255) ggCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCt > 1:767948/1‑140 (MQ=255) aGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTg > 1:138015/1‑140 (MQ=255) tcgtcACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCAGGTAGATTGACATAGCATACATTTGGCTGAGGCTTTGcc > 2:394560/1‑140 (MQ=255) gggAATCGCTTATGTTGAAGCTCAATAGTGCAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTgg < 1:289538/140‑1 (MQ=255) ggAATCGCTTATGTTCAAGCTGCACCGTACAGGATTGATCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGa < 1:212329/140‑1 (MQ=255) cGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGAc < 2:760475/140‑1 (MQ=255) cGCTTATGTTGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGAc < 2:767948/140‑1 (MQ=255) ttATGTTGAAGCTGAACCGTACGGGATTGTTAGTCTCGACTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTg < 2:235725/140‑1 (MQ=255) tGAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCtat < 2:670399/140‑1 (MQ=255) gAAGCTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCtata > 2:94344/1‑140 (MQ=255) cTGAACAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGcc < 2:500990/140‑1 (MQ=255) aaCAGTACAGGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGaa > 2:593885/1‑140 (MQ=255) | CTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAA > NC_000913/291807‑292072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |