Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 981,683 | Δ1 bp | 17.2% | coding (11/533 nt) | ydfQ ← | predicted lysozyme |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 981,683 | 0 | C | . | 17.2% | 80.2 / 10.5 | 29 | coding (11/533 nt) | ydfQ | predicted lysozyme |
Reads supporting (aligned to +/- strand): ref base C (10/14); new base . (3/2); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.32e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01 |
AAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATAC > minE/981614‑981742 | aaGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATc < 1:315486/70‑1 (MQ=255) aaGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCt < 1:775100/71‑1 (MQ=255) gTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAtttt > 1:3578152/1‑70 (MQ=255) gTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATtttt < 1:297082/71‑1 (MQ=255) aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTtgt > 1:1955803/1‑71 (MQ=255) cTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTtgtg > 1:1202278/1‑71 (MQ=255) aggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTtgt < 1:129069/67‑1 (MQ=255) aggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTtgt < 1:2310866/67‑1 (MQ=255) aggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGtt < 1:3671268/69‑1 (MQ=255) gagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCa > 1:2096577/1‑70 (MQ=255) gagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCata < 1:2951577/71‑3 (MQ=255) gCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAccc < 1:3781851/71‑1 (MQ=255) gCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAccc < 1:1287965/71‑1 (MQ=255) gCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCa > 1:194978/1‑67 (MQ=255) aTGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCtttt < 1:2766284/57‑1 (MQ=255) tGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccg > 1:1015000/1‑71 (MQ=255) gcgcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGtt > 1:1477326/1‑59 (MQ=255) ccAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTc < 1:2929815/71‑1 (MQ=255) cAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCg < 1:150026/71‑1 (MQ=255) gCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTc > 1:883700/1‑71 (MQ=255) ccAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTtgtg < 1:3629350/49‑1 (MQ=255) tCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATtttt < 1:2030248/41‑1 (MQ=255) gcCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATtca > 1:3135063/1‑71 (MQ=255) aGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgccc > 1:2773503/1‑48 (MQ=255) cAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACtt > 1:118382/1‑71 (MQ=255) aGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACttt < 1:548852/71‑1 (MQ=255) aGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCa < 1:294943/71‑1 (MQ=255) aGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCa > 1:2383498/1‑71 (MQ=255) ccGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATAc < 1:2453139/71‑1 (MQ=255) tATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTAc < 1:2573675/39‑1 (MQ=255) tATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATtcatc > 1:2946711/1‑49 (MQ=255) tCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGa > 1:2026919/1‑39 (MQ=255) gATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACt < 1:3058801/51‑1 (MQ=255) | AAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCACCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATAC > minE/981614‑981742 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |