| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 226,010 | (A)6→5 | coding (357/804 nt) | dkgB → | 2,5‑didehydrogluconate reductase DkgB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 226,005 | 0 | A | . | 96.4% | 112.5 / ‑1.6 | 28 | coding (352/804 nt) | dkgB | 2,5‑didehydrogluconate reductase DkgB |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (16/11); total (17/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
CGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGG > NZ_CP009273/225917‑226084 | cGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCaa < 1:863819/90‑1 (MQ=255) gATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCaaa > 1:328820/1‑90 (MQ=255) aTGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCaaaaaca < 1:77223/90‑3 (MQ=255) aTCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGc > 2:678812/1‑90 (MQ=255) aTCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCt > 2:628597/1‑90 (MQ=255) tCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTg < 1:628597/90‑1 (MQ=255) gACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGAcg > 2:441611/1‑90 (MQ=255) tAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTga > 1:668703/1‑90 (MQ=255) cTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGt < 1:800582/90‑1 (MQ=255) ggCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTAt > 2:465201/1‑90 (MQ=255) ggCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTAt > 1:185163/1‑90 (MQ=255) ggCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTCAGATCGGTAt > 2:618334/1‑90 (MQ=255) ccGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTAttt > 1:479693/1‑90 (MQ=255) gTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTcc > 2:659041/1‑90 (MQ=255) tCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCa < 1:879700/90‑1 (MQ=255) aCGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCAc < 1:733071/90‑1 (MQ=255) aCGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCAc < 2:668703/90‑1 (MQ=255) aCGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCAc < 2:723505/90‑1 (MQ=255) tctGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGa > 1:865080/1‑90 (MQ=255) gAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCa < 1:357909/66‑1 (MQ=255) gAAGAGTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCa > 2:357909/1‑66 (MQ=255) agagTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGaaaaa < 2:735513/90‑1 (MQ=255) gTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGt < 2:378565/76‑1 (MQ=255) gTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGt > 1:378565/1‑76 (MQ=255) gTTTATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGGGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCg > 1:424893/1‑90 (MQ=255) ttATGCAGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGAt > 2:423861/1‑90 (MQ=255) aGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGAt > 1:66956/1‑42 (MQ=39) aGGCGCTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGAt < 2:66956/42‑1 (MQ=39) gcgcTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGctgctg > 1:754222/1‑90 (MQ=255) gcgcTGCTGGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGctgctg > 1:604233/1‑90 (MQ=255) tgctgGAAGCC‑AAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTgg < 1:332207/90‑1 (MQ=255) | CGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTGACGCGTGAGATCGGTATTTCCAACTTCACGATCCCGTTGATGGAAAAAGCGATTGCTGCTGTTGG > NZ_CP009273/225917‑226084 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |