| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,941,604 | 0 | A | C | 62.5% | 16.4 / 12.9 | 16 | G10G (GGT→GGG) | ruvC | crossover junction endodeoxyribonuclease RuvC |
| Reads supporting (aligned to +/- strand): ref base A (4/2); new base C (5/5); total (9/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC > NZ_CP009273/1941517‑1941686 | atcCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAAcc > 2:331960/1‑90 (MQ=255) tGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCaa < 2:34104/90‑1 (MQ=255) gCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGcc > 1:314203/1‑90 (MQ=255) cGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCg < 2:314203/90‑1 (MQ=255) gATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCgag < 2:276261/90‑1 (MQ=255) aTCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGaataa > 2:109529/1‑90 (MQ=255) tCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGaataat < 1:124916/90‑1 (MQ=255) ggTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCAcgc > 1:313321/1‑90 (MQ=255) gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt > 2:96855/1‑90 (MQ=255) gTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGtt > 2:205698/1‑90 (MQ=255) acctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGt < 2:106449/90‑1 (MQ=255) tacctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAg < 1:184674/62‑1 (MQ=255) tacctGGCGGATGACGCCGTAGCCGGTCACGCGCGACCCCGGATCAATGCCGAGAATAATAg > 2:184674/1‑62 (MQ=255) tacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTa > 2:80604/1‑90 (MQ=255) ggTCACGCGCGACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTg > 1:173359/1‑90 (MQ=255) cgcgcgACCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGAcc < 1:105827/90‑1 (MQ=255) | ATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACC > NZ_CP009273/1941517‑1941686 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |