| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,348,085 | G→A | 31.7% | E89K (GAA→AAA) | npr → | phosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,348,085 | 0 | G | A | 31.7% | 54.3 / 30.8 | 41 | E89K (GAA→AAA) | npr | phosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr) |
| Reads supporting (aligned to +/- strand): ref base G (11/17); new base A (4/9); total (15/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.34e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAGCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCA > W3110S.gb/3348016‑3348150 | gaagCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTtgatg > 1:1104333/1‑70 (MQ=255) aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTa < 1:1600198/71‑1 (MQ=255) aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTa < 1:1631409/71‑1 (MQ=255) aCCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTa < 1:1030198/71‑1 (MQ=255) agaggaagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACttt < 1:2541099/71‑1 (MQ=255) gaggaagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTg < 1:134326/71‑1 (MQ=255) aggaagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGa < 1:2379140/71‑1 (MQ=255) ggaagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGa > 1:145043/1‑70 (MQ=255) aagCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACAt < 1:2551833/71‑1 (MQ=255) gCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:1555755/71‑1 (MQ=255) gCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATAc < 1:138870/71‑1 (MQ=255) gCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATAc < 1:2883717/71‑1 (MQ=255) ggCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTc > 1:2473582/1‑71 (MQ=255) gTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCa > 1:15112/1‑70 (MQ=255) gTTATCGCCCTCTTTAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCa < 1:1492348/70‑1 (MQ=255) ttATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTcc < 1:1166730/64‑1 (MQ=255) ttATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTcc < 1:586639/64‑1 (MQ=255) aTCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACt > 1:9763/1‑71 (MQ=255) aTCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCATCAAACt > 1:683101/1‑71 (MQ=255) tctTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:3607886/49‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:2737500/47‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:3151691/47‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:3695251/47‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATAc < 1:499911/47‑1 (MQ=255) tttAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACt < 1:819948/62‑1 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTcc > 1:97050/1‑62 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTcc > 1:3680703/1‑62 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa > 1:2197068/1‑71 (MQ=255) tAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa > 1:2832954/1‑71 (MQ=255) tAATTCTGGTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa > 1:923695/1‑71 (MQ=255) aaTTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACaa > 1:866155/1‑70 (MQ=255) gTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaaa < 1:3816866/70‑1 (MQ=255) gTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaa < 1:171144/69‑1 (MQ=255) gTTTTGATAAAGATTAATCTTCATCACTTTGACATACaa > 1:1653764/1‑39 (MQ=255) gTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaaa < 1:3816038/70‑1 (MQ=255) gTTTTGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaaa < 1:2908969/70‑1 (MQ=255) ttttGATGAAGATTAATCTTCATCACTTTGACATACa < 1:640360/37‑1 (MQ=255) ttttGATGAAGATTAATCTTCATCACTTTGACATACa < 1:3567145/37‑1 (MQ=255) ttttGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGaaa < 1:2779836/69‑1 (MQ=255) ttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAAt > 1:3607890/1‑69 (MQ=255) ttGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCa > 1:63072/1‑71 (MQ=255) ttGATAAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCa > 1:2001206/1‑71 (MQ=255) | GAAGCGACCGGTCCACAGGAAGAGGAAGCACTGGCCGCCGTTATCGCCCTCTTTAATTCTGGTTTTGATGAAGATTAATCTTCATCACTTTGACATACAACTCCCTTCAAACTCCCCCCGACAATAAGAAAATCA > W3110S.gb/3348016‑3348150 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |