Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 227,251 | A→G | 23.1% | K15K (AAA→AAG) | crl → | DNA‑binding transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 227,251 | 0 | A | G | 23.1% | 41.7 / 10.7 | 26 | K15K (AAA→AAG) | crl | DNA‑binding transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base A (11/9); new base G (3/3); total (14/12) | |||||||||||
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.12e-01 |
TGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT > minE/227185‑227319 | tGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTa > 1:2098637/1‑71 (MQ=255) aTCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCg > 1:902363/1‑71 (MQ=255) aTCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCg < 1:806920/71‑1 (MQ=255) cacaACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCa > 1:3232931/1‑71 (MQ=255) acaacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCa < 1:1061945/70‑1 (MQ=255) caacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCAc > 1:534566/1‑70 (MQ=255) gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCg > 1:1869500/1‑70 (MQ=255) gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCa < 1:1749346/61‑1 (MQ=255) aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGta < 1:2390814/71‑1 (MQ=255) tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt > 1:1508111/1‑71 (MQ=255) gACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGt > 1:1147237/1‑70 (MQ=255) aCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTg < 1:485637/70‑1 (MQ=255) ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGaa < 1:1462802/69‑1 (MQ=255) ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTaa < 1:1617735/71‑1 (MQ=255) cGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGACCGTATATTCGTGAAGGTAAg < 1:1492310/71‑1 (MQ=255) gACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGaaa < 1:3453113/70‑1 (MQ=255) cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCa > 1:1613836/1‑67 (MQ=255) gAAGAGCAGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTg > 1:2989335/1‑50 (MQ=255) gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATc > 1:1558451/1‑71 (MQ=255) aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATc > 1:2136965/1‑70 (MQ=255) agagCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGa > 1:1975703/1‑71 (MQ=255) cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCt > 1:863557/1‑71 (MQ=255) aGATTGATCAAAAAGTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCtt < 1:388054/71‑1 (MQ=255) gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt < 1:3254169/71‑1 (MQ=255) gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt < 1:1445585/71‑1 (MQ=255) gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCttt > 1:3751720/1‑71 (MQ=255) aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTg < 1:1551411/71‑1 (MQ=255) aaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTg > 1:1786757/1‑70 (MQ=255) aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCaa < 1:2821710/45‑1 (MQ=255) aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTc > 1:3615053/1‑71 (MQ=255) aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTc > 1:536249/1‑71 (MQ=255) aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTc > 1:1801330/1‑71 (MQ=255) aaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt > 1:3658840/1‑71 (MQ=255) aTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGAt < 1:1602857/46‑1 (MQ=255) | TGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCT > minE/227185‑227319 |
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: ‑ |