| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,314 | A→T | 41.2% | Q591L (CAG→CTG) ‡ | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,314 | 0 | A | T | 41.2% | 32.1 / 57.8 | 51 | Q591L (CAG→CTG) ‡ | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base A (14/16); new base T (8/13); total (22/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.79e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
CGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTC > minE/833246‑833378 | cGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGc > 1:5440272/1‑71 (MQ=255) tCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGcccc < 1:3377181/71‑3 (MQ=255) aaagCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTg < 1:616480/69‑1 (MQ=255) agCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGc < 1:4208992/71‑1 (MQ=255) gCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGcc < 1:2457101/71‑2 (MQ=255) cGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGcct < 1:2909029/71‑3 (MQ=255) tggcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTg < 1:244647/70‑1 (MQ=255) gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg < 1:2734165/71‑2 (MQ=255) gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg < 1:1126363/71‑2 (MQ=255) gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg < 1:590398/71‑2 (MQ=255) gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg > 1:4792371/1‑70 (MQ=255) gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGc > 1:96120/1‑70 (MQ=255) ggcTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCt > 1:5004806/1‑71 (MQ=255) gcTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt < 1:3890155/71‑1 (MQ=255) cTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCtcg < 1:1484523/71‑3 (MQ=255) cTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCtc > 1:385212/1‑69 (MQ=255) ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc < 1:1988872/70‑1 (MQ=255) ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc < 1:21541/70‑1 (MQ=255) ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACt > 1:2654840/1‑71 (MQ=255) cacACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTgg < 1:1095538/69‑1 (MQ=255) cacACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTgg < 1:1788640/69‑1 (MQ=255) gcACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTgg > 1:3041872/1‑70 (MQ=255) acACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTgg > 1:5470204/2‑70 (MQ=255) aCTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTggaac < 1:1860594/71‑4 (MQ=255) aCTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCtt > 1:1999443/1‑62 (MQ=255) aCTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAg < 1:2991167/71‑1 (MQ=255) tATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCg > 1:1920352/1‑71 (MQ=255) tAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTcc < 1:1339780/70‑1 (MQ=255) gAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGgc < 1:3114799/65‑1 (MQ=255) gAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCtc > 1:3578442/1‑47 (MQ=16) aaaaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTcaca > 1:2860713/1‑71 (MQ=255) aaaCTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCtc < 1:922362/45‑2 (MQ=12) aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:3706132/71‑1 (MQ=255) aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:5646108/71‑1 (MQ=255) aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:3900721/71‑1 (MQ=255) aaCTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATg < 1:1513747/71‑1 (MQ=255) aCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTc > 1:5407000/1‑56 (MQ=255) gTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAg > 1:2884735/1‑70 (MQ=255) tCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTca < 1:2669784/60‑1 (MQ=255) tCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCa > 1:327492/1‑71 (MQ=255) cacAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAAcaccac > 1:657669/3‑70 (MQ=255) cGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCAct > 1:2537324/1‑71 (MQ=255) acAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCAcc > 1:5181827/2‑69 (MQ=255) aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:1340339/59‑1 (MQ=255) aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct < 1:5119343/71‑1 (MQ=255) aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct < 1:5112745/71‑1 (MQ=255) aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct < 1:1205984/71‑1 (MQ=255) aCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGaa > 1:1047678/1‑59 (MQ=255) aCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTa > 1:697915/1‑71 (MQ=255) gCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTc > 1:4419780/1‑71 (MQ=255) gCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTc > 1:17889/1‑71 (MQ=255) | CGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTC > minE/833246‑833378 |
| 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: ‑ |