Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 600,547 | T→C | 22.1% | C145C (TGT→TGC) | 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,547 | 0 | T | C | 22.1% | 93.8 / 20.0 | 72 | C145C (TGT→TGC) | dmsB | dimethyl sulfoxide reductase, anaerobic, subunit B |
Reads supporting (aligned to +/- strand): ref base T (36/20); new base C (10/6); total (46/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01 |
TGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT > minE/600480‑600598 | tGACCAAATGGGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc < 1:2575951/69‑1 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:104665/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc < 1:984901/69‑1 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:732957/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:725863/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:375619/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:3005265/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:2778517/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:2748287/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc < 1:2648700/69‑1 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:2540920/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:1913020/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:1652910/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:1009705/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:1288943/1‑69 (MQ=255) tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc > 1:1484668/1‑69 (MQ=255) gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg > 1:2721789/1‑69 (MQ=255) gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg > 1:782895/1‑69 (MQ=255) gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg > 1:2644615/1‑69 (MQ=255) tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc > 1:566454/1‑67 (MQ=255) tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg > 1:1511350/1‑69 (MQ=255) tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg > 1:1386871/1‑69 (MQ=255) gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg < 1:1778081/68‑1 (MQ=255) gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg < 1:1613414/68‑1 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGt > 1:2549692/1‑44 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:811681/1‑68 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:1443125/1‑68 (MQ=255) aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc > 1:234379/1‑68 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:3182408/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:1825213/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:3057055/69‑1 (MQ=255) gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:1454588/69‑1 (MQ=255) gtTGCGGAAGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc < 1:2990846/69‑1 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:148853/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:817901/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:937652/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:1787895/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:372127/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:3034058/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:1918571/1‑68 (MQ=255) aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg > 1:3142960/1‑68 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2944792/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:638323/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:1654183/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:1677454/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2563538/1‑69 (MQ=255) aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa > 1:2137095/1‑69 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:1201173/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:660962/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2517649/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2293559/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:3245806/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2056622/69‑1 (MQ=255) ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:2052734/69‑1 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGCTTTCGGACCTATCGACGAGCTGCGTAAAAAACa > 1:2233302/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGAt > 1:2156433/1‑36 (MQ=19) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:1741115/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:976327/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:3232381/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:3133524/1‑68 (MQ=255) aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac > 1:1206470/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) cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa < 1:2754409/62‑1 (MQ=255) cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa < 1:1151059/62‑1 (MQ=255) ttGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa < 1:3000979/59‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1202253/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2793834/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1324400/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:1734795/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2467035/52‑1 (MQ=255) ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa < 1:2096578/52‑1 (MQ=255) | TGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT > minE/600480‑600598 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |