Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,034,572 | T→G | Q74H (CAA→CAC) | mobB ← | molybdopterin‑guanine dinucleotide biosynthesis protein MobB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,034,572 | 0 | T | G | 97.1% | 109.9 / ‑3.6 | 35 | Q74H (CAA→CAC) | mobB | molybdopterin‑guanine dinucleotide biosynthesis protein MobB |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (13/21); total (13/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01 |
TCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATG > NZ_CP009273/4034486‑4034658 | tCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgc < 1:31191/90‑1 (MQ=255) cAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctg > 2:392850/1‑90 (MQ=255) aGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctgg > 2:450400/1‑90 (MQ=255) tttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctggctg < 2:302220/90‑1 (MQ=255) tttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctggctg < 1:392850/90‑1 (MQ=255) tttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgctggctg < 1:75133/90‑1 (MQ=255) tGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCtggctggc < 2:22351/90‑1 (MQ=255) aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgc > 2:542988/1‑76 (MQ=255) aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTgc < 1:542988/76‑1 (MQ=255) aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt > 1:357196/1‑90 (MQ=255) aTCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGAt > 1:554533/1‑90 (MQ=255) aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGtt < 1:116872/90‑1 (MQ=255) ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGttt > 2:580347/1‑90 (MQ=255) tGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCcgc > 2:146651/1‑90 (MQ=255) gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc < 1:450400/90‑1 (MQ=255) aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCt < 1:71238/90‑1 (MQ=255) aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCt < 1:151494/90‑1 (MQ=255) aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc < 1:494710/90‑1 (MQ=255) aGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc < 2:337598/90‑1 (MQ=255) gATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCa < 1:146651/90‑1 (MQ=255) aGCTCTTCTTCGTCTGGGGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCa > 1:192485/1‑90 (MQ=255) ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATc < 2:13927/90‑1 (MQ=255) ttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc > 2:90508/1‑64 (MQ=255) ttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc < 1:90508/64‑1 (MQ=255) ttCCGTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt < 1:305665/90‑1 (MQ=255) tCCGTCATAAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTa > 2:45595/1‑90 (MQ=255) gTCATCAAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCa < 2:203705/90‑1 (MQ=255) tcatcaAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCaa < 2:420518/90‑1 (MQ=255) caAGGCCCATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc < 1:280334/90‑1 (MQ=255) cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCa > 2:245740/1‑90 (MQ=255) cATCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCa > 2:498199/1‑90 (MQ=255) aTCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCa < 2:20613/77‑1 (MQ=255) aTCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCa > 1:20613/1‑77 (MQ=255) tCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg < 1:226612/90‑1 (MQ=255) tCGGTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg < 2:59109/90‑1 (MQ=255) | TCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATG > NZ_CP009273/4034486‑4034658 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |