| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,475,322 | 0 | T | G | 53.3% | 15.4 / 17.0 | 15 | E149D (GAA→GAC) | pepA | leucyl aminopeptidase |
| Reads supporting (aligned to +/- strand): ref base T (4/3); new base G (8/0); total (12/3) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.69e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACG > NZ_CP009273/4475241‑4475411 | gTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTGAAAACCATCTTACGCAGCGGGCGACGCGGGTCGGTCtt > 1:241264/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGGTCGCTCTTGTTCGttttt > 1:80461/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGGTCGCTCTTGTTCGTTTTc > 1:65005/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGGGCGTTTAACACCCCCTTTCGCCGCGGGCGACGCGGGTCGCTCTTGTTCGttttt > 2:355345/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGGACGTTGAACACCCCCTTTCGCCGCGGGCGACGCGGGTCGCTCTTGTTTGttttt > 2:425408/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCCTCTTACGCCGCGGGCGACGCGGGTTGCCCTTGTTTGttttt > 2:5276/1‑89 (MQ=255) gTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGGTCGCTCTTGTTTGTTTTTAGCTGATCGAAACTATaga > 2:240451/1‑90 (MQ=255) gggcGGCACGTTTAACACCATCTTACGCCGCGGGCGACGCGGGTCGCTCTTTTTTGTTTTTCGCTGATCGAAACTGTAGAGGGGCTCttt > 2:199643/4‑90 (MQ=255) acacCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAg > 2:393827/1‑90 (MQ=255) aCGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCAc < 2:168214/90‑1 (MQ=255) aCGCAGCGGACGACGCGGTTCGCTCTTGTCCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCAc < 1:116489/90‑1 (MQ=255) gCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTcc > 2:92565/1‑90 (MQ=255) gCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTcc > 2:345362/1‑90 (MQ=255) acgacgCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCagta < 1:23720/90‑1 (MQ=255) cgGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAAAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAgttgtt > 1:483147/1‑90 (MQ=255) ttCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACg < 1:355345/90‑1 (MQ=255) ttCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACg < 1:318719/90‑1 (MQ=255) | GTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTCGTTTTCAGCTGATCGAAACTGTAGAGCGTCTCTTTTGCCGTCTCGACAGCCTGACGCACTTTCCAGTAGTTGTTACG > NZ_CP009273/4475241‑4475411 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |