Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 19,348 | G→A | 15.7% | G335D (GGT→GAT) | gdpP → | Cyclic‑di‑AMP phosphodiesterase GdpP |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 19,348 | 0 | G | A | 15.7% | 73.3 / 7.4 | 32 | G335D (GGT→GAT) | gdpP | Cyclic‑di‑AMP phosphodiesterase GdpP |
Reads supporting (aligned to +/- strand): ref base G (17/10); new base A (2/3); total (19/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01 |
CGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTA > USA300TCH1516_ALE/19209‑19487 | cgcgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGt < 1:353845/141‑1 (MQ=255) cgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGa > 1:27552/1‑141 (MQ=255) aCCAAGTTTCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCa < 1:259373/141‑1 (MQ=255) cAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCAtt < 1:122413/141‑1 (MQ=255) aaGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCTCATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTa < 1:19340/141‑1 (MQ=255) aTTAAAAGTATTAATGGTAAAGTGCGTTTCTATGGCGGTAAGACGGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAAAGGGTGACAAGGTCATTATCATGGGa > 1:349302/1‑141 (MQ=255) ttAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTc > 1:256028/1‑141 (MQ=255) gTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt > 1:253205/1‑141 (MQ=255) tATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGCCATAAACGTCCTGACTTAGCTGCAATTGGt > 1:36769/1‑141 (MQ=255) aTGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTCGATGCAATTGGTg > 1:101167/1‑141 (MQ=255) tGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgc > 1:454463/1‑141 (MQ=255) gCGGTAAGACTGACCCGATGGATAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAACAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgcag < 1:22770/141‑1 (MQ=255) gCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACAGGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgcag > 1:243837/1‑141 (MQ=255) gtgTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTgcagc < 1:378178/139‑1 (MQ=255) tAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAAt < 1:443271/141‑1 (MQ=255) gaAACGGACTCGTGTAAGAGCACGTGTGATCTCACATGCATTAAAAGATATCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGc < 1:46980/141‑1 (MQ=255) gTACTGGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAatga < 1:215416/141‑1 (MQ=255) gTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAatga > 1:330912/1‑141 (MQ=255) gtgtAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGaata > 1:516353/1‑141 (MQ=255) gtAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGaataat > 1:526372/1‑141 (MQ=255) cACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGAGGAATAATTTAGAAg < 1:301889/141‑1 (MQ=255) cGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCa > 1:332367/1‑141 (MQ=255) tgATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACa > 1:53012/1‑141 (MQ=255) tgATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACa > 1:326571/1‑141 (MQ=255) cacaTGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTCCCATCGGACCTAAACGGCCATACTTAGATGCAAGTGGGGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTAt < 1:165749/141‑1 (MQ=255) ttAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAATCATACATCGTATTAAATgag > 1:17413/1‑141 (MQ=255) aaaaGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATGTAGAAGCATACATCGTATTAAATGAGAc > 1:306040/1‑141 (MQ=255) atatCCTTGCAGAGGATGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACa < 1:473943/141‑1 (MQ=255) tatCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACAt > 1:251153/1‑141 (MQ=255) ttGCAGAGGGTTACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATc < 1:265361/141‑1 (MQ=255) tGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGGCTGCCATTGATcc > 1:207476/1‑141 (MQ=255) cAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCaa > 1:327450/1‑141 (MQ=255) ggTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTa < 1:323029/141‑1 (MQ=255) | CGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATTGGTGCAGCAATCGGTGTGTCTAGATTTGCAATGATGAATAATTTAGAAGCATACATCGTATTAAATGAGACTGACATTGATCCAACATTA > USA300TCH1516_ALE/19209‑19487 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |