Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 600,553 | G→A | 28.4% | L147L (CTG→CTA) | 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,553 | 0 | G | A | 28.4% | 58.6 / 20.3 | 56 | L147L (CTG→CTA) | dmsB | dimethyl sulfoxide reductase, anaerobic, subunit B |
Reads supporting (aligned to +/- strand): ref base G (23/17); new base A (10/6); total (33/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.73e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT > minE/600495‑600598 | gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg > 1:2644615/1‑69 (MQ=255) gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg > 1:2721789/1‑69 (MQ=255) gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg > 1:782895/1‑69 (MQ=255) tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc > 1:566454/1‑67 (MQ=255) tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg > 1:1386871/1‑69 (MQ=255) tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg > 1:1511350/1‑69 (MQ=255) gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg < 1:1778081/68‑1 (MQ=255) gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg < 1:1613414/68‑1 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:234379/1‑68 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:1443125/1‑68 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:811681/1‑68 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:1454588/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:3057055/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:1825213/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:3182408/69‑1 (MQ=255) gtTGCGGAAGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:2990846/69‑1 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:1918571/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:372127/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:937652/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:817901/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:3142960/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:148853/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:1787895/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:3034058/1‑68 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2944792/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2563538/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:638323/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2137095/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:1677454/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:1654183/1‑69 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2293559/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:3245806/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2517649/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) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGAt > 1:2156433/1‑36 (MQ=19) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:976327/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:3232381/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:3133524/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:2829044/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:2440226/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:2142610/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:1741115/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:1206470/1‑68 (MQ=255) cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa < 1:2754409/62‑1 (MQ=255) cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa < 1:1151059/62‑1 (MQ=255) ttGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:3000979/59‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2793834/52‑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:1202253/52‑1 (MQ=255) | GTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT > minE/600495‑600598 |
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: ‑ |