| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,259,787 | A→G | 27.3% | G379G (GGT→GGC) | torZ ← | trimethylamine N‑oxide reductase system III, catalytic subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,259,787 | 0 | A | G | 27.3% | 79.5 / 33.7 | 55 | G379G (GGT→GGC) | torZ | trimethylamine N‑oxide reductase system III, catalytic subunit |
| Reads supporting (aligned to +/- strand): ref base A (16/24); new base G (10/5); total (26/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.29e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01 | |||||||||||
GCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGC > minE/1259718‑1259854 | gccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAc > 1:5665775/1‑71 (MQ=255) ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTAcc < 1:1749692/71‑1 (MQ=255) cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa < 1:1386230/71‑1 (MQ=255) cgccAACACGCGGCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa < 1:5047619/71‑1 (MQ=255) cAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATc > 1:2102408/1‑71 (MQ=255) aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCt > 1:109426/1‑71 (MQ=255) cgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCt < 1:1088197/70‑1 (MQ=255) cgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa < 1:1946550/71‑1 (MQ=255) tggggTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCaa > 1:3184431/3‑70 (MQ=255) tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTa < 1:742430/69‑1 (MQ=255) gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTa > 1:4083805/1‑54 (MQ=255) gTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACAttg > 1:1004409/1‑69 (MQ=255) ttaccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc > 1:2755281/4‑71 (MQ=255) ttGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgc > 1:4433378/1‑71 (MQ=255) taccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgc > 1:239651/3‑70 (MQ=255) ccaccaTTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGcc < 1:5319517/50‑1 (MQ=255) ggAGTAATTATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa < 1:4634668/43‑1 (MQ=255) gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc < 1:2844667/71‑1 (MQ=255) gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc < 1:3032831/71‑1 (MQ=255) aGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc > 1:4600250/1‑70 (MQ=255) aaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGc < 1:3618007/62‑1 (MQ=255) aaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGAcga < 1:419678/69‑3 (MQ=255) tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCa < 1:1196808/36‑1 (MQ=255) gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc < 1:138006/70‑1 (MQ=255) gtAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATcc < 1:4525914/70‑1 (MQ=255) tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc > 1:3898621/1‑69 (MQ=255) gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc < 1:1149535/63‑1 (MQ=255) gagaAACCAAAACCGCCGCACGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc < 1:248870/63‑1 (MQ=255) gaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACAt > 1:3003638/1‑42 (MQ=255) gaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt > 1:2201941/1‑71 (MQ=255) agaCCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCa > 1:402649/3‑65 (MQ=255) gaCCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGtttt > 1:4128930/2‑70 (MQ=255) aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGc < 1:3615561/50‑1 (MQ=255) aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCa < 1:908069/60‑1 (MQ=255) aaCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc < 1:2122035/71‑1 (MQ=255) aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc > 1:2733240/1‑70 (MQ=255) ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCa > 1:2902326/1‑55 (MQ=255) cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATtgct > 1:3988192/1‑40 (MQ=12) cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGtttt > 1:3137196/1‑67 (MQ=255) cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc < 1:1660741/60‑1 (MQ=255) cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt < 1:3178490/63‑1 (MQ=255) aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa > 1:5384322/1‑71 (MQ=255) aaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa > 1:2015270/1‑71 (MQ=255) cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt < 1:317868/65‑1 (MQ=255) gccgccGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCa > 1:2928180/1‑70 (MQ=255) gccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa > 1:3376499/1‑70 (MQ=255) cgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg < 1:5535675/61‑1 (MQ=255) gccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt < 1:2246028/57‑1 (MQ=255) gccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt < 1:4789950/71‑1 (MQ=255) ccgggTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCa > 1:3004314/4‑63 (MQ=255) cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg > 1:2805669/1‑71 (MQ=255) cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGGATCCAGTGTTTCTGTTCACCATATTGCt < 1:3908737/71‑1 (MQ=255) cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAg > 1:1119392/1‑49 (MQ=255) tgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc > 1:2667483/1‑71 (MQ=255) tACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg < 1:1224449/52‑1 (MQ=255) tACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg < 1:3432571/52‑1 (MQ=255) | GCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGC > minE/1259718‑1259854 |
| 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: ‑ |