Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,461,154 | T→G | 100% | pseudogene (801/4605 nt) pseudogene (801/4104 nt) |
yfaS ← yfaS ← |
ECK2220:JW2222+JW2221:b4500; hypothetical protein ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,461,154 | 0 | T | G | 83.4% | 73.0 / 6.0 | 36 | pseudogene (801/4605 nt) pseudogene (801/4104 nt) | yfaS yfaS | ECK2220:JW2222+JW2221:b4500; hypothetical protein ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment |
Reads supporting (aligned to +/- strand): ref base T (6/0); new base G (21/9); total (27/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.03e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CCACAAGATCTCTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTC > minE/1461088‑1461208 | ccACAAGATCTCTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCttt > 1:1433959/1‑69 (MQ=255) ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca > 1:1103645/1‑68 (MQ=255) ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca > 1:2941412/1‑68 (MQ=255) ttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca > 1:913356/1‑68 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:1682500/69‑1 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:3039060/69‑1 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:3168823/69‑1 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:449701/69‑1 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:1540428/69‑1 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:151370/69‑1 (MQ=255) ggTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:746087/69‑1 (MQ=255) gTCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:1374790/68‑1 (MQ=255) tCCACACCAGAAGCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGAAACAAACACCACcgt < 1:2887718/67‑1 (MQ=255) cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc > 1:2941958/1‑69 (MQ=255) cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc > 1:1137598/1‑69 (MQ=255) cacaCCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCgtcgtc > 1:3074596/1‑69 (MQ=255) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:714551/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:669198/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:985412/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:479860/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:3367665/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:3166815/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:3158919/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1098081/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:2831582/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:2710983/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:210695/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1884499/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1841359/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1782258/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1567560/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1480407/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1301618/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1291267/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:1213874/1‑42 (MQ=37) gCTCGTTACCTGACACTTTGCTAAGCGCCACTGTATCGGaa > 1:3045739/1‑41 (MQ=37) gCTCGTTACCAGACACTTTGCTAAGCGCCACTGTATCGGaaa > 1:486377/1‑42 (MQ=25) | CCACAAGATCTCTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTC > minE/1461088‑1461208 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |