| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 833,349 | G→C | 19.5% | A603P (GCG→CCG) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 833,349 | 0 | G | C | 19.5% | 104.7 / 18.5 | 51 | A603P (GCG→CCG) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (25/16); new base C (5/5); total (30/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.22e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGG > minE/833286‑833418 | gAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGgc < 1:3114799/65‑1 (MQ=255) aaaaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTcaca > 1:2860713/1‑71 (MQ=255) aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:3706132/71‑1 (MQ=255) aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:3900721/71‑1 (MQ=255) aaCTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATg < 1:5646108/71‑1 (MQ=255) aaCTGCGCCCACAAACCGGCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATg < 1:1513747/71‑1 (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) ccGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGTACAGc < 1:3527937/51‑1 (MQ=255) gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCa > 1:209608/1‑50 (MQ=255) gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCAc > 1:2019425/1‑51 (MQ=255) gCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCAc > 1:3417849/1‑51 (MQ=255) cGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagc < 1:1808947/71‑1 (MQ=255) gTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagcc < 1:5274879/71‑1 (MQ=255) ctttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc < 1:3783453/70‑1 (MQ=255) tttGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCa < 1:1780901/62‑1 (MQ=255) tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGaa > 1:5567355/1‑38 (MQ=255) tttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGaa > 1:4972207/1‑38 (MQ=255) ttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagccag > 1:926560/1‑71 (MQ=255) ttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagccag > 1:3296858/1‑71 (MQ=255) tGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc < 1:1474575/68‑1 (MQ=255) tGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc < 1:2288136/68‑1 (MQ=255) tGCCATTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTAt > 1:1832980/1‑49 (MQ=255) ccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGc < 1:1856097/71‑1 (MQ=255) ccTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGc > 1:4594606/1‑71 (MQ=255) ttGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGt > 1:4556130/1‑71 (MQ=255) tGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGt > 1:3218542/1‑70 (MQ=255) gACTGGCAGCGTCCGGCTCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagcca > 1:2850527/1‑62 (MQ=255) gACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCa > 1:436715/1‑53 (MQ=255) gACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagccaa > 1:3442579/1‑62 (MQ=255) gACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTa > 1:1517193/1‑71 (MQ=255) cTGGCAGCGTCCGGCGCGTCCCATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTAt > 1:706975/1‑70 (MQ=255) gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaaa < 1:4814578/71‑1 (MQ=255) gCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATATCTATAACCa > 1:2405054/1‑48 (MQ=255) cAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAAc > 1:416713/1‑71 (MQ=255) cAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAAc > 1:674545/1‑71 (MQ=255) aGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAAcc < 1:3192181/71‑1 (MQ=255) tCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa < 1:3531251/64‑1 (MQ=255) tCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa < 1:1168531/64‑1 (MQ=255) cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa > 1:4357240/1‑58 (MQ=255) ggCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcgg > 1:114987/1‑71 (MQ=255) | GAAAAACTGCGTCCGCAAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGG > minE/833286‑833418 |
| 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: ‑ |