| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,381 | T→C | 11.1% | Y613Y (TAT→TAC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,381 | 0 | T | C | 11.1% | 142.6 / 9.9 | 54 | Y613Y (TAT→TAC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base T (23/25); new base C (4/2); total (27/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01 | |||||||||||
GCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAAT > minE/833312‑833451 | gCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat < 1:1040045/70‑1 (MQ=255) ttgCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACaa < 1:2131717/69‑1 (MQ=255) cGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATaa > 1:2773946/1‑67 (MQ=255) gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACt > 1:775826/1‑71 (MQ=255) tGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTc < 1:3227087/70‑1 (MQ=255) gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag > 1:229897/4‑71 (MQ=255) gcctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccag > 1:3170020/4‑71 (MQ=255) cGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagc < 1:254711/71‑1 (MQ=255) tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat < 1:2309997/57‑1 (MQ=255) tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat < 1:613701/57‑1 (MQ=255) tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTcca > 1:3372621/1‑65 (MQ=255) tGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt > 1:1590333/1‑71 (MQ=255) cccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCa < 1:184800/58‑1 (MQ=255) ttGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat < 1:832556/50‑1 (MQ=255) gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCtat > 1:22360/1‑48 (MQ=255) gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGtt < 1:1186222/70‑1 (MQ=255) gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa > 1:3239879/1‑71 (MQ=255) gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa > 1:2589941/1‑71 (MQ=255) gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa > 1:617901/1‑71 (MQ=255) gCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt < 1:1024171/56‑1 (MQ=255) cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt < 1:3222984/71‑1 (MQ=255) cGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGt < 1:1574199/71‑1 (MQ=255) gTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTc < 1:2144655/71‑1 (MQ=255) cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt > 1:2523068/1‑51 (MQ=255) cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGt > 1:685868/1‑51 (MQ=255) gcgcGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggc < 1:3085662/71‑1 (MQ=255) cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc < 1:529854/45‑1 (MQ=255) cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag > 1:2353032/1‑60 (MQ=255) cgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaaaag > 1:1769060/1‑60 (MQ=255) gcgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg < 1:1408708/71‑1 (MQ=255) tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa < 1:2789547/55‑1 (MQ=255) tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcgg > 1:4022660/1‑68 (MQ=255) tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg < 1:593061/64‑1 (MQ=255) tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAag < 1:3654739/71‑1 (MQ=255) cacaTGAACAGCACTTCTTATTTCTATAACCACTccagcc < 1:1225315/40‑1 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTcc > 1:1666908/1‑35 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTcc > 1:4035666/1‑35 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag < 1:3111514/71‑1 (MQ=255) aTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAaga < 1:3941187/68‑1 (MQ=255) tGAACAGCACTTCTTATTTCTATAACCACTccagc < 1:2052576/35‑1 (MQ=255) tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg < 1:4062324/59‑1 (MQ=255) tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGaa < 1:251319/65‑1 (MQ=255) tGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg > 1:2453286/1‑71 (MQ=255) aaCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa > 1:1258885/1‑44 (MQ=255) aCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa < 1:2891241/47‑1 (MQ=255) aGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc < 1:314178/71‑1 (MQ=255) gCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc < 1:481431/70‑1 (MQ=255) gCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCa < 1:2750492/71‑1 (MQ=255) cACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTg > 1:2171872/1‑64 (MQ=255) cACTTCATATTTCTATAACCACTCCAGCCAGTGGCGTTa > 1:536322/1‑39 (MQ=255) aCTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCAcc < 1:3387854/71‑1 (MQ=255) tcttATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa > 1:2905927/1‑40 (MQ=255) tcttATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa > 1:2813084/1‑40 (MQ=255) cttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATg > 1:635327/1‑71 (MQ=255) cttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATg > 1:837128/1‑71 (MQ=255) cttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATg > 1:1635795/1‑71 (MQ=255) ttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAagag > 1:391556/1‑55 (MQ=255) tttCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCg < 1:1692137/70‑1 (MQ=255) tCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGa > 1:3439424/1‑70 (MQ=255) taACTACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACa < 1:1296452/68‑1 (MQ=255) taACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACa < 1:3205636/68‑1 (MQ=255) taACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACa < 1:1038695/68‑1 (MQ=255) taACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAAt < 1:2571205/71‑1 (MQ=255) taACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTATCACCGATGGCGGACa < 1:792621/68‑1 (MQ=255) | GCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGATGGCGGACAAAT > minE/833312‑833451 |
| 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: ‑ |