| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | W3110S.gb | 942,500 | Δ1 bp | coding (1120/2445 nt) | dmsA → | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 942,498 | 0 | A | . | 97.3% | 141.6 / ‑2.0 | 37 | coding (1118/2445 nt) | dmsA | dimethyl sulfoxide reductase, anaerobic, subunit A |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (17/19); total (18/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01 | |||||||||||
ACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCT > W3110S.gb/942434‑942556 | aCGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGgt < 1:281832/71‑1 (MQ=255)aCGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGgt < 1:694568/71‑1 (MQ=255)aCGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGgt < 1:3890285/71‑1 (MQ=255) aTGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCaa < 1:1160147/71‑1 (MQ=255) aTGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCa > 1:1857840/1‑70 (MQ=255) ggCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGgtag > 1:1503162/1‑62 (MQ=255) ggCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAAcc > 1:1090130/1‑71 (MQ=255) ggCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAAcc > 1:2897189/1‑71 (MQ=255) ggCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAAcc > 1:2243530/1‑71 (MQ=255) aaaTCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTa > 1:51020/1‑71 (MQ=255) aaaTCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTa > 1:2929880/1‑71 (MQ=255) aaTCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAAcc < 1:422746/62‑1 (MQ=255) tCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGgtag < 1:2328221/51‑1 (MQ=255) aCTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGaa > 1:1066519/1‑42 (MQ=255) aCTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGc < 1:3724306/71‑1 (MQ=255) gtgtTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg < 1:914920/71‑1 (MQ=255) gtgtTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg < 1:2806494/71‑1 (MQ=255) gtgtTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg < 1:2536473/71‑1 (MQ=255) gtgtTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg < 1:1418131/71‑1 (MQ=255) gtgtTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg < 1:234320/71‑1 (MQ=255) tgtTCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAggg < 1:1575620/71‑1 (MQ=255) ttCCGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGtt < 1:2652232/57‑1 (MQ=255) ccGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATgg > 1:2045760/1‑71 (MQ=255) ccGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATgg > 1:4216096/1‑71 (MQ=255) ccGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATgg > 1:3859483/1‑71 (MQ=255) ccGGCAGACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATg < 1:3748799/70‑1 (MQ=255) gACAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAgg < 1:1917655/60‑1 (MQ=255) cAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCa > 1:1008255/1‑52 (MQ=255) cAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCa > 1:4438658/1‑52 (MQ=255) cAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGa > 1:527729/1‑60 (MQ=255) cAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCg > 1:4067494/1‑69 (MQ=255) cAAAATCATCAAATTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCg > 1:1357974/1‑69 (MQ=255) aaaTTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATc > 1:2502046/1‑42 (MQ=255) aaaTTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCt < 1:2652329/71‑1 (MQ=255) aaaTTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCt < 1:2646821/71‑1 (MQ=255) aaaTTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGCGATGGGGCCCGCAGCGTCACg > 1:1900537/1‑69 (MQ=255) aaTTGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCt < 1:305183/70‑1 (MQ=255) tGGCTCGTG‑AATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGa > 1:24614/1‑46 (MQ=39) | ACGCCGGAATGGGCCTCGCAAATCACTGGTGTTCCGGCAGACAAAATCATCAAATTGGCTCGTGAAATCGGTAGTACCAAACCGGCGTTTATCAGCCAGGGATGGGGCCCGCAGCGTCACGCT > W3110S.gb/942434‑942556 |
| 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: ‑ |