| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,293,531 | G→A | V123V (GTG→GTA) | nrfF → | putative formate‑dependent nitrite reductase complex subunit NrfF |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,293,531 | 0 | G | A | 96.9% | 111.7 / NA | 32 | V123V (GTG→GTA) | nrfF | putative formate‑dependent nitrite reductase complex subunit NrfF |
| Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (11/20); minor base T (1/0); total (12/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTTAACGGGTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTGAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCT > NC_000913/4293453‑4293606 | gTTAACGGGTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAg < 2:856139/81‑1 (MQ=255) aCGGGTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGc < 1:856140/80‑1 (MQ=255) cGGGTCAGACGCTGGTGTTATGGGCGCTGCCGGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGa < 2:856141/81‑1 (MQ=255) gggTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGaa < 2:856142/81‑1 (MQ=255) gTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGc < 2:856144/81‑1 (MQ=255) gTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGc < 2:856143/81‑1 (MQ=255) cGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGa > 1:856178/1‑80 (MQ=255) gCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGaaa > 1:856191/1‑81 (MQ=255) gtgtTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAg < 2:856145/80‑1 (MQ=255) tgtTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGcc > 2:856187/1‑81 (MQ=255) ggCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAAt < 2:856146/80‑1 (MQ=255) ggtgTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATg > 1:856203/1‑81 (MQ=255) ttgttACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTg < 2:856147/80‑1 (MQ=255) ttACTGCTGATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGaca < 1:856148/80‑1 (MQ=255) gctgATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACt > 2:856193/1‑80 (MQ=255) ctgATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTc < 1:856149/80‑1 (MQ=255) tgATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCAc > 1:856197/1‑81 (MQ=255) gATGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCAc > 1:856194/1‑80 (MQ=255) aTGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCAc < 1:856150/79‑1 (MQ=255) tGGCACTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGat < 2:856151/81‑1 (MQ=255) aCTGATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTg < 1:856152/81‑1 (MQ=255) aCTGATCCTCTGGCGACTTAGCGCGCCTCCATGAAACAGCCCAAAATACCGCTCAAAATGCTTACAACACTCACGATATTg > 2:856200/1‑81 (MQ=255) gATCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATg > 2:856188/1‑81 (MQ=255) aTCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATgg < 2:856154/81‑1 (MQ=255) aTCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATg < 2:856153/80‑1 (MQ=255) tCCTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATgg > 1:856192/1‑80 (MQ=255) ccTCTGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTa > 1:856189/1‑81 (MQ=255) ctGGCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTAttt > 2:856190/1‑81 (MQ=255) gCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTAttt < 1:856156/78‑1 (MQ=255) gCGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTtctc < 1:856155/81‑1 (MQ=255) cGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTctct < 2:856157/81‑1 (MQ=255) cGAGTAAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTctct < 2:856158/81‑1 (MQ=255) | GTTAACGGGTCAGACGCTGGTGTTATGGGCGCTGCCAGTGGTGTTGTTACTGCTGATGGCACTGATCCTCTGGCGAGTGAGGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCT > NC_000913/4293453‑4293606 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |