| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | AP009048 | 3,824,324 | T→C | 100% | pseudogene (459/687 nt) | rph → | ECK3633:JW3618:b3643; ribonuclease PH |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | AP009048 | 3,824,324 | 0 | T | C | 98.3% | 206.5 / ‑4.2 | 60 | pseudogene (459/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (26/33); total (26/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC > AP009048/3824256‑3824388 | aCGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc > 1:2935435/1‑71 (MQ=255) cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAAccc < 1:2645234/70‑1 (MQ=255) cGGGTGCCTGCGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCg < 1:585790/71‑1 (MQ=255) cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg < 1:1415071/69‑1 (MQ=255) cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg < 1:2542651/69‑1 (MQ=255) cGTGGCGCTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAg < 1:174440/68‑1 (MQ=255) gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg < 1:3370479/70‑1 (MQ=255) gcgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt < 1:5098955/71‑1 (MQ=255) cgcTGGCAGATGCGCTACATAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGt > 1:3706329/1‑70 (MQ=255) cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg > 1:981643/1‑69 (MQ=255) cgcTGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa < 1:822376/71‑1 (MQ=255) tGGCAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGc > 1:4384472/1‑70 (MQ=255) caGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGcc < 1:256995/67‑1 (MQ=255) aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCTCAGtt < 1:4543053/71‑1 (MQ=255) aTGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt > 1:2023255/1‑71 (MQ=255) tGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt > 1:3721382/1‑71 (MQ=255) gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGAt < 1:651413/45‑1 (MQ=255) gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGa < 1:2948026/53‑1 (MQ=255) gcTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGtt < 1:1826880/67‑1 (MQ=255) cTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt < 1:4884105/71‑1 (MQ=255) aCAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCg > 1:477730/1‑71 (MQ=255) cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg > 1:2431471/1‑71 (MQ=255) cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg > 1:2369932/1‑71 (MQ=255) cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg > 1:4910736/1‑71 (MQ=255) cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg < 1:1248891/71‑1 (MQ=255) cAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg > 1:4058954/1‑71 (MQ=255) gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCggta < 1:2569743/71‑3 (MQ=255) aaGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAgg < 1:933474/44‑1 (MQ=255) tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa < 1:596347/47‑1 (MQ=255) tggtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTa < 1:3457649/47‑1 (MQ=255) gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgtaa < 1:1414674/70‑3 (MQ=255) gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta < 1:2351791/69‑2 (MQ=255) gtggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATtgta < 1:760980/69‑2 (MQ=255) ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTc < 1:3864843/58‑1 (MQ=255) ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg > 1:5218398/1‑71 (MQ=255) gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTg > 1:2803735/1‑55 (MQ=255) gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg > 1:2909608/1‑71 (MQ=255) gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg < 1:4858929/71‑1 (MQ=255) aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa > 1:2634112/1‑67 (MQ=255) aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt < 1:4445981/70‑2 (MQ=255) aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggt > 1:1718140/1‑69 (MQ=255) aaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACg < 1:3969534/68‑1 (MQ=255) aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg < 1:1766097/37‑1 (MQ=255) aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATg < 1:4326532/36‑1 (MQ=255) aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCg < 1:1356147/43‑1 (MQ=255) aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg > 1:1505201/1‑68 (MQ=255) aaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGa < 1:5118812/71‑1 (MQ=255) aCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa > 1:1355205/1‑64 (MQ=255) cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACggtg < 1:4834124/68‑3 (MQ=255) ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt < 1:3229632/48‑1 (MQ=255) gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGc > 1:5367070/1‑70 (MQ=255) gCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTCTCGCTATTGTAAACGGTGAAGc > 1:1179099/1‑70 (MQ=255) cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAAc < 1:1067627/61‑1 (MQ=255) cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTg < 1:664473/71‑1 (MQ=255) cTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGAGAAGCAATTTg > 1:2355962/1‑71 (MQ=255) gAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGc > 1:3037079/1‑70 (MQ=255) aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa > 1:3451413/1‑71 (MQ=255) aaaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGa > 1:180746/1‑71 (MQ=255) aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg > 1:4660370/1‑55 (MQ=255) aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt < 1:4163764/71‑1 (MQ=255) aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt > 1:1900461/1‑71 (MQ=255) aaccaaccCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATc < 1:4235171/71‑1 (MQ=255) ccAACCCGATGAAAGGGATGGTAGCCGCAGTTTCt > 1:906782/1‑35 (MQ=255) | ACGGGTGCCTGCGTGGCGCTGGTAGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATC > AP009048/3824256‑3824388 |
| 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: ‑ |