Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 4,306,692 | A→T | I115N (ATC→AAC) | phnM ← | alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,306,692 | 0 | A | T | 93.3% | 148.8 / 2.5 | 45 | I115N (ATC→AAC) | phnM | alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase |
Reads supporting (aligned to +/- strand): ref base A (2/1); new base T (20/22); total (22/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01 |
GGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT > NZ_CP009273/4306607‑4306777 | ggtgTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGttcat > 2:184554/1‑90 (MQ=255) gtgtGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGatcatc < 1:608217/90‑1 (MQ=255) aTGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTtctc > 2:975564/1‑90 (MQ=255) aTGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTtctc > 2:594863/1‑90 (MQ=255) gCGGCAGTTCGCAGCGCAGATGGAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCa > 1:124150/1‑90 (MQ=255) gCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCa > 1:741780/1‑90 (MQ=255) gCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCa > 2:696511/1‑90 (MQ=255) gCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGAt < 1:1536275/90‑1 (MQ=255) gcgcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCg > 2:1337382/1‑90 (MQ=255) gcgcagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCg > 1:249214/1‑90 (MQ=255) cagATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATc > 1:396641/1‑90 (MQ=255) aTGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATcgcc < 1:392045/90‑1 (MQ=255) gCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATcgccgc < 1:1415764/90‑1 (MQ=255) gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg > 1:719500/1‑90 (MQ=255) gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg > 1:1316493/1‑90 (MQ=255) gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg > 2:900275/1‑90 (MQ=255) gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAg > 1:47599/1‑63 (MQ=255) gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAg < 2:47599/63‑1 (MQ=255) gTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTcgcg > 1:1069849/1‑90 (MQ=255) gCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGttcatc > 2:856874/1‑53 (MQ=255) gCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGttcatc < 1:856874/53‑1 (MQ=255) cTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCAc > 1:195983/1‑90 (MQ=255) cgcgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGc > 1:702238/1‑90 (MQ=255) cgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa > 1:102963/1‑90 (MQ=255) cgGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCa > 1:1230132/1‑90 (MQ=255) cgccgcGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCa > 2:384599/1‑90 (MQ=255) gccgcGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCAc < 1:1206134/90‑1 (MQ=255) gccgcGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCAc < 1:184554/90‑1 (MQ=255) gccgcGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCAc < 2:741780/90‑1 (MQ=255) gccgcGTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCAc < 1:518659/90‑1 (MQ=255) cccgcCTTTCTGCGTCTCTCCGATGGCGTTCTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCAc < 1:695880/89‑1 (MQ=255) cgggTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 1:796028/87‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 2:720419/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 2:619859/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 2:234243/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 2:1230132/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 2:1168922/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 1:696511/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 1:900275/90‑1 (MQ=255) cgcgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACgg < 1:640467/90‑1 (MQ=255) cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa < 2:336652/90‑1 (MQ=255) cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa < 2:615043/90‑1 (MQ=255) cgTTTCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCa > 1:856807/1‑90 (MQ=255) ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc < 2:856807/90‑1 (MQ=255) ttCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc < 1:933221/90‑1 (MQ=255) tCTGCGTCTCTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCa > 2:1281202/1‑90 (MQ=255) tctcTTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACgg > 1:798612/1‑90 (MQ=255) tctTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc > 1:1432200/1‑81 (MQ=255) tctTCGATGGCGTTGTTCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATcc < 2:1432200/81‑1 (MQ=255) gttgttCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCt < 2:553924/90‑1 (MQ=255) gttgttCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCt > 1:978559/1‑90 (MQ=255) | GGTGTGATGCGGCAGTTCGCAGCGCAGATGCAGACGGTGCTCGGCGCGGTTGACGCCGCGTTTCTGCGTCTCTTCGATGGCGTTGATCATCTTCTCCAGATTCTCCAGCCGATCGCCGCCGTCGCGCACGTCGCCAATTGCCACGGCATCCAGTACGGTGGTGATGCCGCT > NZ_CP009273/4306607‑4306777 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |