| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,076,873 | C→T | 21.3% | H147Y (CAT→TAT) | ebgC → | cryptic beta‑D‑galactosidase, beta subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,076,873 | 0 | C | T | 21.3% | 92.8 / 19.4 | 47 | H147Y (CAT→TAT) | ebgC | cryptic beta‑D‑galactosidase, beta subunit |
| Reads supporting (aligned to +/- strand): ref base C (20/17); new base T (3/7); total (23/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.86e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACG > minE/2076806‑2076942 | tCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCat > 1:2615147/1‑70 (MQ=255) gTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataaca < 1:643099/71‑1 (MQ=255) tCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTataacaaata < 1:56183/70‑1 (MQ=255) tCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaataa > 1:1603395/1‑71 (MQ=255) tCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaata < 1:2101030/70‑1 (MQ=255) tCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaata < 1:21084/70‑1 (MQ=255) tCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaata < 1:2159013/70‑1 (MQ=255) aaCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATaacaac > 1:1072828/1‑67 (MQ=255) aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATaacaac > 1:3531915/1‑66 (MQ=255) aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATaacaac > 1:3527372/1‑66 (MQ=255) aCGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTAcgg > 1:2526387/1‑71 (MQ=255) gTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCa < 1:905060/71‑1 (MQ=255) cAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCaa > 1:2142926/1‑70 (MQ=255) aaaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaataa > 1:958947/1‑54 (MQ=255) aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAg < 1:1042732/71‑1 (MQ=255) aaaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAg > 1:1909046/1‑71 (MQ=255) aaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCat > 1:1922547/1‑43 (MQ=255) aaaaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCat > 1:1063123/1‑43 (MQ=255) aaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAg < 1:2580548/71‑1 (MQ=255) aaGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAg > 1:564180/1‑71 (MQ=255) aGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGt < 1:3198806/71‑1 (MQ=255) gTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGtt < 1:2639232/71‑1 (MQ=255) tGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaata < 1:457978/46‑1 (MQ=255) tGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCataacaaata < 1:2734507/46‑1 (MQ=255) tGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCaaaa > 1:352882/1‑64 (MQ=255) ctcAAAGTCACCATCGAAGATGGTTATTTCTataacaa < 1:3734688/38‑1 (MQ=255) ctcAAAGTCACCATCGAAGATGGTTATTTCCataacaaataac < 1:1696399/43‑1 (MQ=255) ctcAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgccgc > 1:872234/1‑71 (MQ=255) ctcAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgccgc > 1:2381174/1‑71 (MQ=255) ctcAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTgccg < 1:1362710/70‑1 (MQ=255) aaaGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTAcgg < 1:2713038/48‑1 (MQ=255) aaaGTCACCATCGAAGATGGTTATTTCCATAACAAAT‑ACAACTACGGCGGCAAAAGGAGTTTGCCGCCAc < 1:3898771/70‑1 (MQ=255) aaGTCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCAcc > 1:3880586/1‑71 (MQ=255) aaGTCACCATCGAAGATGGTTATTTCCATAACAAATACCACCTAcggcgg > 1:3882316/1‑50 (MQ=255) aaGCCACCATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCAcc > 1:3764629/1‑71 (MQ=255) ccATCGAAGATGGTTATTTCTATAACAAATAACAACTACGGCGGCAAAAGGa < 1:1012022/52‑1 (MQ=255) gAAGATGGTTATTTCCATAACAAATAACAACTAcgg > 1:3096960/1‑36 (MQ=255) gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACt > 1:184178/1‑71 (MQ=255) gAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACt > 1:2423348/1‑71 (MQ=255) aaGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACt > 1:320938/1‑70 (MQ=255) tGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCAtt < 1:3587979/70‑1 (MQ=255) ggTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCAttt > 1:2552671/1‑70 (MQ=255) ttATTTCCATAACAAATAACAACTACGGCGGCaaa < 1:3568200/35‑1 (MQ=255) ttATTTCCATAACAAATAACAACTACGGCGGCaaa < 1:3151313/35‑1 (MQ=255) tATTTCCATAACAAATAACAACTACGGCGGCAAAAgg < 1:1396458/37‑1 (MQ=255) tATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCg < 1:2779604/70‑1 (MQ=255) ttCTATAACAAATAACAACTACGGCGGCAAAAGGa < 1:3236067/35‑1 (MQ=255) ccATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACg > 1:2551321/1‑71 (MQ=255) | TCGGTTTATCTGCAATAACGCGGTCAAAAAAGTGGTTCTCAAAGTCACCATCGAAGATGGTTATTTCCATAACAAATAACAACTACGGCGGCAAAAGGAGTTTGCCGCCACCGCTACCCTACTCATTTTCGGAGACG > minE/2076806‑2076942 |
| 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: ‑ |