| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 1,064,082 | G→T | 19.0% | A136E (GCA→GAA) | USA300HOU_1000 ← | LytR family transcriptional regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 1,064,082 | 0 | G | T | 19.0% | 39.1 / 4.4 | 21 | A136E (GCA→GAA) | USA300HOU_1000 | LytR family transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base G (12/5); new base T (2/2); total (14/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.74e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01 | |||||||||||
GTGTTAGGTTCATTTAAGTCGTATGGCACATCATAGTATATGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTA > CP000730/1063944‑1064210 | gtgtTAGGTTCATTTAAGTCGTATGGCACATCATAGTATATGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGATACATTCATGGTAGCTTCTACCGAATCCATTG‑Cg > 1:66951/1‑141 (MQ=255) tgtTAGGTTCATTTAAGTCGTATGGCACATCATAGTATATGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTT‑Cgg < 2:275472/141‑1 (MQ=255) tgtTAGGTTCATTTAAGTCGTATGGCACATCATAGTATATGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑Cgg > 1:135665/1‑141 (MQ=255) aGTCGTATGGCACATCATAGTATATGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCat > 1:244807/1‑141 (MQ=255) atAGTATATGCCACCACGCTCATTTACATCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCACCAGGTACAATCAAGGTAGCCTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGt < 1:167/141‑1 (MQ=255) atGCCACCAAGCTCATTTACAGCCTCGACAATTGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTa < 2:281282/141‑1 (MQ=255) tGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCAAGTGCATGGGTCACCTTAt > 2:32476/1‑141 (MQ=255) tCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTT‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCa > 1:246938/1‑141 (MQ=255) aTTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCCTGGGCAATCTTTTCGTAATAGCCAAc < 2:334377/141‑1 (MQ=255) gCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGa < 1:5163/141‑1 (MQ=255) gACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTgg > 2:179794/1‑141 (MQ=255) aCTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGa > 2:298752/1‑141 (MQ=255) ctcACGTAATAATCAACGGGTACATTCATGGTAGCTACTACCGAATAAACTGTCGTGGAGT‑GACGGCCATATGCATGTGCATGGGTAATCTTATCGTAGGAGCCTACTTTAGGATTGTGGCTGATTGTAGCACATGGCATa > 2:332747/1‑141 (MQ=255) tcACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTT‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATAc < 2:53979/141‑1 (MQ=255) ataatCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTACGCAt > 2:292755/1‑141 (MQ=255) taatCAACAGGTACCTTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCAtt < 2:165654/141‑1 (MQ=255) ggTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGa > 2:78904/1‑141 (MQ=255) ggTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGa > 1:20614/1‑141 (MQ=255) cATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGACTGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTCATTTGACGCTttga > 2:198424/1‑141 (MQ=255) tGGTAGCTTCTACCGAATCCATTT‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTttgatt > 1:202273/1‑141 (MQ=255) aCCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAagta > 1:238121/1‑141 (MQ=255) | GTGTTAGGTTCATTTAAGTCGTATGGCACATCATAGTATATGCCACCAAGCTCATTTACAGCCTCGACAAATGCTTTCATATTGACTCTCACATAATAATCAACAGGTACATTCATGGTAGCTTCTACCGAATCCATTG‑CGGCAATTGGACCACCATATGCATGTGCATGGGTAATCTTATCGTAATAGCCAACTTTAGGAATGTAGCTGATAGTATCACGTGGAATACTAAGCATTCTAATTTGATGTTTTGATTGATTAAAAGTA > CP000730/1063944‑1064210 |
| 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: ‑ |