| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,642,131 | A→G | 73.0% | Y57H (TAT→CAT) | essQ ← | predicted S lysis protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,642,131 | 0 | A | G | 73.0% | 61.2 / 21.3 | 37 | Y57H (TAT→CAT) | essQ | predicted S lysis protein |
| Reads supporting (aligned to +/- strand): ref base A (4/6); new base G (20/7); total (24/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.18e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGT > W3110S.gb/1642065‑1642198 | aTAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGa < 1:1005050/71‑1 (MQ=255) ttGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTg > 1:224403/1‑69 (MQ=255) ttGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGt > 1:1479127/1‑70 (MQ=255) tGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGAttt > 1:1787437/1‑59 (MQ=255) tcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGAtt > 1:797637/1‑48 (MQ=255) tcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCa > 1:2310937/1‑70 (MQ=255) tcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCa > 1:1710647/1‑70 (MQ=255) tcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCa > 1:1707325/1‑70 (MQ=255) ctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCaaa > 1:1716442/1‑70 (MQ=255) ctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCa < 1:1381742/51‑1 (MQ=255) ctctcGCAGACTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATATGTCAGCAGGCCaaa > 1:972571/1‑70 (MQ=17) tctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGa < 1:141076/52‑1 (MQ=255) ctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGa > 1:1172865/1‑51 (MQ=255) aCGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTAcc > 1:2251484/1‑69 (MQ=255) aCGCCTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTAcc > 1:2444813/1‑69 (MQ=255) ccTGTCTTCTTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGca < 1:589820/71‑1 (MQ=255) tGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAAcc < 1:1737032/56‑1 (MQ=255) gTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATAc < 1:1447695/37‑1 (MQ=255) gTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATAc < 1:386897/37‑1 (MQ=255) tcttctTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCacac < 1:2621239/70‑1 (MQ=255) tcttctTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACAcc > 1:65495/1‑71 (MQ=255) tcttctTTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACAcc > 1:1838972/1‑71 (MQ=255) ttcttTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGa > 1:2774901/1‑71 (MQ=255) cttTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCaaaaa > 1:1174221/1‑48 (MQ=255) cttTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCaaaaa > 1:47926/1‑48 (MQ=255) cttTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCa > 1:326271/1‑51 (MQ=255) cttTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCg < 1:1668458/68‑1 (MQ=255) cttTAATCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGa < 1:1535573/69‑1 (MQ=255) tAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGa > 1:1359093/1‑66 (MQ=255) tAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGa > 1:1187384/1‑66 (MQ=255) tAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTg < 1:1262284/69‑1 (MQ=255) aTCTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAAcc < 1:1419179/45‑1 (MQ=255) cTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAg > 1:337139/1‑45 (MQ=255) cTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAg > 1:40994/1‑45 (MQ=255) cTTGAAATGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAg > 1:467568/1‑45 (MQ=255) ttGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCa > 1:1425930/1‑66 (MQ=255) aaTGAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACtgt < 1:1342728/71‑1 (MQ=255) | ATAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGT > W3110S.gb/1642065‑1642198 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |