| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,933,819 | C→G | 12.5% | T208T (ACG→ACC) | lysS ← | lysine tRNA synthetase, constitutive |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,933,819 | 0 | C | G | 12.5% | 143.7 / 11.4 | 56 | T208T (ACG→ACC) | lysS | lysine tRNA synthetase, constitutive |
| Reads supporting (aligned to +/- strand): ref base C (27/22); new base G (3/4); total (30/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.93e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 | |||||||||||
TCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCA > minE/1933751‑1933886 | tCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCg < 1:214493/70‑1 (MQ=255) cAGAGCATTATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGttt < 1:3709489/71‑1 (MQ=255) aGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGttt > 1:1721954/1‑70 (MQ=255) aGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTc < 1:3296067/71‑1 (MQ=255) cgcgTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCaa > 1:2246622/1‑71 (MQ=255) attATGATGGGTAATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACtt > 1:3654041/2‑71 (MQ=255) tATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTc > 1:3135061/1‑70 (MQ=255) tggtggGCGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCAt < 1:3850904/71‑1 (MQ=255) tgatgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTcc > 1:2541611/4‑69 (MQ=255) gtgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATa > 1:2177956/1‑70 (MQ=255) atgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATa < 1:2810000/69‑1 (MQ=255) atgggtaATAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATa < 1:3696467/69‑1 (MQ=255) tgggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATa > 1:1411182/1‑69 (MQ=255) gggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAg > 1:1414807/1‑71 (MQ=255) ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAg < 1:2314223/70‑1 (MQ=255) ggtgATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGc > 1:3260724/1‑71 (MQ=255) gtgATAAACGGACGCGCAGCGGCTCCGCCAGGGATCACCTGCATCATCGGCGTTTCAACt > 1:1454886/1‑60 (MQ=255) aTAAACGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgc < 1:3420029/70‑1 (MQ=255) tAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCgg > 1:3834172/1‑71 (MQ=255) tAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCgg > 1:2073661/1‑71 (MQ=255) aaaCGGGCGAGCAGATGCCCCACCTGGAATTACCTGCATCATCGGGGTTTCTACTTCCATAAAGCcgcgcg > 1:269631/1‑69 (MQ=255) aaaCGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGt > 1:333916/1‑71 (MQ=255) aaCGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGtt < 1:3052754/71‑1 (MQ=255) gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc < 1:1877772/71‑1 (MQ=255) gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc < 1:1859456/71‑1 (MQ=255) gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc < 1:1816317/71‑1 (MQ=255) gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc < 1:1664655/71‑1 (MQ=255) gACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc < 1:3683957/71‑1 (MQ=255) aCGCGCAGCGGCACCGCCAGGGATCACGTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCa > 1:1597963/1‑71 (MQ=255) cgcgCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTcc > 1:2770204/1‑51 (MQ=255) gcgcAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCAt < 1:1095457/70‑1 (MQ=255) gCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTg < 1:550273/71‑1 (MQ=255) gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGc < 1:1325918/70‑1 (MQ=255) gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg < 1:2587624/71‑1 (MQ=255) gCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg < 1:565835/71‑1 (MQ=255) ccGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGa > 1:585602/1‑61 (MQ=255) cGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATa > 1:625992/1‑71 (MQ=255) gCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGaa > 1:1156056/1‑60 (MQ=255) gCCAGGGATCACCTGCATCATCGGCGTTTCAACATCCATAAAGCCGCGGTTCACCATGaa > 1:1745595/1‑60 (MQ=255) aGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATAc > 1:2608047/1‑68 (MQ=255) gATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCg < 1:154930/60‑1 (MQ=255) gATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagaga < 1:1872459/71‑1 (MQ=255) aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAcc > 1:3484387/1‑48 (MQ=255) aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGa > 1:854294/1‑60 (MQ=255) aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagagag < 1:3152919/71‑1 (MQ=255) aTCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagagag < 1:3835040/71‑1 (MQ=255) tCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCagagag > 1:1154825/1‑70 (MQ=255) tCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAgg < 1:2071003/71‑1 (MQ=255) aCCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCAc > 1:989205/1‑44 (MQ=255) aCCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGAt > 1:2978667/1‑71 (MQ=255) tGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTg > 1:3429732/1‑71 (MQ=255) catcatCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTgg > 1:2739514/1‑49 (MQ=255) catcatCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGc > 1:14197/1‑70 (MQ=255) tcGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAgcgc < 1:3439968/71‑1 (MQ=255) cGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCa > 1:2287857/1‑71 (MQ=255) ggCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCa < 1:2028368/70‑1 (MQ=255) | TCCAGCGCGTTATGGTGGGTGATAAACGGACGCGCAGCGGCACCGCCAGGGATCACCTGCATCATCGGCGTTTCAACTTCCATAAAGCCGCGGTTCACCATGAACTGGCGAATACCAGAGAGGATCTGCGAGCGCA > minE/1933751‑1933886 |
| 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: ‑ |