Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,013,998 | T→C | 51.8% | A2A (GCA→GCG) | exbD ← | membrane spanning protein in TonB‑ExbB‑ExbD complex |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,013,998 | 0 | T | C | 51.8% | 3.2 / 72.3 | 83 | A2A (GCA→GCG) | exbD | membrane spanning protein in TonB‑ExbB‑ExbD complex |
Reads supporting (aligned to +/- strand): ref base T (20/20); new base C (34/9); total (54/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.46e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 |
GATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTG > minE/2013932‑2014061 | gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:1137346/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:1973233/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:1968940/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:2429462/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:2446569/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:2675969/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:1664826/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:3220841/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:1514358/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:517858/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:722471/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:117405/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:699122/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:594665/1‑69 (MQ=255) gATAAACGGCGTCACGTTGATATAATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATCGc > 1:538436/1‑69 (MQ=255) cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta < 1:1856769/69‑1 (MQ=255) cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta < 1:1916164/69‑1 (MQ=255) cACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta < 1:2033038/69‑1 (MQ=255) aCGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtattat > 1:607003/1‑69 (MQ=255) tCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGt < 1:552137/69‑1 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:577900/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:3005809/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:2973099/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:3416119/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:377021/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:2285271/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:2196410/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:521292/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:877743/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:2185360/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:1645949/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatta > 1:1307452/1‑50 (MQ=255) tCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGAtatt > 1:3139784/1‑49 (MQ=255) gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc > 1:2319648/1‑69 (MQ=255) gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc > 1:1322779/1‑69 (MQ=255) gtcgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAAc > 1:322360/1‑69 (MQ=255) cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg < 1:2249768/69‑1 (MQ=255) cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg < 1:692533/69‑1 (MQ=255) cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg < 1:2439881/69‑1 (MQ=255) cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg < 1:3423653/69‑1 (MQ=255) cgAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg < 1:3407145/69‑1 (MQ=255) gAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCg < 1:1876626/68‑1 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACg > 1:1626561/1‑62 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:2850078/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1347871/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:2886301/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1658140/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:3131535/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:3162929/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:2843410/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:51551/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1694793/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1436332/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:2556039/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:2551921/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1341580/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:234267/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1750292/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:2236872/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTAATTTTTGTGCGACACGa > 1:1001794/1‑63 (MQ=255) aGGTTTTCGTTAAGATGCATCGCCATCGGAAATTATc > 1:3025965/1‑37 (MQ=37) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:311711/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:1662003/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:1140723/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:1189063/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:379593/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:228451/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:2679952/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:2587099/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:137561/68‑1 (MQ=255) ggTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGa < 1:2562140/68‑1 (MQ=255) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:2139418/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:76530/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:122149/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:243333/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:467569/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:1563156/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:2889443/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:2955377/41‑1 (MQ=25) ttCGTTAAGATGCATCGCCATCGGACATTATCCTGCACGTa < 1:2143550/41‑1 (MQ=25) tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg > 1:1731907/1‑69 (MQ=255) tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg > 1:1696644/1‑69 (MQ=255) tGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTg > 1:796282/1‑69 (MQ=255) | GATAAACGGCGTCACGTTGATATCATGCATTTCGCCGTTATCGTCGAGGTTTTCGTTAAGATGCATTGCCATCGGATATTATCCTGCGCGTAATTTTTGTGCGACACGAACCGGATGCGCAGCGGCGCTG > minE/2013932‑2014061 |
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: ‑ |