Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,025,079 | C→A | 11.8% | Y12* (TAC→TAA) | wbpZ → | glycosyl transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,025,079 | 0 | C | A | 11.8% | 83.9 / 2.9 | 34 | Y12* (TAC→TAA) | wbpZ | glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base C (13/17); new base A (2/2); total (15/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01 |
GTCTGCTGCTGTGCAGCGGGGGCTGCTGCATGCCGGTGGCTTTTCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAA > NC_002947/2024945‑2025228 | gTCTGCTGCCGTGCAGCGGGGGCTGCTGCATGCCGGTGGCTTTTCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGt < 1:299539/150‑1 (MQ=255) tgctgctgTGCAGCGGGGGCTGCTGCATGCCGGTGGCTTTTCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGc > 1:68594/1‑150 (MQ=255) gtgCAGCGGGGGCTGCTGCATGCCGGTGGCTTTTCTTGGGAGCGTTGCGCTCAAGGTATAGTCGAACTATGCAAGCTGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAAc < 1:277790/150‑1 (MQ=255) gtgCAGCGGGGGCTGCTGCATGCCGGTGGCTTTTCTTGCGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAAc < 2:695572/150‑1 (MQ=255) ggggCTGCTGCATGCCGGTGGCTTTTCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCAtt > 1:434628/1‑150 (MQ=255) gTGGCTTTTCTTGGGCGCGTTGCGCTCAGGGTACAGGCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCAGTTTTCAGATCGCCGAAgg < 2:160125/150‑1 (MQ=255) gTGGCTTTTCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAgg > 1:561172/1‑150 (MQ=255) ttttCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCAccc > 1:272486/1‑150 (MQ=255) ttGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGaa > 1:256752/1‑149 (MQ=255) gAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACg > 1:148170/1‑150 (MQ=255) aGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTg < 2:219781/101‑1 (MQ=255) aGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTg > 1:219781/1‑101 (MQ=255) ttGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTg < 2:434628/150‑1 (MQ=255) aGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCaga < 2:580657/150‑1 (MQ=255) tACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTAACCg > 2:184442/1‑72 (MQ=255) tACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTAACCg > 2:184402/1‑72 (MQ=255) tACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTAACCg < 1:184442/72‑1 (MQ=255) tACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTAACCg < 1:184402/72‑1 (MQ=255) aCAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGc < 2:561172/150‑1 (MQ=255) aCAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTACATGGGTGGCATCGAACAGTTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGc < 1:531189/150‑1 (MQ=255) gTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCtat > 1:311174/1‑150 (MQ=255) aaGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTg > 2:741212/1‑150 (MQ=255) ggAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGGCTCCATGGGt > 1:409515/1‑57 (MQ=255) ggAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGGCTCCATGGGt < 2:409515/57‑1 (MQ=255) cGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCa < 2:232470/150‑1 (MQ=255) gatCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCtata < 1:171933/120‑1 (MQ=255) gatCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCtata > 2:171933/1‑120 (MQ=255) gCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTg > 2:460278/1‑103 (MQ=255) gCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTg < 1:460278/103‑1 (MQ=255) tttttttAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCa < 1:741212/150‑1 (MQ=255) ttttAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGcc < 1:644777/110‑1 (MQ=255) ttttAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGcc > 2:644777/1‑110 (MQ=255) ttAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGa < 2:311174/150‑1 (MQ=255) ttAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGa < 2:68594/150‑1 (MQ=255) cccGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTa < 2:600477/124‑1 (MQ=255) cccGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTa > 1:600477/1‑124 (MQ=255) cccGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCaa > 2:416680/1‑150 (MQ=255) | GTCTGCTGCTGTGCAGCGGGGGCTGCTGCATGCCGGTGGCTTTTCTTGGGAGCGTTGCGCTCAGGGTACAGTCGAGGTTTACAAGTCGGTCAAGGAACGATGATCAAGGTTCTGCATTTTTTTAAAACGTATTACCCGGACTCCATGGGTGGCATCGAACAGGTCATTTTTCAGATCGCCGAAGGCACCCGTGAGCACGGCGTGCAATCAGAGGTGCTATACCTCAGCCCGCGTGGGGCTGCGCGCAATGAACCGGTAGGCAGCCACCTGACTCACCGCTCCAA > NC_002947/2024945‑2025228 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |