| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 837,728 | C→G | 7.8% | F67L (TTC→TTG) | narI → | nitrate reductase 1, gamma (cytochrome b(NR)) subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 837,728 | 0 | C | G | 7.8% | 154.0 / 3.4 | 51 | F67L (TTC→TTG) | narI | nitrate reductase 1, gamma (cytochrome b(NR)) subunit |
| Reads supporting (aligned to +/- strand): ref base C (23/24); new base G (2/2); total (25/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01 | |||||||||||
CGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAA > minE/837660‑837798 | cGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtctt > 1:1498950/1‑71 (MQ=255) aaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGt < 1:3770314/71‑1 (MQ=255) aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt < 1:709286/71‑1 (MQ=255) aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt < 1:1074279/71‑1 (MQ=255) aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt > 1:1116994/1‑71 (MQ=255) aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt < 1:703249/71‑1 (MQ=255) gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATg > 1:3353311/1‑67 (MQ=255) aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa < 1:2689415/70‑1 (MQ=255) aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa > 1:1468098/1‑70 (MQ=255) aCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcg > 1:3558725/1‑71 (MQ=255) ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa > 1:473151/1‑68 (MQ=255) tGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTa < 1:767103/60‑1 (MQ=255) ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg > 1:2619216/1‑70 (MQ=255) tggTCGAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAAcgccgc < 1:3267859/68‑1 (MQ=255) tCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg < 1:4073534/71‑1 (MQ=255) gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGa > 1:1150103/1‑71 (MQ=255) aaCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccgc < 1:712000/64‑1 (MQ=255) aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGAt > 1:3365325/1‑70 (MQ=255) aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGAt > 1:3345214/1‑70 (MQ=255) aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTAGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATg > 1:3738402/1‑71 (MQ=255) tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt > 1:2807950/1‑71 (MQ=255) gTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTg > 1:4077518/1‑52 (MQ=255) ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg < 1:1268234/44‑1 (MQ=255) ttCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATGt < 1:2914066/67‑1 (MQ=255) atatCGGGATTCTGGGGATTTTTGTCGGTCACttcttc < 1:446137/38‑1 (MQ=255) atatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAg < 1:2977024/69‑1 (MQ=255) atatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGcc < 1:2957944/71‑1 (MQ=255) tatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccgc < 1:744199/53‑1 (MQ=255) tCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGc < 1:2446692/67‑1 (MQ=255) tCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg < 1:402644/70‑1 (MQ=255) cGGAATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGc > 1:618234/1‑71 (MQ=255) gggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAc > 1:1741810/1‑44 (MQ=255) gggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt > 1:2940257/1‑71 (MQ=255) ggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg < 1:198432/67‑1 (MQ=255) ttGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGa > 1:252109/2‑49 (MQ=16) ttGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGa > 1:1597727/2‑49 (MQ=16) cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa > 1:3290693/1‑37 (MQ=255) tGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGt < 1:2581900/51‑1 (MQ=255) gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTg < 1:352018/70‑1 (MQ=255) gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa < 1:3092479/71‑1 (MQ=255) gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa < 1:1051889/71‑1 (MQ=255) tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAtt > 1:135447/1‑65 (MQ=255) ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa > 1:217656/1‑71 (MQ=255) ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg > 1:3533683/1‑63 (MQ=255) ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGa > 1:1408935/1‑64 (MQ=255) ggTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGa > 1:935723/1‑70 (MQ=255) gTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaa < 1:1068878/71‑1 (MQ=255) tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa > 1:1883278/1‑71 (MQ=255) tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa < 1:966850/71‑1 (MQ=255) aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa > 1:3648140/1‑51 (MQ=255) aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt < 1:3360085/70‑1 (MQ=255) aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt < 1:2728200/70‑1 (MQ=255) aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATg < 1:133515/71‑1 (MQ=255) cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGaa < 1:752778/52‑1 (MQ=255) cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaaa > 1:3100994/1‑65 (MQ=255) cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa > 1:291927/1‑71 (MQ=255) | CGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAA > minE/837660‑837798 |
| 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: ‑ |