| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,297 | T→C | 47.1% | R585R (CGT→CGC) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,297 | 0 | T | C | 47.1% | 7.8 / 63.2 | 51 | R585R (CGT→CGC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base T (11/16); new base C (6/18); total (17/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.72e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-01 | |||||||||||
TTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCA > minE/833231‑833367 | tttCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCg > 1:1400618/1‑70 (MQ=255)cttCTGCGGTTGTGTCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGccca > 1:949659/2‑69 (MQ=255) tCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCa < 1:5222416/70‑1 (MQ=255) gtgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg > 1:5529711/1‑71 (MQ=255) tgtCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg < 1:1003885/70‑1 (MQ=255) tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTggt < 1:4121250/71‑2 (MQ=255) tgtCGGGCAAAGCGGTGGCGGCTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTgg < 1:1119725/70‑1 (MQ=255) tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa < 1:1294849/59‑1 (MQ=255) tCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCaa < 1:2999885/59‑1 (MQ=255) 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:4792371/1‑70 (MQ=255) gcggcTGGGCACACTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGcg < 1:590398/71‑2 (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) ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTgg > 1:907848/1‑46 (MQ=255) ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc < 1:21541/70‑1 (MQ=255) ggCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGAc < 1:1988872/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) cTATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTg < 1:4761113/39‑1 (MQ=21) tATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTg < 1:4127228/38‑1 (MQ=21) tATGTCGGCCAGGAAAAACTGCGCCCACAAACCGGCTg < 1:4313406/38‑1 (MQ=21) tATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCg > 1:1920352/1‑71 (MQ=255) tGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTg < 1:1626937/36‑1 (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) | TTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATGTAGGTCAGGAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCA > minE/833231‑833367 |
| 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: ‑ |