Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 335,162 | G→A | 78.3% | G130G (GGC→GGT) | ybaK ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 335,162 | 0 | G | A | 78.3% | 116.7 / 32.0 | 60 | G130G (GGC→GGT) | ybaK | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (5/8); new base A (17/30); total (22/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 |
AAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGT > minE/335096‑335222 | aaaTTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGccaccg > 1:2195544/1‑70 (MQ=255) gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCTATATCCAGTCCGCGCTTGCCACCGGAAACAt < 1:2608908/70‑1 (MQ=255) gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa < 1:2166512/71‑1 (MQ=255) gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACAt < 1:1219997/70‑1 (MQ=255) gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACAt < 1:354985/70‑1 (MQ=255) cATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa > 1:4053717/1‑70 (MQ=255) cATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa > 1:735941/1‑70 (MQ=255) tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt < 1:2848182/71‑1 (MQ=255) tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt < 1:2924940/71‑1 (MQ=255) tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt < 1:2959134/71‑1 (MQ=255) tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt < 1:3297988/71‑1 (MQ=255) tCGAGGATCTTTGCCAGATCGCCTGACGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt < 1:1138807/71‑1 (MQ=255) aTCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGaa < 1:61388/57‑1 (MQ=255) tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt < 1:2532597/70‑1 (MQ=255) tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt < 1:3610426/70‑1 (MQ=255) tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt < 1:3178271/70‑1 (MQ=255) gCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAtt < 1:1378834/70‑1 (MQ=255) gCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAtt < 1:3481334/70‑1 (MQ=255) ccAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCt < 1:1803517/71‑1 (MQ=255) ccAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCt < 1:305597/71‑1 (MQ=255) cAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCtt < 1:166438/71‑1 (MQ=255) cAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCtt < 1:180308/71‑1 (MQ=255) gATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTtgt < 1:1450490/71‑1 (MQ=255) cGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:3282782/1‑56 (MQ=255) cGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:156035/1‑56 (MQ=255) gCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCgg > 1:3481094/1‑71 (MQ=255) gCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCgg > 1:1205851/1‑71 (MQ=255) gccgccAGTTCGATATCCATTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc < 1:1421727/59‑1 (MQ=255) gccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc < 1:442779/59‑1 (MQ=255) gccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc < 1:3427281/59‑1 (MQ=255) ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAg > 1:87097/1‑48 (MQ=255) ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAg > 1:3260895/1‑48 (MQ=255) ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAg > 1:843015/1‑48 (MQ=255) ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc < 1:3218345/58‑1 (MQ=255) aGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGAt < 1:462743/71‑1 (MQ=255) aGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGAt < 1:1674630/71‑1 (MQ=255) gTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTc < 1:2464489/52‑1 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa < 1:286851/36‑1 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa < 1:574340/36‑1 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:3797925/1‑43 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:1249907/1‑43 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:1789061/1‑43 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:405130/1‑43 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:1881157/1‑43 (MQ=255) ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg > 1:188440/1‑43 (MQ=255) atatCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCg < 1:3982755/70‑1 (MQ=255) atCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTg > 1:3326903/1‑71 (MQ=255) tCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGt < 1:1441922/35‑1 (MQ=255) tCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGt < 1:2511775/35‑1 (MQ=255) tCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCg < 1:896342/67‑1 (MQ=255) tCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCg < 1:1746902/67‑1 (MQ=255) tCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGt < 1:1312870/68‑1 (MQ=255) tCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt < 1:335603/37‑1 (MQ=255) tCCGCGCTTGCCACCGGAAAAATAAATAGTGGCAAAt < 1:3291729/37‑1 (MQ=38) ccGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt < 1:17040/36‑1 (MQ=255) ccGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc > 1:2394663/1‑38 (MQ=255) cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGc > 1:1420698/1‑65 (MQ=255) cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGt > 1:2693409/1‑71 (MQ=255) cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGt > 1:2344861/1‑71 (MQ=255) cTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGAc > 1:1972735/1‑65 (MQ=255) cTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGAc > 1:1454264/1‑65 (MQ=255) | AAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGT > minE/335096‑335222 |
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: ‑ |