Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,827,438 | C→G | 54.5% | V227V (GTC→GTG) | 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,438 | 0 | C | G | 54.5% | 5.7 / 59.9 | 55 | V227V (GTC→GTG) | mutS | methyl‑directed mismatch repair protein |
Reads supporting (aligned to +/- strand): ref base C (11/14); new base G (8/22); total (19/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01 |
ATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCA > minE/1827373‑1827500 | aTCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgag < 1:1160925/69‑1 (MQ=255) tCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgaga > 1:1591806/1‑69 (MQ=255) tCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCgaga > 1:88971/1‑69 (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:2119034/1‑69 (MQ=255) tGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCg > 1:2961219/1‑69 (MQ=255) tGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCg > 1:800779/1‑69 (MQ=255) cgcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:2774249/1‑69 (MQ=255) cgcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg < 1:1160070/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:507020/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:3295733/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:2370520/69‑1 (MQ=255) gcgATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGt < 1:2645584/69‑1 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGCACGCACCACGCGGACTTtg > 1:527084/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:144784/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:2341533/2‑45 (MQ=255) agATCTGGTCGGTTTTGGTGTGGAGAACGCACCACGCGGACTTtg > 1:2366894/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) agATCTGGTCGGTTTTGGTGTGGAGAACCCACcacgcg > 1:1091939/2‑38 (MQ=16) cTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:1964172/1‑63 (MQ=255) cTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAg > 1:3353190/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:797977/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1303624/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:968293/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1439114/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:558665/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3231235/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3336366/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:3288434/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2670197/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2058237/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1627914/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2477186/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1992132/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1566501/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2886735/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2897094/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:1513068/68‑1 (MQ=255) ttttGGTGTGGAGAACGCACCACGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATa < 1:2984517/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) | ATCGACACCGCGCGCCAGCAGTTGAATCTGCAATTTGGGACCCGCGATCTGGTCGGTTTTGGCGTCGAGAACGCGCCGCGCGGACTTTGTGCTGCCGGTTGTCTGTTGCAGTATGCGAAAGATACCCA > minE/1827373‑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: ‑ |