Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 744,026 | G→C | 38.6% | G226A (GGT→GCT) | nagZ → | beta N‑acetyl‑glucosaminidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 744,026 | 0 | G | C | 38.6% | 38.0 / 48.7 | 57 | G226A (GGT→GCT) | nagZ | beta N‑acetyl‑glucosaminidase |
Reads supporting (aligned to +/- strand): ref base G (29/6); new base C (9/13); total (38/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.61e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGA > minE/743958‑744074 | aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:1052595/69‑1 (MQ=255) aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:2038634/69‑1 (MQ=255) aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:581824/69‑1 (MQ=255) aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:2555192/69‑1 (MQ=255) aaaTAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:498184/69‑1 (MQ=255) aTAAACTCGACTCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:2657376/67‑1 (MQ=255) cTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCgg < 1:143267/62‑1 (MQ=255) cTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTccc < 1:710391/69‑1 (MQ=255) cTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTccc < 1:2687083/69‑1 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:1537655/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:2110800/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:600667/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:1962219/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:1937982/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:1002015/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:446154/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:2774724/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:2875372/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:3405107/1‑69 (MQ=255) aTGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGaaa > 1:758119/1‑69 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:871351/68‑1 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:2839318/68‑1 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:3286625/68‑1 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:1479140/68‑1 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:2300647/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:1319731/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:3357488/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:1759735/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:1516556/1‑69 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:656377/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:3399937/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:3027927/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:1421375/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:1487598/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:1621789/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:670702/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:637812/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:201991/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:203900/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:2341893/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:2083788/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTAATGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:62069/65‑1 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:57566/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:3311075/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:2604424/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:2465361/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:1040451/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTAATGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:1646673/1‑69 (MQ=255) ttGATGCGCGTCCGGCGAGCGGTTCTCACTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:2590086/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:3387991/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:306386/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:855217/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:1511340/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:3257963/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:2150274/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:313918/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:1110897/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:1638043/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:1751925/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:235777/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:2493275/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:2663203/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:3027576/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGGTTTGa > 1:3112835/1‑69 (MQ=255) | AAATAAACTCGACGCCATTATGCCTGCGCATGTGATCTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGA > minE/743958‑744074 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |