Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 833,372 | T→G | 18.5% | S610S (TCT→TCG) | narG → | nitrate reductase 1, alpha subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 833,372 | 0 | T | G | 18.5% | 96.4 / 24.1 | 54 | S610S (TCT→TCG) | narG | nitrate reductase 1, alpha subunit |
Reads supporting (aligned to +/- strand): ref base T (27/17); new base G (5/5); total (32/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01 |
AAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGAT > minE/833302‑833440 | aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct < 1:5119343/71‑1 (MQ=255) aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct < 1:1205984/71‑1 (MQ=255) aaaCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCActtct < 1:5112745/71‑1 (MQ=255) aCCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTa > 1:697915/1‑71 (MQ=255) gCTGGTTGCCACTGGCCTTTGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTc > 1:4419780/1‑71 (MQ=255) gCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTc > 1:17889/1‑71 (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) ttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagccag > 1:926560/1‑71 (MQ=255) ttGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTccagccag > 1:3296858/1‑71 (MQ=255) tGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc < 1:2288136/68‑1 (MQ=255) tGCGCTCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagcc < 1:1474575/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:674545/1‑71 (MQ=255) cAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAAc > 1:416713/1‑71 (MQ=255) aGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAAcc < 1:3192181/71‑1 (MQ=255) tCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa < 1:1168531/64‑1 (MQ=255) tCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGaa < 1:3531251/64‑1 (MQ=255) caccgcgCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGaaa < 1:4525448/60‑1 (MQ=255) cGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTa > 1:4357240/1‑58 (MQ=255) ggCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcgg > 1:114987/1‑71 (MQ=255) gcgcCAAATGAACAGCACCTCGTTTTTCTACAATCATTccagccaa > 1:4170556/1‑45 (MQ=255) gcgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTAt < 1:2253670/55‑1 (MQ=255) gcgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcggcg > 1:1445627/1‑70 (MQ=255) cgTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGa < 1:3839467/71‑1 (MQ=255) tCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAag < 1:2050216/71‑1 (MQ=255) cacaTGAACAGCACTTCTTATTTCTATAACCACTccag < 1:5281818/38‑1 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTccagcca > 1:5628058/1‑40 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTccagcca > 1:5499084/1‑40 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTccagcca > 1:547798/1‑40 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTccagcca > 1:1170920/1‑40 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTccagcca > 1:1177539/1‑40 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAAc > 1:3051946/1‑55 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGaa > 1:2999371/1‑68 (MQ=255) acaTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGaa > 1:5126154/1‑68 (MQ=255) caTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCAcg < 1:3812713/61‑1 (MQ=255) caTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGt < 1:3012132/71‑1 (MQ=255) aaTGAACAGCACCTCGTTTTTCTACAATCATTccagcca > 1:2136937/2‑39 (MQ=12) aTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGtt < 1:1978814/71‑1 (MQ=255) aTGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCCTTATGaa < 1:657668/51‑1 (MQ=255) aCAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTg < 1:1617511/71‑1 (MQ=255) aGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTc < 1:2841665/71‑1 (MQ=255) gCACTTCTTATTTCTATAACCACTACAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGtt > 1:4711858/1‑69 (MQ=255) aCTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCg > 1:3057029/1‑46 (MQ=255) cttcttATTTCTATAACCACTCCAGCCAGTGGCGt < 1:3248070/35‑1 (MQ=255) ttcttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGa > 1:1530185/1‑71 (MQ=255) tcttATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAag < 1:113318/55‑1 (MQ=255) tcttATTTCTATAACCACTCCAGCCAGTGGCGATATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGAt > 1:2557337/1‑71 (MQ=255) | AAACCGGCTGGCAGCCGCTGGCGTTTGCCCTTGACTGGCAGCGTCCGGCGCGTCACATGAACAGCACTTCTTATTTCTATAACCACTCCAGCCAGTGGCGTTATGAAACCGTCACGGCGGAAGAGTTGCTGTCACCGAT > minE/833302‑833440 |
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: ‑ |