| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 732,825 | C→T | 7.9% | S1081S (TCC→TCT) | rhsC → | Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 732,825 | 0 | C | T | 7.9% | 108.5 / 3.0 | 51 | S1081S (TCC→TCT) | rhsC | Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor |
| Reads supporting (aligned to +/- strand): ref base C (29/17); new base T (2/2); total (32/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.29e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01 | |||||||||||
CGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAA > NC_000913/732699‑732973 | cGACAGACCCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAACCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAgtg < 2:255298/149‑1 (MQ=17) gACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCGCCAGCCCGGTGGCGAAGACGGTGGCAgtgt > 1:726094/1‑149 (MQ=12) gACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAgtga > 1:595293/1‑148 (MQ=17) ccGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTcc > 1:426630/1‑149 (MQ=32) ccGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCAACAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGAGCCGCAGCCTGGCGGATACCCTTCAGCAGGCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTcc > 2:566010/1‑149 (MQ=12) cGCATCCAGATGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTccc < 1:651602/149‑1 (MQ=17) ggACGATTTATCAGCCGGGGAGCTTCACGCCGCTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTg > 2:572307/2‑149 (MQ=17) gACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAATACGGTGGCAGTGTGGTGTTCCCGCCGGTGc > 2:454729/1‑149 (MQ=17) gACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGc > 1:98026/1‑149 (MQ=32) cAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGAt > 2:195821/1‑149 (MQ=32) cAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGGGTGGTGTTCCCGCCGGTGCTGGTGCGGAt > 1:197539/1‑149 (MQ=17) cAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTAGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGCTGTTCCCGCCAGTGCTGGTGCAGAt > 1:476422/1‑149 (MQ=255) aGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATg > 2:87840/1‑149 (MQ=32) gCCGGGGGGCTTCACGCCACGCAGCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGc < 2:60893/149‑1 (MQ=12) ggggAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGGGTTCCCCCCGGGGCTGGTGCAGATGCCCg > 2:380441/1‑149 (MQ=255) gggAGCTTCACGCCACTCATCAGGGTTGAAACCGCCGCCGGTGAGCTGGCGAAAGCGCAGCGGGGCGGCCTGGCGGATACGCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGGGTTCCTGCCGGTGCTGGTGCGGGTGCCCGa > 2:774399/1‑149 (MQ=11) gggAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCTGCCGGTGCTGGTGCAGATGCTCGa > 1:723188/1‑149 (MQ=17) cTTCACGCCACTCAGCCGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGc < 1:187357/149‑1 (MQ=17) tCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAACCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGTTGGTGCAGATGCTCGACCGGCTg < 1:412009/149‑1 (MQ=17) tCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTg < 1:236721/149‑1 (MQ=32) ccGCCACGCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGa > 2:485222/2‑149 (MQ=17) aCGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGa > 2:458819/1‑149 (MQ=32) gCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaaa > 2:725638/1‑149 (MQ=32) gCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGCTGCTGGTCCAGATGCACGACCGGCTGGaaa > 1:207987/1‑149 (MQ=12) tcTCATCAGGGTTACAACCGCCACCGTCGAGCTGGCGAAACCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGTCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTg < 2:412305/148‑1 (MQ=11) aCTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTg < 1:33485/149‑1 (MQ=32) aCTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTg < 1:285996/149‑1 (MQ=32) tcatcaGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGaa > 1:770637/1‑149 (MQ=32) caGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATcc > 1:430374/1‑149 (MQ=32) ttGACACCGCCCGCGGTCATGTGGCCAAAGCGCTGCGGCGCAGCCTGGCGGGTATCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt < 2:532943/149‑1 (MQ=11) ttGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt < 2:56418/149‑1 (MQ=32) ttGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt < 2:568673/149‑1 (MQ=32) ttGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt > 1:594674/1‑149 (MQ=32) ttGAAACCGCCACCGGGGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt < 2:369680/149‑1 (MQ=17) aCCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGCCCg > 1:107237/1‑149 (MQ=17) ccGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAATACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg < 1:454729/149‑1 (MQ=17) ccGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg < 1:172190/149‑1 (MQ=32) aGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCTGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGaa > 2:29130/1‑149 (MQ=17) aGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCTGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGaa > 2:28781/1‑149 (MQ=17) tGGCGAAAACGCAACGCCGCAGCCTGGCGGATACCCGTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGGGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgc > 1:535855/1‑149 (MQ=12) gCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAACgccg > 1:640317/1‑149 (MQ=17) gAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGATCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAgccgcc > 2:565081/1‑149 (MQ=17) aaaaCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg < 1:329632/149‑1 (MQ=32) aaaCGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGACGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCcgccgc > 2:87261/1‑149 (MQ=17) cGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGGAGATGCTCGACCGGCTGGAAAGGGGAATCGTGGCTGGCCGGGTGAGTGAGGGAAGCCGCGgcagg > 2:57728/1‑149 (MQ=11) cgcAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATc < 2:430374/149‑1 (MQ=32) gcAGCCTGTCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCg > 1:571563/1‑149 (MQ=17) gcAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCGGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCACTGGCTGGCATCg > 2:604102/1‑149 (MQ=17) gcAGCCTGGCGGATACCCTTCAGCAGTCAGATTGGTAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCg > 2:363716/1‑149 (MQ=11) cTGGCGGATACCCTTCAGCAGTCTGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg < 1:29130/149‑1 (MQ=17) cTGGCGGATACCCTTCAGCAGTCTGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg < 1:28781/149‑1 (MQ=17) cGGTGGCGAAGACAGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCaaa > 1:533371/1‑149 (MQ=25) cGGGGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCaaa < 2:402012/149‑1 (MQ=17) | CGACAGAACCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTTGAAACCGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATACCCTTCAGCAGTCCGGTGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAA > NC_000913/732699‑732973 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |