| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,304,455 | T→C | 13.4% | R389R (CGT→CGC) | sbcB → | exonuclease I |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,304,455 | 0 | T | C | 13.4% | 84.7 / 6.8 | 38 | R389R (CGT→CGC) | sbcB | exonuclease I |
| Reads supporting (aligned to +/- strand): ref base T (13/19); new base C (3/2); total (17/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.34e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01 | |||||||||||
AGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAAT > minE/1304385‑1304515 | aGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGt < 1:130336/71‑1 (MQ=255) ataACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAAttt < 1:2535908/70‑1 (MQ=255) taACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTa < 1:3556971/70‑1 (MQ=255) tttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACt < 1:2295151/69‑1 (MQ=255) tttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTg < 1:2285595/70‑1 (MQ=255) aGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGata < 1:1925776/70‑1 (MQ=255) tGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTTATTTAc > 1:1450618/1‑56 (MQ=255) aCGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTAc > 1:479186/1‑54 (MQ=255) cAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTTATTTACCGGCACTGGAAATCACttt > 1:1808952/1‑70 (MQ=255) cAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGCAATTTACCGGCACTGGATATCACtttt < 1:2288073/71‑1 (MQ=255) aGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTttg > 1:3456236/1‑71 (MQ=255) aTCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCa < 1:2937259/53‑1 (MQ=255) tGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGata < 1:3408889/56‑1 (MQ=255) tGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGata < 1:1292094/56‑1 (MQ=255) tGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGata < 1:3609106/56‑1 (MQ=255) gcaATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGc < 1:1903967/44‑1 (MQ=255) aTGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGTATATCACTTTTGTTGATAAACGGAt < 1:979969/71‑1 (MQ=255) tGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACtttt < 1:2352772/56‑1 (MQ=255) tGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACAGGCACTGGATATCACtttt < 1:2068725/56‑1 (MQ=255) gAAAATTGTGCTTGAAACCGAGCCGCGTAATTTACCGGCa > 1:3522289/1‑40 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGTTATTTACCGGCa > 1:3363314/1‑39 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGTAATTTACCg > 1:3404852/1‑36 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGTAATTTACCg > 1:1921467/1‑36 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCa > 1:2654384/1‑39 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGa < 1:1159956/71‑1 (MQ=255) aaaaTTGTGCTGGAAACCGAGCCGCGGAATTTACCGGCa > 1:1825201/1‑39 (MQ=255) ttGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGaaaa > 1:3366048/1‑70 (MQ=255) ttGTGCTGGAAACCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAg > 1:2625318/1‑71 (MQ=255) tgCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtg < 1:1793637/71‑1 (MQ=255) cTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgtt > 1:655698/1‑71 (MQ=255) tGGAAACCGAGCCGCGTAATTTTCCGGCACTGGATATCa > 1:3677562/1‑39 (MQ=255) gAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTgttgtt < 1:2454119/71‑1 (MQ=255) aaaCCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgtt > 1:3564544/1‑67 (MQ=255) aaaCCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTg > 1:1002749/1‑57 (MQ=255) aaaCCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTg > 1:1120558/1‑57 (MQ=255) aaaCCGAGCCGCGCAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTc < 1:3403952/71‑1 (MQ=255) ccGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAAt < 1:3933186/71‑1 (MQ=255) ccGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAt < 1:1552623/47‑1 (MQ=255) ccGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAt < 1:1522701/47‑1 (MQ=255) tAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAtt > 1:1119269/1‑43 (MQ=255) tAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGAtt > 1:230869/1‑43 (MQ=255) | AGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAAT > minE/1304385‑1304515 |
| 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: ‑ |