Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 624,924 | T→C | 11.5% | D395G (GAT→GGT) | tagE_2 ← | Poly(glycerol‑phosphate) alpha‑glucosyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 624,924 | 0 | T | C | 11.5% | 81.0 / 4.5 | 37 | D395G (GAT→GGT) | tagE_2 | Poly(glycerol‑phosphate) alpha‑glucosyltransferase |
Reads supporting (aligned to +/- strand): ref base T (22/10); new base C (2/2); total (25/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.88e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01 |
ATCTGAAGTAATCAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTATTTTTGGGACA > USA300TCH1516_ALE/624794‑625061 | aTCTGAAGTAATCAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTAAATTc > 2:224492/1‑141 (MQ=255) ttataaGTTCGTTAGGTACATATTTAATATCATAACTAATAAATGCAAAACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGTATCGCTATATTCT‑TGGTCTa < 1:268421/136‑1 (MQ=255) cAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGAtt > 1:94681/1‑141 (MQ=255) tCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCtt > 2:5268/1‑141 (MQ=255) aGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATAGCTATATTCT‑TGGTCTAGATTGCTTaaaaa > 2:241951/1‑141 (MQ=255) tatTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgt > 2:25637/1‑141 (MQ=255) tatTTAATATCATAACTAATAACTGGAAAACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgt > 1:378317/1‑141 (MQ=255) tAATATCATAACTAATAACTGGCTGACCAAGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCCCTATATTGT‑TGGTCTAGATTGCTTAAAAAGACGAgtaatt > 1:185638/1‑139 (MQ=255) tAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgtaagt > 1:107106/1‑141 (MQ=255) tAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCgtaagt > 1:165386/1‑141 (MQ=255) aaaaCTAATAACTCGAACACCATGTGTAAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATg < 2:179640/139‑1 (MQ=255) aTAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATg < 2:327668/141‑1 (MQ=255) cTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGAttt < 1:335714/141‑1 (MQ=255) cTGGAACACCATGTGCTAAAGATTAAATTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATc > 2:271389/1‑141 (MQ=255) acacCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCatat > 1:289683/1‑141 (MQ=255) cTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCa > 1:312919/1‑141 (MQ=255) ttCAAGTAGCGTTAAAGAGCAACCTACCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCAtt < 1:86199/141‑1 (MQ=255) tCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCAttt > 2:350191/1‑141 (MQ=255) tCAAGTAGCGATAAAGAGAAACCTGACATGCTACTTGTGATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAACTCAttt > 2:88412/1‑141 (MQ=255) cAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAAATCAAATAGGCATCACTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCACTCCATATTTTTGTATCAATTCATTTa > 1:233113/1‑141 (MQ=255) cAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATAAAAATCAAATAAACAACAC‑ATTTGGTGTTGTCTGGATTAGTAAAAAAACCACGTACGTAAACCTGGATTTCCAAGCCATAATTTTGTAGCACTTCAttga > 2:95754/1‑139 (MQ=255) aaGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTaa > 2:55313/1‑141 (MQ=255) aaGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCAGCGTAAGTAAACATGATTTTCCAATCCATATTGTTGTATCAATTCTTTTaa > 2:286535/1‑141 (MQ=255) aGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAAtttt > 1:348339/1‑141 (MQ=255) tacagcgaAACCTTCCATGTTACTTGTTATTAAACTCAAACAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTCATTTTTTACt < 1:56047/135‑1 (MQ=255) aaaGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCACTGCATTTGATTTGTTCCtt > 1:33916/1‑141 (MQ=255) aaGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATGAATTAATTTAATTTTTTACttt > 2:211782/1‑141 (MQ=255) aaGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGCCATCGCTATATTCT‑TGGTCTAGCTTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCACTTCATTTAATTTTTTACttt > 2:45583/1‑141 (MQ=255) tGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACcat < 1:288671/141‑1 (MQ=255) taAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGctc < 1:25637/141‑1 (MQ=255) taAACTCAAATAGGCACCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGctc < 2:378317/141‑1 (MQ=255) aCTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAt > 1:376015/1‑141 (MQ=255) cTCAAATATGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAtt < 1:1666/141‑1 (MQ=255) cTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAtt > 2:97685/1‑141 (MQ=255) aaaTAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTAttttt < 1:154915/141‑1 (MQ=255) tAGGCATCGCCATATTCG‑TGGTCTAGATTGCTCAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGATGCTCTATTTTTggg < 2:3522/141‑1 (MQ=255) gCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTATTTTTGGGACa > 1:376547/1‑141 (MQ=255) | ATCTGAAGTAATCAGTTCGTTAGGTCCATATTTAATATCATAACTAATAACTGGAACACCATGTGCTAAAGATTCAAGTAGCGCTAAAGAGAAACCTTCCATGTTACTTGTTATTAAACTCAAATAGGCATCGCTATATTCT‑TGGTCTAGATTGCTTAAAAAGCCGCGTAAGTAAACATGATTTTCCAATCCATATTTTTGTATCAATTCATTTAATTTTTTACTTTCAGAACCAAAACCATACATATGAAGCTCTATTTTTGGGACA > USA300TCH1516_ALE/624794‑625061 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |