| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,826,261 | C→A | 20.6% | A21A (GCG→GCT) | hydN ← | formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,826,261 | 0 | C | A | 20.6% | 208.1 / 27.9 | 82 | A21A (GCG→GCT) | hydN | formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit |
| Reads supporting (aligned to +/- strand): ref base C (37/28); new base A (13/4); total (50/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.71e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01 | |||||||||||
GATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTC > minE/1826192‑1826330 | gATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCAc < 1:1998606/70‑1 (MQ=255) cgcgGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTAc < 1:1741691/71‑1 (MQ=255) ggTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACtt > 1:2227259/1‑70 (MQ=255) ggTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACtt > 1:2232043/1‑70 (MQ=255) ggTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACtt > 1:1351630/1‑70 (MQ=255) tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTcac < 1:1640991/71‑1 (MQ=255) tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTca > 1:1351575/1‑70 (MQ=255) tAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTcac < 1:3690796/71‑1 (MQ=255) aaaaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTca > 1:3104207/1‑69 (MQ=255) aaaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAgg > 1:209274/1‑71 (MQ=255) aaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg > 1:2554580/1‑70 (MQ=255) aaaGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAgg > 1:1732410/1‑70 (MQ=255) tttCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACgg < 1:3952329/71‑1 (MQ=255) cGGGGGCAGCGATGCACAGTCCTGATTTTCCAGATGAGAAACCGCGCACGCTACTTc < 1:3913386/57‑1 (MQ=38) gggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCaa > 1:892809/1‑68 (MQ=255) gggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCa < 1:3188542/71‑1 (MQ=255) tCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATa < 1:1548548/71‑1 (MQ=255) tCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATa < 1:1551969/71‑1 (MQ=255) aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACa < 1:3711314/71‑1 (MQ=255) aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGTCAACCAATAc < 1:257914/70‑1 (MQ=255) aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa < 1:2733904/71‑1 (MQ=255) aGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa > 1:456722/1‑71 (MQ=255) aGCGATGCAAAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACa > 1:3773237/1‑71 (MQ=255) cGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACAtt > 1:1524184/1‑71 (MQ=255) aTGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGTAACCAATACATTTa > 1:1004968/1‑71 (MQ=255) tGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTAc > 1:3501009/1‑71 (MQ=255) tGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTAc < 1:2996749/71‑1 (MQ=255) cacaGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACg < 1:3234888/52‑1 (MQ=255) cattCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTc < 1:3592228/66‑1 (MQ=255) aatCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTcgca < 1:3886121/69‑2 (MQ=255) tCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTc < 1:2760197/71‑1 (MQ=255) tCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAAccga > 1:3632359/1‑53 (MQ=255) tCCTGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAAccga > 1:457613/1‑53 (MQ=255) tGATTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg > 1:3303453/1‑68 (MQ=255) gATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGt > 1:1131247/1‑40 (MQ=255) aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcgc > 1:3706805/1‑63 (MQ=255) aTTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCaa > 1:2801771/1‑71 (MQ=255) aTTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCaa > 1:1252820/1‑71 (MQ=255) aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATccg > 1:2860902/1‑66 (MQ=255) aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGc < 1:1301137/69‑1 (MQ=255) aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaa > 1:2380664/1‑71 (MQ=255) aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaa > 1:2018472/1‑71 (MQ=255) aTTCTCATGATGCGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCaa > 1:2167744/1‑71 (MQ=255) ttttCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAat < 1:1880420/71‑1 (MQ=255) tttCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcg > 1:1066757/1‑60 (MQ=255) tttCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTcg > 1:783581/1‑60 (MQ=255) tttCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTc < 1:1525199/64‑1 (MQ=255) ccTGATGAGAAACCGCGCACGCTACTTCACAGGTACg > 1:1191393/1‑37 (MQ=255) cTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGc < 1:184737/63‑1 (MQ=255) cTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAAtgatga < 1:585801/71‑1 (MQ=255) tgatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaa > 1:3912037/1‑71 (MQ=255) gatgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac > 1:829433/1‑70 (MQ=255) atgcGACACTGCGCAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATaaac < 1:3390769/70‑2 (MQ=255) atgaGAAACCGCGCACGCTACTTCACAGGTACGGCAAc > 1:3402681/1‑38 (MQ=255) atgaGAAACCGCGCACGCTACTTCACAGGTACGGCAAc > 1:1973530/1‑38 (MQ=255) tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACAttt < 1:3449270/47‑1 (MQ=255) tgaGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACg > 1:3370607/1‑71 (MQ=255) gagaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCaa > 1:1160774/1‑60 (MQ=255) agaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt > 1:1502972/1‑71 (MQ=255) agaAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt < 1:3194102/71‑1 (MQ=255) gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat > 1:529211/1‑62 (MQ=255) gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat > 1:2146661/1‑62 (MQ=255) gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAatgat > 1:3120871/1‑62 (MQ=255) gaAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGt > 1:2359820/1‑70 (MQ=255) aaaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGtt > 1:2441630/1‑70 (MQ=255) aaaCCACGCACGCTACTTCACAGGTACGGCAACCAATACa < 1:3453150/40‑1 (MQ=255) aaaCCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTc < 1:2575076/47‑1 (MQ=255) aaCCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAAtgatga < 1:3015079/62‑1 (MQ=255) actgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcat < 1:1392119/68‑1 (MQ=255) aCCACGCACGCTACTTCACAGGTACGGCAACCAATACAtt > 1:777117/1‑40 (MQ=255) ccGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcat > 1:221195/1‑70 (MQ=255) tgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCa > 1:729198/2‑52 (MQ=255) tgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCa > 1:3328465/2‑52 (MQ=255) tgcgcAAGCCACTTCACAGGTACGGCAACCGATACATTTCGTCGCATCCGCAATAATAAACCGGTtcatga > 1:3730334/2‑69 (MQ=255) gcgcAAGCCACTTCACAGGTACGGCAACCGATACa > 1:571261/1‑35 (MQ=25) cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcat > 1:135977/1‑67 (MQ=255) cgcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCaaa < 1:201223/71‑1 (MQ=255) gcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAt < 1:384788/71‑1 (MQ=255) gcACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAt > 1:2888500/1‑71 (MQ=255) cACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAAtt < 1:2533518/71‑1 (MQ=255) aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTtcat > 1:850961/1‑64 (MQ=255) aCGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc > 1:3439196/1‑71 (MQ=255) aCGCTACTTCACAAGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATtc < 1:768890/71‑1 (MQ=255) | GATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCGCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTC > minE/1826192‑1826330 |
| 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: ‑ |