| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,259,805 | C→T | 16.4% | A373A (GCG→GCA) | 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,805 | 0 | C | T | 16.4% | 149.9 / 19.4 | 61 | A373A (GCG→GCA) | torZ | trimethylamine N‑oxide reductase system III, catalytic subunit |
| Reads supporting (aligned to +/- strand): ref base C (23/28); new base T (6/4); total (29/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.96e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 | |||||||||||
GCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTG > minE/1259739‑1259872 | gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccg > 1:1662212/1‑71 (MQ=255)gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccg > 1:3842005/1‑71 (MQ=255) ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAg < 1:744832/71‑1 (MQ=255) ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTc < 1:1265038/70‑1 (MQ=255) gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCa < 1:1543201/70‑1 (MQ=255) gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc < 1:423356/71‑1 (MQ=255) gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc > 1:311156/1‑71 (MQ=255) aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc < 1:3902408/71‑1 (MQ=255) tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa > 1:1484681/1‑55 (MQ=255) gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc < 1:1773659/63‑1 (MQ=255) aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt > 1:1623018/1‑65 (MQ=255) ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTATCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg < 1:1866248/71‑1 (MQ=255) cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAg < 1:1463004/46‑1 (MQ=16) cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt > 1:237533/1‑71 (MQ=255) cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt > 1:1587592/1‑71 (MQ=255) aaaaCCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccgc > 1:855298/1‑42 (MQ=255) aaaCCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgct > 1:1019307/1‑68 (MQ=255) aaaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa > 1:3627077/1‑43 (MQ=255) cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc > 1:2249374/1‑46 (MQ=255) cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa < 1:2071491/71‑1 (MQ=255) cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa > 1:956166/1‑71 (MQ=255) ccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc > 1:1219358/1‑44 (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:3121549/1‑65 (MQ=255) tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg > 1:3416274/1‑65 (MQ=255) gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc < 1:3303665/44‑2 (MQ=12) gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc < 1:3576895/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:3514365/71‑1 (MQ=255) ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa < 1:843594/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:3115781/35‑1 (MQ=255) aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt < 1:1159233/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) | GCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTG > minE/1259739‑1259872 |
| 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: ‑ |