Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 5,628 | A→G | 5.5% | intergenic (+98/+55) | yaaX → / ← yaaA | hypothetical protein/conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 5,628 | 0 | A | G | 5.5% | 209.1 / 2.6 | 73 | intergenic (+98/+55) | yaaX/yaaA | hypothetical protein/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (28/41); new base G (2/2); total (30/43) | |||||||||||
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.53e-01 |
TCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTA > minE/5560‑5697 | tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATg > 1:2321021/1‑71 (MQ=255) cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTTATTCTGCAATATATTGAATCTGCATGc > 1:772863/1‑71 (MQ=255) cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGc > 1:1961496/1‑71 (MQ=255) aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt > 1:3159980/1‑71 (MQ=255) gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt > 1:1413923/1‑71 (MQ=255) ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt < 1:1291823/70‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg < 1:2561009/71‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg < 1:376284/70‑1 (MQ=255) tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg > 1:3251649/1‑70 (MQ=255) gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg < 1:2569631/67‑1 (MQ=255) cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa > 1:3668316/1‑70 (MQ=255) cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt < 1:2723348/71‑1 (MQ=255) gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt > 1:3482613/1‑70 (MQ=255) gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa < 1:1389208/71‑1 (MQ=255) ggCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa < 1:841959/66‑1 (MQ=255) gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGc > 1:3199098/1‑71 (MQ=255) gcgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGt > 1:3160887/1‑71 (MQ=255) cgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt < 1:1265795/71‑1 (MQ=255) tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt < 1:2891438/69‑1 (MQ=255) cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa < 1:1661665/70‑1 (MQ=255) cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc > 1:2892450/1‑71 (MQ=255) tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc > 1:2386255/1‑71 (MQ=255) tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc > 1:2381545/1‑71 (MQ=255) aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg < 1:3764262/61‑1 (MQ=255) aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg < 1:3352438/61‑1 (MQ=255) aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg < 1:731706/69‑1 (MQ=255) aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg < 1:2396046/69‑1 (MQ=255) aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg < 1:1736853/69‑1 (MQ=255) aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc < 1:3514893/68‑1 (MQ=255) cAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg > 1:479535/1‑61 (MQ=255) ggCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCAt < 1:1797126/71‑1 (MQ=255) gCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc > 1:569159/1‑71 (MQ=255) cTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTa > 1:2274993/1‑42 (MQ=255) tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG > 1:2444461‑M1/1‑71 (MQ=255) cgctaatTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA < 1:556490‑M1/68‑1 (MQ=255) gTTAATTCTGCAATATATTGAATCTGCATGCtttt < 1:2291463/35‑1 (MQ=255) gTTAATTCTGCAATATATTGAATCTGCATGCtttt < 1:1661963/35‑1 (MQ=255) gTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa < 1:2316662/48‑1 (MQ=255) gTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa < 1:545881/48‑1 (MQ=255) ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt > 1:2268892/1‑71 (MQ=255) ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt < 1:2171837/71‑1 (MQ=255) aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt < 1:1609033/69‑1 (MQ=255) aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTg < 1:3008581/70‑1 (MQ=255) aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTg < 1:325016/70‑1 (MQ=255) ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACt < 1:1625082/71‑1 (MQ=255) tCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTg > 1:443700/1‑71 (MQ=255) tGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc > 1:816602/1‑70 (MQ=255) tGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc > 1:1338375/1‑70 (MQ=255) gCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt < 1:2632500/63‑1 (MQ=255) aaTATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGc < 1:3248275/40‑1 (MQ=255) aaTATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAAc < 1:2682855/71‑1 (MQ=255) atatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG < 1:1018048‑M1/56‑1 (MQ=255) tatatTGATTCTGCATTCTTTTGTAGGCAGTATAAGGCGtt < 1:3766705/41‑1 (MQ=255) tatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcac > 1:2025025‑M1/1‑58 (MQ=255) atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAT > 1:1835342‑M1/1‑57 (MQ=255) tatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa < 1:274357/64‑1 (MQ=255) tatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTa < 1:745164/70‑1 (MQ=255) tatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG < 1:2988205‑M1/53‑1 (MQ=255) ttGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactca < 1:3227803‑M1/70‑17 (MQ=255) ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg < 1:3882410/35‑1 (MQ=255) ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg < 1:1069140/35‑1 (MQ=255) ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt > 1:320899/1‑37 (MQ=255) ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc > 1:2448972/1‑45 (MQ=255) ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACg > 1:2275236/1‑71 (MQ=255) ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACg < 1:976455/71‑1 (MQ=255) gAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact > 1:1493863‑M1/1‑52 (MQ=255) gAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGAc < 1:727554/56‑1 (MQ=255) gAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgc < 1:710464/70‑1 (MQ=255) aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaaca < 1:1098384‑M1/70‑20 (MQ=255) atttgcgTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaacaat < 1:258357‑M1/68‑22 (MQ=255) tCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACt < 1:583103/61‑1 (MQ=255) tCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctgct > 1:727503/1‑71 (MQ=255) cTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTc > 1:2309488/1‑71 (MQ=255) tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa > 1:1755539/1‑55 (MQ=255) aTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGAc > 1:2661154/1‑48 (MQ=255) aTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCg > 1:2256088/1‑68 (MQ=255) aTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTa > 1:372398/1‑70 (MQ=255) | TCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTA > minE/5560‑5697 |
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: ‑ |