| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 2,529,057 | G→C | 39.0% | A545P (GCG→CCG) | ptsI → | PEP‑protein phosphotransferase of PTS system (enzyme I) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 2,529,057 | 0 | G | C | 39.0% | 6.0 / 18.0 | 18 | A545P (GCG→CCG) | ptsI | PEP‑protein phosphotransferase of PTS system (enzyme I) |
| Reads supporting (aligned to +/- strand): ref base G (5/6); new base C (2/5); total (7/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.37e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01 | |||||||||||
AAATGGACTGGCATGTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGAAAAAACAATCTGCTAATCCACGAGATGCGGCCCAATTTACTGCTTAGGAGAAGATCATGGGT > CP009273/2528913‑2529198 | aaaTGGACTGGCATGTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATCCGaa > 2:62307/1‑149 (MQ=255) ctgaccccgCATGTGTGGTGAGCTTGCTGGCGATGACCGTGCTACACCTCTGTTGCTGGGGATGGGTCCCGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGc < 1:90104/136‑1 (MQ=255) tGGACTGGCATTTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGATCGCTATTTCTATCTCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGgt < 2:263493/149‑1 (MQ=255) gCATGTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATCCGAAGGTGTTagca > 2:53425/1‑149 (MQ=255) tgtgGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCTCGCGCATTAAGAAGATTTTTTGTAACACGAACTTCGAAGATGCGAAGGTGTTagcagagca < 2:9174/149‑1 (MQ=255) gtgAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGttgcgagaga > 1:200977/1‑136 (MQ=255) gCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATCCGAAGGTGTTAGCAGAGCAGGc < 1:254965/144‑1 (MQ=255) tACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGt > 2:219295/1‑149 (MQ=255) tGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATCCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACg < 2:125070/149‑1 (MQ=255) tGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAAAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGAGGCGGGtt > 2:16155/1‑149 (MQ=255) ggTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATCCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACGCTGGTTAAc < 1:57424/144‑1 (MQ=255) tgggACGAATTCTCGATGAGCGCCATTTCTATGCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATCCGAAGGTGTTAGCAGAGCAGTCTCTTGCTCAACCGTCAATGGACGAGTTAATGACGCTGGTTAACAAg < 1:11029/142‑1 (MQ=255) ctctATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCATGCTCTTGCTTAACCGACAACGTATGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGaaaaaa < 2:7125/149‑1 (MQ=255) aTTTCTATCCCGCCCATTAAGAAGATTATCCTTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCAGGCTCTTTCTCAACCGACAACTGTCTAGTTAATTATGCTTGTTAACAATTTCATTGAAGAAAAAACAATCTGCTAATc < 2:209859/149‑1 (MQ=255) ttAAGAAGATTATCCGTAACACGCACTTCGAAGATGCGAAGATATTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGAAAAAACAACATACTGACCCACGAGATGCGGCCCa > 2:3274/1‑149 (MQ=255) aagaagATTCTCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGAAAAAACAATTTGCTTATCCACGTGATGCGGc < 1:40543/144‑1 (MQ=255) aagaagATTATCCGTAACACGAACTTCGAAGATCCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGAAAAAACAATCTGCTAATCCACGGGATGCGGc < 1:244702/144‑1 (MQ=255) cGAAGATGCGAAGGTGTTAGCAGAGCAGGATCTTGCTCAACCGACAACGAACGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGAAAAAACAATCTGCTAATCCACGAGATGCGGCCCAATTTACTGCTTAGGAGAAGATCATGGGt > 2:158946/1‑149 (MQ=255) | AAATGGACTGGCATGTGTGGTGAGCTTGCTGGCGATGAACGTGCTACACTTCTGTTGCTGGGGATGGGTCTGGACGAATTCTCTATGAGCGCCATTTCTATCCCGCGCATTAAGAAGATTATCCGTAACACGAACTTCGAAGATGCGAAGGTGTTAGCAGAGCAGGCTCTTGCTCAACCGACAACGGACGAGTTAATGACGCTGGTTAACAAGTTCATTGAAGAAAAAACAATCTGCTAATCCACGAGATGCGGCCCAATTTACTGCTTAGGAGAAGATCATGGGT > CP009273/2528913‑2529198 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |