| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 973,796 | C→T | 38.6% | G10S (GGC→AGC) | ycbC ← | conserved inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 973,796 | 0 | C | T | 38.6% | 47.0 / 61.8 | 57 | G10S (GGC→AGC) | ycbC | conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (17/18); new base T (16/6); total (33/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01 | |||||||||||
AAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCG > W3110S.gb/973728‑973862 | aaaCCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTaa > 1:723961/1‑71 (MQ=255)aaaCCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTaa > 1:4477859/1‑71 (MQ=255) aaCCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAAt < 1:3623276/71‑1 (MQ=255) aCCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAAt > 1:2187776/1‑70 (MQ=255) agcGCCAGGCCAGCGCCTATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttt < 1:1177935/69‑1 (MQ=255) ccaggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAgt > 1:4664811/1‑71 (MQ=255) aggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAg > 1:2919817/1‑68 (MQ=255) aggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAgt > 1:4937163/1‑69 (MQ=255) aggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAgt > 1:2121259/1‑69 (MQ=255) aggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAgt > 1:2532866/1‑69 (MQ=255) gcgcCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAg < 1:1758614/71‑1 (MQ=255) cgcCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGc < 1:1960437/71‑1 (MQ=255) cgcCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGc < 1:2373122/71‑1 (MQ=255) ccAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCAt > 1:564581/1‑71 (MQ=255) cAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATc > 1:4784598/1‑71 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTc > 1:3723176/1‑51 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTc < 1:2706682/71‑1 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaaaa > 1:2994537/1‑60 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaaaa > 1:1809814/1‑60 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCttt > 1:2065372/1‑70 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCttt > 1:3397321/1‑70 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTc < 1:398630/71‑1 (MQ=255) aacaTCAGCGGTAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTgg < 1:1681623/68‑1 (MQ=255) aacaTCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTgg < 1:935647/68‑1 (MQ=255) caTCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGtt > 1:2361412/1‑68 (MQ=255) caTCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGtt > 1:821794/1‑68 (MQ=255) caTCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGtt > 1:4756286/1‑68 (MQ=255) tCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAg < 1:1439382/52‑1 (MQ=255) tCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAg < 1:4288796/52‑1 (MQ=255) cAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttt > 1:2847922/1‑71 (MQ=255) aGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCttttt > 1:3453444/1‑71 (MQ=255) aGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCttttt > 1:2869475/1‑71 (MQ=255) aGCGGAAGTGG‑AGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCttttt > 1:2391185/1‑70 (MQ=255) cGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaa < 1:3465597/44‑1 (MQ=255) ggAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCt > 1:666681/1‑57 (MQ=255) ggAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTaa < 1:4105884/71‑1 (MQ=255) gAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCa > 1:2049400/1‑48 (MQ=255) gAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCa > 1:4546087/1‑48 (MQ=255) aGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCa < 1:3749200/71‑1 (MQ=255) gTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCt < 1:2940758/49‑1 (MQ=255) tGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCata < 1:2127075/71‑1 (MQ=255) ggCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCata > 1:547921/1‑70 (MQ=255) gcagcaACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCatat < 1:4282297/70‑1 (MQ=255) gcagcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCAtata < 1:1000721/71‑1 (MQ=255) cagcaACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAAc > 1:87284/1‑64 (MQ=255) cagcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCa < 1:2313379/66‑1 (MQ=255) gcaACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAAc > 1:3652865/1‑62 (MQ=255) gcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGc < 1:4052222/70‑1 (MQ=255) gcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCa < 1:3718872/71‑1 (MQ=255) aCATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTc < 1:2211866/49‑1 (MQ=255) aCATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTc < 1:2305977/49‑1 (MQ=255) atatTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGa > 1:764464/1‑71 (MQ=255) atatTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGa > 1:2594616/1‑71 (MQ=255) atatTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAggg > 1:1963026/1‑70 (MQ=255) tatTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCt < 1:4459723/37‑1 (MQ=255) ttGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCg > 1:1927018/1‑71 (MQ=255) ttGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATc > 1:2657327/1‑70 (MQ=255) ccAATCACTTTTTTCAGTGTATAAAGCATCCTTTc < 1:1636186/35‑1 (MQ=255) | AAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCG > W3110S.gb/973728‑973862 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |