Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,524,293 | A→G | intergenic (+98/‑105) | BW25113_RS07645 → / → BW25113_RS26145 | RHS repeat protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,524,293 | 0 | A | G | 96.7% | 103.3 / ‑3.5 | 30 | intergenic (+98/‑105) | BW25113_RS07645/BW25113_RS26145 | RHS repeat protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (13/16); total (13/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01 |
AAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCA > NZ_CP009273/1524157‑1524420 | aaagcaaagcAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATgaa < 2:850322/139‑1 (MQ=255) aagcaaagcaAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATgaaa < 1:269157/139‑1 (MQ=255) aagcaaagcaAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATgaaa < 2:797734/139‑1 (MQ=255) agcaaGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCaa < 2:321835/139‑1 (MQ=255) gcaaGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCaaa < 2:870778/139‑1 (MQ=255) ggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCg < 1:537033/139‑1 (MQ=255) ttCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATgaa < 2:750301/122‑1 (MQ=255) ttCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATgaa > 1:750301/1‑122 (MQ=255) cAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGtt > 2:1193039/1‑135 (MQ=255) aaTTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCt < 2:590017/139‑1 (MQ=255) gCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAatat > 2:326245/1‑139 (MQ=255) tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa > 2:897752/1‑139 (MQ=255) gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg > 1:722103/1‑138 (MQ=255) cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACACGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTc < 1:359691/139‑1 (MQ=255) tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAaca > 2:1000626/1‑52 (MQ=255) tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAaca < 1:1000626/52‑1 (MQ=255) tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt > 1:1151874/1‑139 (MQ=255) gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGt < 1:761044/62‑1 (MQ=255) gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGt > 2:761044/1‑62 (MQ=255) gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGt > 1:285525/1‑60 (MQ=255) gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGt < 2:285525/60‑1 (MQ=255) gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt > 2:804045/1‑138 (MQ=255) tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat < 1:306479/138‑1 (MQ=255) tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat > 1:293500/1‑138 (MQ=255) tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat < 2:293500/138‑1 (MQ=255) tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat < 2:1151874/138‑1 (MQ=255) aaTTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa < 1:897752/139‑1 (MQ=255) aTTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt > 1:252855/1‑139 (MQ=255) aTTCAAATTCTGACATCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt > 1:159388/1‑139 (MQ=255) aTCGTATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa < 1:326245/139‑1 (MQ=255) tATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc < 2:642626/139‑1 (MQ=255) tATTTTATGGACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc < 1:804045/139‑1 (MQ=255) ttATGGACATGAAACACCAAAGATCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCa < 2:185778/138‑1 (MQ=255) ggACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat > 1:287801/1‑120 (MQ=255) ggACATGAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat < 2:287801/120‑1 (MQ=255) | AAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCA > NZ_CP009273/1524157‑1524420 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |