| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,824,533:1 | +G | 100% | pseudogene (666/687 nt) | rph → | ECK3633:JW3618:b3643; ribonuclease PH |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,824,533 | 1 | . | G | 92.0% | 73.0 / 0.0 | 25 | pseudogene (666/687 nt) | rph | ECK3633:JW3618:b3643; ribonuclease PH |
| Reads supporting (aligned to +/- strand): ref base . (2/0); new base G (12/11); total (14/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.87e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
AGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGC > W3110S.gb/3824464‑3824593 | aGTTCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:2574685/71‑1 (MQ=255)aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:3077997/71‑1 (MQ=255)aGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:3073790/71‑1 (MQ=255) gTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑gg > 1:3015130/1‑71 (MQ=255) gACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGGTgaatc > 1:906219/1‑70 (MQ=255) cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcgaatc > 1:2589909/1‑70 (MQ=255) cAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGACCGAGGGG‑gaatcgaatc < 1:1088597/70‑1 (MQ=255) aaGGCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCa > 1:1135389/1‑69 (MQ=255) ggCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTg > 1:2955416/1‑70 (MQ=255) gCGAGCCGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑ggg‑gga < 1:603375/56‑3 (MQ=255) ccGTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGAc < 1:3128187/71‑1 (MQ=255) gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTg < 1:3146124/62‑1 (MQ=255) gTTCACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGa > 1:2925693/1‑68 (MQ=255) cACCCATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑GGG‑GAATCGAATCCAt > 1:1964553/1‑56 (MQ=255) ccATGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcg < 1:2815257/47‑1 (MQ=255) tGAAGAGCTACTCACCTTGTTGGCGCTGGCCCGA‑g < 1:3672158/35‑1 (MQ=37) gAAGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑gaatcgaat > 1:704675/1‑46 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa > 1:1070846/1‑64 (MQ=255) aaGAGCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa > 1:1095182/1‑64 (MQ=255) agagCTACTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGa > 1:1263837/1‑56 (MQ=255) cTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGc < 1:131486/47‑1 (MQ=255) cTCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCa < 1:2816232/53‑1 (MQ=255) tCACCTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGa > 1:168848/1‑54 (MQ=255) caccTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACt > 1:2908440/4‑71 (MQ=255) accTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa > 1:1220052/3‑68 (MQ=255) accTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCaaa > 1:1913107/3‑68 (MQ=255) cTTGTTGGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGAc < 1:3096261/45‑1 (MQ=255) ttgttgGCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGaa < 1:2587878/50‑1 (MQ=255) ggCGCTGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTa < 1:3049618/70‑1 (MQ=255) gcgctGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAg < 1:2726289/68‑1 (MQ=255) gctGGCCCGAGGGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGc < 1:3541063/70‑1 (MQ=255) | AGTGCAGGGGACGGCAGAAGGCGAGCCGTTCACCCATGAAGAGCTACTCATCTTGTTGGCTCTGGCCCGA‑GGG‑GAATCGAATCCATTGTAGCGACGCAGAAGGCGGCGCTGGCAAACTGATTTTTAAGGC > W3110S.gb/3824464‑3824593 |
| 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: ‑ |