| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,259,822 | G→T | 19.7% | L368M (CTG→ATG) ‡ | 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,822 | 0 | G | T | 19.7% | 200.2 / 26.3 | 76 | L368M (CTG→ATG) ‡ | torZ | trimethylamine N‑oxide reductase system III, catalytic subunit |
| Reads supporting (aligned to +/- strand): ref base G (29/32); new base T (8/7); total (37/39) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.77e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01 | |||||||||||
GATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTC > minE/1259756‑1259892 | gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc > 1:311156/1‑71 (MQ=255) aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc < 1:3902408/71‑1 (MQ=255) gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc < 1:1773659/63‑1 (MQ=255) aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt > 1:1623018/1‑65 (MQ=255) ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTATCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg < 1:1866248/71‑1 (MQ=255) cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt > 1:237533/1‑71 (MQ=255) cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt > 1:1587592/1‑71 (MQ=255) aaaCCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgct > 1:1019307/1‑68 (MQ=255) cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa > 1:956166/1‑71 (MQ=255) cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa < 1:2071491/71‑1 (MQ=255) cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc < 1:814624/49‑1 (MQ=255) cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt < 1:507695/70‑1 (MQ=255) cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt < 1:1074327/54‑1 (MQ=255) cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc > 1:589207/1‑70 (MQ=255) cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt > 1:1012164/1‑71 (MQ=255) cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt > 1:1824658/1‑71 (MQ=255) ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg < 1:1504740/71‑1 (MQ=255) ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg > 1:1156444/1‑71 (MQ=255) gtgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt > 1:2807651/1‑57 (MQ=255) gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg > 1:3950268/1‑71 (MQ=255) tACCAATCTGCCCTAACATCGCCTCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg < 1:1955581/70‑1 (MQ=255) ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt < 1:332172/52‑1 (MQ=255) ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt < 1:1372389/51‑1 (MQ=255) atttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt > 1:1223111/4‑71 (MQ=255) aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt > 1:233247/1‑42 (MQ=255) aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt < 1:2919854/49‑1 (MQ=255) aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt < 1:3436580/49‑1 (MQ=255) tCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa > 1:2617609/1‑50 (MQ=255) tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg > 1:3416274/1‑65 (MQ=255) tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg > 1:3121549/1‑65 (MQ=255) gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc < 1:3576895/44‑2 (MQ=12) gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc < 1:3303665/44‑2 (MQ=12) cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt < 1:3882319/53‑1 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt < 1:2709702/36‑1 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt < 1:1739990/36‑1 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt < 1:3317346/42‑1 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGacgc > 1:2492214/1‑60 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa < 1:843594/71‑1 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa < 1:3514365/71‑1 (MQ=255) tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg < 1:1012511/51‑1 (MQ=255) caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTg < 1:2524134/59‑1 (MQ=255) caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGc > 1:3033525/2‑71 (MQ=255) aaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCaaa > 1:2259757/1‑45 (MQ=12) aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt < 1:1163994/35‑1 (MQ=255) aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt < 1:1159233/35‑1 (MQ=255) aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt < 1:3115781/35‑1 (MQ=255) aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt > 1:3326263/1‑40 (MQ=255) aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTgca > 1:3921167/1‑58 (MQ=255) aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCt > 1:2520160/1‑71 (MQ=255) aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg > 1:3673769/1‑58 (MQ=255) cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt > 1:1182637/1‑35 (MQ=255) cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg < 1:2655954/71‑1 (MQ=255) gctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCa < 1:3094874/67‑1 (MQ=255) gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc > 1:3002658/1‑45 (MQ=255) gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc > 1:1904084/1‑45 (MQ=255) gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt > 1:2221078/1‑71 (MQ=255) gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt > 1:2216343/1‑71 (MQ=255) ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt < 1:2241161/68‑1 (MQ=255) ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt > 1:3379490/3‑70 (MQ=255) gccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAg > 1:1137738/1‑71 (MQ=255) gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg > 1:2827844/1‑71 (MQ=255) ccAGCCTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg < 1:1938514/63‑1 (MQ=255) cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa > 1:3621450/1‑71 (MQ=255) cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa < 1:3841247/71‑1 (MQ=255) cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa < 1:1910719/71‑1 (MQ=255) cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa < 1:1507521/71‑1 (MQ=255) cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa < 1:462589/71‑1 (MQ=255) cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa < 1:235307/71‑1 (MQ=255) aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt > 1:792378/1‑71 (MQ=255) gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg < 1:1881162/71‑2 (MQ=255) gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg < 1:480872/71‑2 (MQ=255) cGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg > 1:1084140/1‑71 (MQ=255) gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtgt < 1:415203/71‑1 (MQ=255) gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg < 1:371696/70‑1 (MQ=255) cACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc > 1:2726561/1‑71 (MQ=255) aGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGAtt > 1:3160833/1‑71 (MQ=255) gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc > 1:3191134/1‑71 (MQ=255) gCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTc < 1:2928067/71‑1 (MQ=255) | GATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGCGATTC > minE/1259756‑1259892 |
| 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: ‑ |