| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 600,571 | C→A | 32.4% | G153G (GGC→GGA) | dmsB → | dimethyl sulfoxide reductase, anaerobic, subunit B |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 600,571 | 0 | C | A | 32.4% | 60.3 / 24.9 | 49 | G153G (GGC→GGA) | dmsB | dimethyl sulfoxide reductase, anaerobic, subunit B |
| Reads supporting (aligned to +/- strand): ref base C (16/17); new base A (9/7); total (25/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.62e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCT > minE/600504‑600634 | aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:234379/1‑68 (MQ=255)aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:811681/1‑68 (MQ=255)aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:1443125/1‑68 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:3182408/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:1454588/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:3057055/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:1825213/69‑1 (MQ=255) gtTGCGGAAGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:2990846/69‑1 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:372127/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:3142960/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:3034058/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:1918571/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:1787895/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:148853/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:817901/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:937652/1‑68 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:638323/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:1677454/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:1654183/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2563538/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2137095/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2944792/1‑69 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:3245806/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2517649/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2293559/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2056622/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2052734/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:660962/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:1201173/69‑1 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGCTTTCGGACCTATCGACGAGCTGCGTAAAAAACa > 1:2233302/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:3133524/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:3232381/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:976327/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:2829044/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:1206470/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:1741115/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:2142610/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:2440226/1‑68 (MQ=255) cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa < 1:1151059/62‑1 (MQ=255) cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa < 1:2754409/62‑1 (MQ=255) ttGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:3000979/59‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2467035/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2096578/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1734795/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1324400/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2793834/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1202253/52‑1 (MQ=255) tacGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1405285/42‑1 (MQ=38) ttCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGgccgttgcgccgttg > 1:2016841/1‑60 (MQ=255) ttCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCt < 1:2515862/69‑1 (MQ=255) ttCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCt < 1:90452/69‑1 (MQ=255) ttCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCt > 1:287806/1‑69 (MQ=255) | ACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCT > minE/600504‑600634 |
| 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: ‑ |