Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,259,319 | T→G | 6.7% | Q535H (CAA→CAC) ‡ | 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,319 | 0 | T | G | 6.7% | 194.2 / 7.8 | 75 | Q535H (CAA→CAC) ‡ | torZ | trimethylamine N‑oxide reductase system III, catalytic subunit |
Reads supporting (aligned to +/- strand): ref base T (41/29); new base G (2/3); total (43/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.45e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01 |
TCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCA > minE/1259249‑1259388 | tctcTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAAt > 1:687857/1‑71 (MQ=255) ctctTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAAtt > 1:2357291/1‑71 (MQ=255) ttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGc < 1:1408938/68‑1 (MQ=255) ttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCg > 1:3826915/1‑69 (MQ=255) tCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGa < 1:556606/70‑1 (MQ=255) cGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCa < 1:1027478/71‑1 (MQ=255) cAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCaa > 1:167296/1‑69 (MQ=255) aGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACa > 1:1735057/1‑70 (MQ=255) gTTTGAGTAATTCCGCAAGATCGGCAAACACGACAAAATCGTTACGCGCTTCAAAt < 1:1053596/56‑1 (MQ=255) tttGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCt < 1:1595917/71‑1 (MQ=255) tGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGc > 1:3257345/1‑71 (MQ=255) tGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGc > 1:1436213/1‑71 (MQ=255) gAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCt > 1:3658702/1‑71 (MQ=255) aGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCaa < 1:381077/61‑1 (MQ=255) tAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGc < 1:1690565/57‑1 (MQ=255) tAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGc > 1:1026603/1‑63 (MQ=255) tAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCa < 1:3612087/71‑1 (MQ=255) tAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCa > 1:3782297/1‑71 (MQ=255) aTTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATc > 1:80364/1‑71 (MQ=255) ttCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCg < 1:943781/71‑1 (MQ=255) ttCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCg > 1:430903/1‑71 (MQ=255) tCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCgg > 1:798592/1‑71 (MQ=255) cGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTg < 1:1684399/60‑1 (MQ=255) cGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGc > 1:68286/1‑61 (MQ=255) gCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATc > 1:2846398/1‑66 (MQ=255) cAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGcaca < 1:1373352/71‑1 (MQ=255) cAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGcac > 1:3463489/1‑70 (MQ=255) gATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAg < 1:3289119/51‑1 (MQ=255) tCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCAt < 1:1929874/59‑1 (MQ=255) tCGGCAAACACGGCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATg < 1:1939879/71‑1 (MQ=255) cGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAAc > 1:1313376/1‑46 (MQ=255) cGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCg > 1:3119788/1‑60 (MQ=255) ggCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGc < 1:3449898/70‑1 (MQ=255) gCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTg < 1:1539223/71‑1 (MQ=255) gCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTg > 1:576774/1‑71 (MQ=255) cAAACACGTCAAAATCGTTACGCGCTTCAAATTGCg < 1:3364963/36‑1 (MQ=255) cAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTg > 1:2655832/1‑49 (MQ=255) aaCACGTCAAAATCGTTACGCGCTTCAAATTGCGGa > 1:99912/1‑36 (MQ=255) aaCACGTCAAAATCGTTACGCGCTTCAAATTGCGGa > 1:2197614/1‑36 (MQ=255) acacGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTg < 1:2929217/71‑1 (MQ=255) cacGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGc < 1:544249/71‑1 (MQ=255) gTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCt > 1:2098813/1‑69 (MQ=255) gTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTg > 1:3067295/1‑70 (MQ=255) gTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGt < 1:504385/71‑1 (MQ=255) atCAAAATCATTACGCGCTTCATAGCGTGGCGGCACCACTTGCTTCATCGGCACCAGATGCTGATTACTGt < 1:1807957/70‑1 (MQ=255) tCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGc > 1:1157128/1‑38 (MQ=255) tCAAAATCATTACGCGCTTCATAGCGTGGCGGCACCACTTGCTTCATCGGCACCAGATGCTGATTACTGTa < 1:188086/71‑1 (MQ=255) cAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGc > 1:44713/1‑59 (MQ=255) cAAAATCATTACGCGCTTCATAGCGTGGCGGCACCACTTGCTTCATCGGcacc > 1:903052/1‑52 (MQ=255) cAAAATCATTACGCGCTTCATAGCGTGGCGGCACCACTTGCTTCATCGGCACCAGATGCTGATTACTGTaa > 1:165247/1‑71 (MQ=255) aaaaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAAt < 1:2587043/71‑1 (MQ=255) aaaaTCATTACGCGCTTCATAGCGTGGCGGCACCACTTGCTTCATCGGCACCAGATGCTGATTACTGTAAt < 1:337808/71‑1 (MQ=255) aaaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATg > 1:2273071/1‑56 (MQ=255) aaaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGc > 1:1347719/1‑57 (MQ=255) aaaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAAt > 1:1249524/1‑70 (MQ=255) aaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACaa < 1:2557413/51‑1 (MQ=255) aTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCAc < 1:1778016/71‑1 (MQ=255) aTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCAc < 1:470536/71‑1 (MQ=255) tCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCAcc > 1:3372617/1‑71 (MQ=255) tCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCAcc > 1:3105599/1‑71 (MQ=255) tCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCAcc > 1:513725/1‑71 (MQ=255) cGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCg > 1:3768493/1‑71 (MQ=255) cGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCg > 1:3200162/1‑71 (MQ=255) cGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCg > 1:151080/1‑71 (MQ=255) gTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCgg > 1:455261/1‑71 (MQ=255) gTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCg > 1:128098/1‑70 (MQ=255) gTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCg > 1:746990/1‑70 (MQ=255) ttACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATg < 1:2040405/50‑1 (MQ=255) aCGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTc < 1:2202681/70‑1 (MQ=255) cgcgCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCAt < 1:1215971/71‑1 (MQ=255) gcgcTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATc > 1:3400514/1‑71 (MQ=255) gcTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGt > 1:1666713/1‑71 (MQ=255) gcTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGt > 1:3452958/1‑71 (MQ=255) cTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATc < 1:397559/58‑1 (MQ=255) tCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGc > 1:1738389/1‑49 (MQ=255) aaaTTGCTGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAg < 1:921377/70‑1 (MQ=255) aaTTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAg > 1:1433139/1‑69 (MQ=255) aTTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCa < 1:1517939/71‑1 (MQ=255) | TCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCA > minE/1259249‑1259388 |
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: ‑ |