| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 726,200 | 0 | A | C | 53.8% | 11.0 / 16.4 | 13 | T388P (ACC→CCC) | rhsC | rhs element protein RhsC |
| Reads supporting (aligned to +/- strand): ref base A (4/2); new base C (0/7); total (4/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.10e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAA > NZ_CP009273/726122‑726280 | gACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGAc < 2:111556/90‑1 (MQ=40) aGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGc > 2:503980/1‑90 (MQ=38) gAACAGCAAACCCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTc < 2:98632/90‑1 (MQ=14) aCAGCAAACCCCGGCGGGTTTAAGTTACCCGTATCAGTAGGAGAAGGCCCCCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCt < 1:143781/90‑1 (MQ=14) aGCAAACCCGGGCGGGTTAAATTTACCCGTATCATTAGGAGAGGGCCCCCACCCCCCCCCCCGACAGCCTAAACCGCCGTGAAGTCCTGc < 1:309801/90‑1 (MQ=255) cccgggCGGGCTTAAGTTACACGTATCAGTAGGAGAAGGCCCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGc < 1:610714/86‑1 (MQ=14) ccccGGCGGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCACCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGc < 1:63543/89‑1 (MQ=14) ccGGCAGGCTTAATTTACACGTATCAGTATGAGAAGGCCCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAg < 2:192985/90‑1 (MQ=17) gCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcac < 1:150451/44‑1 (MQ=255) gCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcac > 2:150451/1‑44 (MQ=255) gCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGt > 1:192985/1‑90 (MQ=255) aTCAGTAGGAAAAGGCCCGCCTCCCCCTCCCCGACAGCCTGACCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCggg < 2:345345/90‑1 (MQ=17) gACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAgg > 1:43985/1‑62 (MQ=255) gACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAgg < 2:43985/62‑1 (MQ=255) gACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGaa > 2:290879/1‑90 (MQ=255) | GACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAA > NZ_CP009273/726122‑726280 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |