| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,855 | Δ1 bp | 39.3% | coding (573/1539 nt) | narH → | nitrate reductase 1, beta (Fe‑S) subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 835,855 | 0 | A | . | 39.3% | 149.1 / 139.2 | 89 | coding (573/1539 nt) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base A (25/29); new base . (22/13); total (47/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.36e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01 | |||||||||||
AGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACC > minE/835790‑835922 | aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtgt > 1:1436486/1‑71 (MQ=255) aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc > 1:2746859/1‑70 (MQ=255) gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGa > 1:4997169/1‑71 (MQ=255) gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCa < 1:2105302/61‑1 (MQ=255) gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCa < 1:370947/61‑1 (MQ=255) gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg < 1:4585217/70‑1 (MQ=255) gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg < 1:1012600/70‑1 (MQ=255) gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCg < 1:4811165/70‑1 (MQ=255) ttgAAAACAACTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGAcc > 1:3261578/3‑71 (MQ=255) tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg > 1:2505592/1‑63 (MQ=255) tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGCACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt > 1:2191295/2‑71 (MQ=255) tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt > 1:664389/2‑71 (MQ=255) tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt > 1:727956/2‑71 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgtg > 1:5626662/1‑62 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCAtgt < 1:2833587/61‑1 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt > 1:139118/1‑71 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt > 1:4604916/1‑71 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCt > 1:83831/1‑71 (MQ=255) gAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTTCCTCAATCCCAGC‑TGCGTGGCGACCTg > 1:2362692/1‑71 (MQ=255) gAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTg > 1:4134574/1‑71 (MQ=255) aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGc < 1:2847320/71‑1 (MQ=255) aaaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGc > 1:559860/1‑71 (MQ=255) aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGc > 1:3584047/1‑69 (MQ=255) acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCg < 1:1157043/71‑1 (MQ=255) acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCg < 1:3126183/71‑1 (MQ=255) acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGccca < 1:3879190/71‑2 (MQ=255) acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGccca < 1:3838918/71‑2 (MQ=255) cacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGccc > 1:2062836/1‑69 (MQ=255) acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg < 1:2239291/71‑1 (MQ=255) acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg < 1:1122157/71‑1 (MQ=255) acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg < 1:2812730/71‑1 (MQ=255) acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg < 1:3395560/71‑1 (MQ=255) ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAg < 1:5009832/70‑1 (MQ=255) ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGc > 1:1982968/1‑71 (MQ=255) tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGt > 1:4422061/1‑71 (MQ=255) tCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggc > 1:1641382/1‑70 (MQ=255) cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCggc < 1:4020327/70‑2 (MQ=255) cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCgg < 1:860058/69‑1 (MQ=255) atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCgg < 1:73443/68‑1 (MQ=255) atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgc > 1:3207135/1‑71 (MQ=255) tgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg < 1:5125726/71‑1 (MQ=255) tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgct < 1:5492776/71‑2 (MQ=255) tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCACAGC‑TGCGTGGCGACCTGCCCAAGCGgcgc > 1:3424952/1‑70 (MQ=255) atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt < 1:5454100/71‑1 (MQ=255) atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt < 1:4593942/71‑1 (MQ=255) atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgcta < 1:1889141/70‑3 (MQ=255) cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc > 1:2959457/3‑70 (MQ=255) cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc > 1:965612/3‑70 (MQ=255) tGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAg > 1:4295982/1‑71 (MQ=255) gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac > 1:1914269/1‑69 (MQ=255) gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac > 1:5474769/1‑69 (MQ=255) ccGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg > 1:4131440/1‑71 (MQ=255) ccGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGT‑GCGACCTGCCCGAGCGGTGCGATTTACAAGCGt > 1:4347390/1‑71 (MQ=255) gcgcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGc > 1:3401348/1‑39 (MQ=255) gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa > 1:3909760/1‑66 (MQ=255) gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa > 1:2057523/1‑66 (MQ=255) gcgcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa > 1:4792923/1‑66 (MQ=255) gcCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTga < 1:1424948/70‑1 (MQ=255) gcCTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaa < 1:4008683/71‑1 (MQ=255) cTGTGCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGgcgct < 1:1982515/54‑2 (MQ=255) tgtgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCt > 1:863524/1‑39 (MQ=255) tgtgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaag < 1:5114094/69‑1 (MQ=255) gtgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa > 1:3075944/1‑67 (MQ=255) tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc < 1:2743499/57‑1 (MQ=255) tgCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:5674081/69‑1 (MQ=255) tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCg > 1:2189811/1‑46 (MQ=12) tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa < 1:5420271/69‑1 (MQ=255) tgCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa < 1:1690358/69‑1 (MQ=255) gCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg > 1:2820836/1‑50 (MQ=255) gCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCC‑AGCGGTGCg > 1:1934623/1‑49 (MQ=255) gCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa > 1:298998/1‑58 (MQ=255) gCGAACACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa > 1:4783294/1‑58 (MQ=255) cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg < 1:4920128/60‑1 (MQ=255) cGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg > 1:1109417/1‑62 (MQ=255) gAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga > 1:5406577/1‑68 (MQ=255) aaCACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg > 1:2365420/1‑70 (MQ=255) aaCACTGCCTCAATCCCAGC‑TGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAagaagacg < 1:809713/69‑3 (MQ=255) cTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:4544790/61‑1 (MQ=255) cTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCg < 1:3433784/71‑1 (MQ=255) tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATTGCATCGTCCt > 1:854048/1‑70 (MQ=255) tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc > 1:1163641/1‑65 (MQ=255) tGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt > 1:5067144/1‑70 (MQ=255) aa‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag > 1:3137770/1‑52 (MQ=255) aa‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGGGAAGAAGATGGCATCGTCCTGAt < 1:4859008/71‑1 (MQ=255) cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg < 1:4446315/71‑1 (MQ=255) ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa < 1:998695/51‑1 (MQ=255) cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag < 1:3814880/51‑1 (MQ=255) cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc < 1:2858466/71‑1 (MQ=255) ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag < 1:2387615/50‑1 (MQ=255) ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc < 1:4683953/71‑1 (MQ=255) cATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag < 1:4593058/48‑1 (MQ=255) | AGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAA‑CCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACC > minE/835790‑835922 |
| 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: ‑ |