| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 834,731 | T→C | 9.6% | R1063R (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 | 834,731 | 0 | T | C | 9.6% | 222.4 / 11.8 | 83 | R1063R (CGT→CGC) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base T (27/48); new base C (6/2); total (33/50) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.38e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
GTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAG > minE/834667‑834798 | gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg < 1:4491433/71‑1 (MQ=255) tcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc < 1:5683528/69‑1 (MQ=255) cAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCTGTGAAAGCCTGCTGGTTTATCGTccgcc < 1:3695165/68‑1 (MQ=255) agcagcTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATc < 1:136506/71‑1 (MQ=255) agcagcTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATc < 1:2868745/71‑1 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGt > 1:5052966/1‑62 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa > 1:3815079/1‑69 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa > 1:3327716/1‑69 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa > 1:1724886/1‑69 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGAt > 1:154109/1‑70 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGAt > 1:2261944/1‑70 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc < 1:247980/71‑1 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc < 1:4205993/71‑1 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc < 1:5529276/71‑1 (MQ=255) aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc > 1:2032962/1‑71 (MQ=255) gCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg < 1:4420269/71‑1 (MQ=255) cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg < 1:837594/70‑1 (MQ=255) cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa < 1:4325118/71‑1 (MQ=255) cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa < 1:4825885/71‑1 (MQ=255) aCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca > 1:4505128/1‑71 (MQ=255) gTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACAccc > 1:4966644/1‑71 (MQ=255) ctATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACAccc > 1:3472373/2‑71 (MQ=255) ctATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACAcc > 1:3260613/2‑70 (MQ=255) tATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGa < 1:2438894/61‑1 (MQ=255) tCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTAt > 1:1388478/1‑60 (MQ=255) cAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGtag > 1:557796/1‑69 (MQ=255) aGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGtagc > 1:2117312/1‑68 (MQ=255) aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCg < 1:577294/71‑1 (MQ=255) aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCg < 1:334152/71‑1 (MQ=255) aGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCAGTTCg < 1:2738612/71‑1 (MQ=255) ggATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGac < 1:12983/61‑1 (MQ=255) gATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGCATATcgcccgcc < 1:898825/53‑1 (MQ=255) gATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGt < 1:3629170/71‑1 (MQ=255) tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa < 1:3104137/71‑1 (MQ=255) tCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGa < 1:1917789/71‑1 (MQ=255) aCCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGaa < 1:3637375/70‑1 (MQ=255) cAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaaga < 1:3081619/71‑1 (MQ=255) cAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaaga < 1:120565/71‑1 (MQ=255) cAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaag > 1:3904650/1‑70 (MQ=255) tGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc < 1:2354310/48‑1 (MQ=255) gATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACAcc > 1:4474786/1‑52 (MQ=255) gATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg < 1:4011875/70‑1 (MQ=255) gATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg < 1:929108/70‑1 (MQ=255) aTGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGAt < 1:2580975/71‑1 (MQ=255) aTGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGAt < 1:4214503/71‑1 (MQ=255) gCGTGCTTTTGGTGAAAGCCTGGTGGCTTATcgcccgcc < 1:3254970/39‑1 (MQ=19) cGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:3909307/65‑1 (MQ=255) cGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:2438277/65‑1 (MQ=255) gTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca > 1:1074618/1‑45 (MQ=255) gTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGc < 1:4730299/71‑1 (MQ=255) gTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGc > 1:5599783/1‑71 (MQ=255) tGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg < 1:415935/37‑1 (MQ=255) tGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg < 1:2063339/37‑1 (MQ=255) cttttGGTGAAAGCCTGGTGGCTTATCGCCCGCCTATCGACACCCGta > 1:868645/2‑47 (MQ=255) tttCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGa < 1:2058021/71‑1 (MQ=255) ttCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg > 1:4734083/1‑37 (MQ=255) gTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATcc < 1:3319256/71‑1 (MQ=255) tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa < 1:4452795/56‑1 (MQ=255) tGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa < 1:4493381/56‑1 (MQ=255) aaaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAaagaa > 1:3234401/1‑50 (MQ=255) aaGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:4145517/51‑1 (MQ=255) gCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGc < 1:4618459/71‑1 (MQ=255) ctgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCaa > 1:5376213/1‑65 (MQ=255) ctgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCaa > 1:2972975/1‑65 (MQ=255) ctgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCaa > 1:2953691/1‑65 (MQ=255) tgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:4062973/46‑1 (MQ=255) tgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:39515/46‑1 (MQ=255) tgctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:4365230/46‑1 (MQ=255) gctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc > 1:5526011/1‑70 (MQ=255) gctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc > 1:5248953/1‑70 (MQ=255) ctgGTTTATCTTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa > 1:2626553/1‑46 (MQ=255) ctgGTTTATCTTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa > 1:3998311/1‑46 (MQ=255) ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGaaa < 1:5093938/57‑1 (MQ=255) ctgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAAc < 1:1172881/69‑1 (MQ=255) tgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGt < 1:5152740/43‑1 (MQ=255) tgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa < 1:3383669/45‑1 (MQ=255) tgGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa < 1:2605016/45‑1 (MQ=255) tttATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTg < 1:1546638/41‑1 (MQ=255) ttATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaa > 1:3051587/1‑70 (MQ=255) tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa > 1:4625639/1‑40 (MQ=255) tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGa > 1:4965810/1‑40 (MQ=255) tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa < 1:5466708/70‑1 (MQ=255) tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaaa > 1:2027858/1‑70 (MQ=255) tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaa > 1:1439398/1‑69 (MQ=255) tATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGaa > 1:1148456/1‑69 (MQ=255) tCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAg < 1:668103/71‑1 (MQ=255) | GTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAG > minE/834667‑834798 |
| 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: ‑ |