Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,235,757 | G→T | 16.0% | intergenic (+350/‑149) | pgi → / → yjbE | glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,235,757 | 0 | G | T | 16.0% | 72.8 / 6.5 | 25 | intergenic (+350/‑149) | pgi/yjbE | glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein |
Reads supporting (aligned to +/- strand): ref base G (11/10); new base T (2/2); total (13/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01 |
TAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAAT > NC_000913/4235654‑4235901 | tAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTtat > 1:913331/1‑112 (MQ=255) tAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTtat < 2:913331/112‑1 (MQ=255) gCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGa < 4:442236/79‑1 (MQ=255) gCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGa > 3:442236/1‑79 (MQ=255) gCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTgtgt > 3:1393030/1‑149 (MQ=255) aGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGta < 4:826437/148‑1 (MQ=255) aGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGta > 3:826437/1‑148 (MQ=255) cTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAAATcccc > 1:1212884/1‑108 (MQ=255) cTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAAATcccc < 2:1212884/108‑1 (MQ=255) ttttGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAg < 4:532738/130‑1 (MQ=255) ttttGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAg > 3:532738/1‑130 (MQ=255) ttATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTtgtg > 3:1402955/1‑116 (MQ=255) ttATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTtgtg < 4:1402955/116‑1 (MQ=255) tatcaaaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt < 3:153706/123‑1 (MQ=255) tatcaaaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt > 4:153706/5‑127 (MQ=255) tttAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTaa > 4:1137269/1‑64 (MQ=255) tttAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTaa < 3:1137269/64‑1 (MQ=255) tttAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAAt > 3:425409/1‑149 (MQ=255) aaaTATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt < 3:597063/123‑1 (MQ=255) aaaTATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttttt > 4:597063/1‑123 (MQ=255) aaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAAtctc > 3:456339/1‑133 (MQ=255) aaTATAAGTTTCAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAAtctc < 4:456339/133‑1 (MQ=255) aGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAAt < 4:1393030/149‑1 (MQ=255) gTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCa < 4:1134603/74‑1 (MQ=255) gTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCa > 3:1134603/1‑74 (MQ=255) | TAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCAATGAAGGGAAAT > NC_000913/4235654‑4235901 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |