Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 744,046 | A→G | 37.3% | T233A (ACC→GCC) | nagZ → | beta N‑acetyl‑glucosaminidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 744,046 | 0 | A | G | 37.3% | 27.5 / 35.7 | 59 | T233A (ACC→GCC) | nagZ | beta N‑acetyl‑glucosaminidase |
Reads supporting (aligned to +/- strand): ref base A (26/11); new base G (9/13); total (35/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.23e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-02 |
CTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGATGGA > minE/743994‑744107 | cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:3286625/68‑1 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:871351/68‑1 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:1479140/68‑1 (MQ=255) cTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAgg < 1:2839318/68‑1 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:1516556/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:1759735/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:2300647/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:3357488/1‑69 (MQ=255) cAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACt > 1:1319731/1‑69 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:3399937/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:637812/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:2341893/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:2083788/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:656377/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:3027927/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:203900/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:201991/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:670702/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:1421375/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:1621789/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:1487598/65‑1 (MQ=255) gTGATGTTGATCCGCGTCCGGCGAGCGCTTCTCCCTAATGGCTGAAAGCCGTTTTGCGTCAGGaa < 1:62069/65‑1 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:3311075/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:57566/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:1040451/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:2465361/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:2604424/1‑69 (MQ=255) aTGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTAATGGCTGAAAACCGTTTTGCGTCAGGAACTGGGtt > 1:1646673/1‑69 (MQ=255) ttGATGCGCGTCCGGCGAGCGGTTCTCACTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:2590086/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:855217/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:3257963/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:306386/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:3387991/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:1511340/1‑69 (MQ=255) ttGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTg > 1:2150274/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:2493275/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:2663203/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:235777/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:3027576/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:1638043/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:313918/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:1751925/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGTTTTGa > 1:1110897/1‑69 (MQ=255) tGATCCGCGTCCGGCGAGCGCTTCTCCCTACTGGCTGAAAGCCGTTTTGCGTCAGGAATTGGGGTTTGa > 1:3112835/1‑69 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:992559/69‑1 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:2546332/69‑1 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:2471601/69‑1 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:2183539/69‑1 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:1972416/69‑1 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:1968764/69‑1 (MQ=255) ctcCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGAtt < 1:1444102/69‑1 (MQ=255) cTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTg > 1:2202536/1‑46 (MQ=255) cTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGa > 1:403161/1‑68 (MQ=255) cTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGa > 1:2193645/1‑68 (MQ=255) cTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGa > 1:206589/1‑68 (MQ=255) gCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGATGGa > 1:3230358/1‑69 (MQ=255) gCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGATGGa > 1:284141/1‑69 (MQ=255) gCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGATGGa > 1:2009447/1‑69 (MQ=255) gCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGATGGa > 1:698391/1‑69 (MQ=255) | CTACAGTGATGTTGATCCGCGTCCGGCGAGCGGTTCTCCCTACTGGCTGAAAACCGTTTTGCGTCAGGAACTGGGTTTTGACGGCGTGATTTTCTCTGACGATTTATCGATGGA > minE/743994‑744107 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |