| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,586,248 | T→C | 52.7% | T343T (ACA→ACG) | mdtF ← | multidrug transporter, RpoS‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,586,248 | 0 | T | C | 52.7% | ‑2.0 / 29.9 | 38 | T343T (ACA→ACG) | mdtF | multidrug transporter, RpoS‑dependent |
| Reads supporting (aligned to +/- strand): ref base T (5/13); new base C (14/6); total (19/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.17e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAGGATAAACCGTCTTCAG > minE/2586188‑2586316 | gcGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGaa < 1:2921638/68‑1 (MQ=18)gcGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGaa < 1:2735621/68‑1 (MQ=18)gcGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGaa < 1:409547/68‑1 (MQ=18)gcGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGaa < 1:2458647/68‑1 (MQ=18)gcGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGaa < 1:1698671/68‑1 (MQ=18) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAg > 1:3070/1‑68 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAg > 1:521153/1‑68 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1058496/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:2530450/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1975862/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1792700/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1720448/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1705871/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1143074/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1332896/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:750595/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:711075/1‑69 (MQ=255) aaaCAGATACATGACCAGGAAGACGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTACTGAATAGa > 1:2648902/1‑69 (MQ=255) aaaCAGATACATGACCAGAAAGGCGAGGATGATAGCCTCAACCAGCGTTTTGAAAACTTCCTGAATAGa > 1:1523810/1‑69 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:1789365/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:3430073/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:2438821/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:608564/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:2021312/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:615388/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:19643/69‑1 (MQ=255) aGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAAtt < 1:1848019/69‑1 (MQ=255) gATACATGACCAGGAAGACGAGGATGATAGCCTCAAGCAGCGTTTTGaa < 1:1333247/49‑1 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGa < 1:3216816/38‑1 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGa < 1:1092345/38‑1 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGa < 1:2630052/38‑1 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGa < 1:2423263/38‑1 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGa < 1:1665799/38‑1 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAgg > 1:384735/1‑68 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAgg > 1:13962/1‑68 (MQ=255) atAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAgg > 1:1389028/1‑68 (MQ=255) cTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAGGATAAAc > 1:3051947/1‑69 (MQ=255) cTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAGGATAAAc > 1:888105/1‑69 (MQ=255) tgtTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAGGATAAACCGTCTTCAg < 1:2567144/69‑1 (MQ=255) tgtTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAGGATAAACCGTCTTCAg < 1:2304559/69‑1 (MQ=255) | ACGGAAATTCTGCAAAAACAGATACATGACCAGGAAGACTAGGATGATAGCCTCAACCAGTGTTTTGAAAACTTCCTGAATAGAAATTTCGATAAACGGCGTGGTGTCGTAAGGATAAACCGTCTTCAG > minE/2586188‑2586316 |
| 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: ‑ |