| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,827,435 | C→T | 51.7% | G226G (GGC→GGT) | mutS → | methyl‑directed mismatch repair protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,827,435 | 0 | C | T | 51.7% | ‑3.9 / 51.1 | 58 | G226G (GGC→GGT) | mutS | methyl‑directed mismatch repair protein |
| Reads supporting (aligned to +/- strand): ref base C (14/14); new base T (8/22); total (22/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.04e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01 | |||||||||||
TTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCA > minE/1827367‑1827500 | tttGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGc > 1:3138469/1‑69 (MQ=255)tttGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGc > 1:1614118/1‑69 (MQ=255)tttGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGc > 1:1636380/1‑69 (MQ=255)tttGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGc > 1:1719810/1‑69 (MQ=255)tttGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGc > 1:1832648/1‑69 (MQ=255)tttGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGc > 1:2053020/1‑69 (MQ=255) tGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCg > 1:563802/1‑68 (MQ=255) tGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCg > 1:1088973/1‑68 (MQ=255) tGAAATCGACACCGCGCGCCAGCAGTTGAATCGGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCg > 1:2141631/1‑68 (MQ=255) aTCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgag < 1:1160925/69‑1 (MQ=255) tCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgaga > 1:88971/1‑69 (MQ=255) tCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgaga > 1:1591806/1‑69 (MQ=255) tCGACACCGCGCGCCAGCAGATGAATCTGCATTTTGGGACCCGCGATCTGGTCGGTTTTGGctt > 1:255794/1‑62 (MQ=255) ccAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCgg < 1:1290479/69‑1 (MQ=255) ccAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCgg < 1:583740/69‑1 (MQ=255) cagcagTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCgg < 1:1144777/68‑1 (MQ=255) gAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGaa > 1:997124/1‑48 (MQ=255) tGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGcc > 1:1294568/1‑68 (MQ=255) tGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCg > 1:800779/1‑69 (MQ=255) tGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCg > 1:2119034/1‑69 (MQ=255) tGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCg > 1:2961219/1‑69 (MQ=255) cgcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg < 1:1160070/69‑1 (MQ=255) cgcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:2774249/1‑69 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:507020/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:3295733/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:2645584/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:2370520/69‑1 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGCACGCACCACGCGGACTTtg > 1:527084/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:1004352/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:737654/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:2529267/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:2366894/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:2341533/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:144784/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACCCACcacgcg > 1:1091939/2‑38 (MQ=16) cTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:3353190/1‑63 (MQ=255) cTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:1964172/1‑63 (MQ=255) cTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:3048850/1‑63 (MQ=255) ttttGGTTTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1455769/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTTTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3276514/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3336366/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2058237/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1992132/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:558665/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1513068/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:797977/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1566501/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:968293/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1627914/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2477186/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3288434/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3231235/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1303624/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2984517/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1439114/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2897094/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2886735/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2670197/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCAAGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2380528/68‑1 (MQ=255) tttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:693709/67‑1 (MQ=255) ttGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTg < 1:1370694/43‑1 (MQ=255) ggCGTCGAGAACGCGCCGCGCGGACTTTGTTCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCa < 1:2820505/68‑1 (MQ=255) ggCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCa < 1:3411074/68‑1 (MQ=255) ggCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCa < 1:2611801/68‑1 (MQ=255) ggCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCa < 1:91638/68‑1 (MQ=255) | TTTGAAATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCA > minE/1827367‑1827500 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |