| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,533,727 | C→T | 17.8% | A45T (GCG→ACG) | yiiM ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,533,727 | 0 | C | T | 17.8% | 61.5 / 7.4 | 28 | A45T (GCG→ACG) | yiiM | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (17/6); new base T (3/2); total (20/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.06e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01 | |||||||||||
AGAGATAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTGCG > W3110S.gb/3533661‑3533797 | agagATAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCt > 1:1210205/1‑69 (MQ=255) gaTAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTc > 1:2212985/1‑70 (MQ=255) gaTAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCg > 1:3118208/1‑71 (MQ=255) tAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGTCTGCTCGTc < 1:940124/71‑1 (MQ=255) tAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGTCTGCTCGTc < 1:2907616/71‑1 (MQ=255) gCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCt > 1:133101/1‑71 (MQ=255) gCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCt > 1:2579437/1‑71 (MQ=255) gCTCACGAGGATAATGACACAGCGCTCTCTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCt > 1:3095518/1‑71 (MQ=255) aCGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCTTc > 1:2687930/1‑69 (MQ=255) gACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTc > 1:711068/1‑49 (MQ=255) gACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTcc > 1:3202323/1‑71 (MQ=255) gACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTcc > 1:3050446/1‑71 (MQ=255) gACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTc > 1:2277463/1‑70 (MQ=255) ccgTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc < 1:2427305/71‑1 (MQ=255) ccgTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc < 1:2838599/71‑1 (MQ=255) ccgTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc < 1:1363337/71‑1 (MQ=255) ccgTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc < 1:1999157/71‑1 (MQ=255) ccgTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc < 1:1869344/71‑1 (MQ=255) cgTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc < 1:1424913/70‑1 (MQ=255) aaCCTTTTTCTCCGTCTGCTCGTCACCTTCCAGccc > 1:2382637/1‑36 (MQ=255) aaCCTTTTTCTCCGTCTGCTCGTCACCTTCCAGccc > 1:1023054/1‑36 (MQ=255) aaCCTTTTTCTCCGTCTGCTCGTCACCTTCCAGccc > 1:168740/1‑36 (MQ=255) aaCCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGa > 1:1764559/1‑71 (MQ=255) ttttCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAAc > 1:466307/1‑61 (MQ=255) ttttCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGAttt > 1:2083022/1‑69 (MQ=255) ttttCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGAttt > 1:2076248/1‑69 (MQ=255) ttttCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGAttt > 1:625135/1‑69 (MQ=255) ttttCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGAttt > 1:1959149/1‑69 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCa > 1:2364376/1‑36 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:114797/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:53295/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:379755/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:3050279/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:2842782/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:2730398/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTgcg > 1:2508639/1‑71 (MQ=255) ccTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCAACAACCTGGATTTTAGCAAATgc > 1:2021793/1‑70 (MQ=255) | AGAGATAATGCTCACGAGGATAATGACACAGCGCTCTGTCTGGCCCGCCGTGAACCTTTTTCTCCGCCTGCTCGTCACCTTCCAGCCCCAGCTCCGTCAGCATCAACTCACCATCAACCTGGATTTTAGCAATTGCG > W3110S.gb/3533661‑3533797 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |