| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,297,298 | T→C | 18.3% | E258G (GAA→GGA) | katG ← | catalase/hydroperoxidase HPI(I) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,297,298 | 0 | T | C | 18.3% | 101.1 / 17.0 | 49 | E258G (GAA→GGA) | katG | catalase/hydroperoxidase HPI(I) |
| Reads supporting (aligned to +/- strand): ref base T (19/21); new base C (4/5); total (23/26) | |||||||||||
| 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.97e-01 | |||||||||||
CATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTCGCCGC > minE/2297230‑2297367 | cATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttct > 1:1013898/1‑71 (MQ=255)cATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttct > 1:3903140/1‑71 (MQ=255)cATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttc < 1:866875/70‑1 (MQ=255) tttGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttcttc > 1:1053312/1‑71 (MQ=255) tttGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttctt > 1:1915477/1‑70 (MQ=255) tGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTtcg < 1:1720557/70‑1 (MQ=255) aTGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgtcg > 1:1134820/1‑71 (MQ=255) aTGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgtcg > 1:169687/1‑71 (MQ=255) cGGACCGGCACCCTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCa > 1:3505490/1‑71 (MQ=255) ggACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCa < 1:2189018/70‑1 (MQ=255) aCCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttctt > 1:1112215/1‑59 (MQ=255) cGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCAt < 1:1547899/67‑1 (MQ=255) ggCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTc < 1:1736932/64‑1 (MQ=255) gCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCAt > 1:1968086/1‑71 (MQ=255) cACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATg < 1:3505656/71‑1 (MQ=255) aCCGTGGGTTTTACCCAGCGTATGACCCCCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGt > 1:3633673/1‑71 (MQ=255) gTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgtcg > 1:2618579/1‑55 (MQ=255) ggTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCg > 1:1093523/1‑70 (MQ=255) ggTTTTACCCAGCGTATGACCACCCGAAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGa < 1:3720498/71‑1 (MQ=255) ttttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCa < 1:1054863/60‑1 (MQ=255) tttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAgg > 1:643425/1‑71 (MQ=255) tttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAgg > 1:813026/1‑71 (MQ=255) ttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGt > 1:333733/1‑71 (MQ=255) ttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGt > 1:53705/1‑71 (MQ=255) tACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGt < 1:1940261/70‑1 (MQ=255) aCCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcg > 1:2235251/1‑71 (MQ=255) ccAGCGTATGACCACCGGCAATCTGCGCCACGGTTTCTTCGTCGtt < 1:1427249/46‑1 (MQ=39) ccAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgcg > 1:3585521/1‑71 (MQ=255) cAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgtc > 1:1961669/1‑42 (MQ=255) aGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgcgcg > 1:2035605/1‑71 (MQ=255) cGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa < 1:127650/71‑1 (MQ=255) cGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa < 1:3510414/71‑1 (MQ=255) gTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa > 1:3161269/1‑70 (MQ=255) tATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGAt < 1:3721668/70‑1 (MQ=255) tATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGAt < 1:1271502/70‑1 (MQ=255) aTGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGccc > 1:165503/1‑46 (MQ=255) accaccCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAgctg < 1:869108/71‑1 (MQ=255) ccaccCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgcgc > 1:2427975/1‑61 (MQ=255) caccCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGtt < 1:2276148/46‑1 (MQ=255) accCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGctgctg > 1:1957630/1‑71 (MQ=255) ccGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGcc > 1:1225455/1‑71 (MQ=255) gCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCg > 1:3930299/1‑70 (MQ=255) tCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCg < 1:3571705/66‑1 (MQ=255) tCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCg < 1:3567164/66‑1 (MQ=255) tCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGaa < 1:457556/71‑1 (MQ=255) tCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGt < 1:1412227/36‑1 (MQ=255) cAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGt > 1:3244956/1‑45 (MQ=255) aGCGGCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCTCGCGGATAGCTGCTGCCGCAGAAAg < 1:1070678/71‑1 (MQ=255) aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGaa > 1:3079101/1‑41 (MQ=255) aGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGaaa < 1:3640486/70‑1 (MQ=255) gcCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGTAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGc < 1:1625175/69‑1 (MQ=255) ccACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGa < 1:2792836/64‑1 (MQ=255) ccACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGt < 1:1882748/71‑1 (MQ=255) ccACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGt < 1:567458/71‑1 (MQ=255) cACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGtt < 1:583320/71‑1 (MQ=255) tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTcgccg < 1:925577/71‑1 (MQ=255) tcttcGTCGTTCATGCCCATGTTGCCTAAGGTCGTGCGGATAGCTGCTGCCGCAGAAAGCGGTTcgccgc < 1:1557675/70‑1 (MQ=255) | CATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTCGCCGC > minE/2297230‑2297367 |
| 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: ‑ |