Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,370,630 | 0 | A | C | 57.1% | 14.1 / 21.8 | 21 | Y148D (TAT→GAT) | menH | 2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase |
Reads supporting (aligned to +/- strand): ref base A (4/5); new base C (0/12); total (4/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.11e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGC > NZ_CP009273/2370544‑2370708 | aaGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCa > 2:1218171/1‑90 (MQ=255) gCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGt < 1:1521368/90‑1 (MQ=255) attattGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATAc > 2:538189/1‑90 (MQ=255) tttgttgCGGCGCAGCGCCCCCAGCCCCCGGGGTTGATCGTCATTGAGTGAGGAAAAAACGGGCTGTTGATCCCAGTCGGAAAATAccgc < 1:890040/89‑1 (MQ=255) tttgttgCGGCGCAGCCCCACCAGCCCCCGGCTTTGATCGTCATGGAGTGAGGAAAAAACAGGCTGTTGATCCCAGTCGGCAAATAccgc < 1:1498893/89‑1 (MQ=255) cggcgcAGCCCCACCAGCCCCCGGCGTTGCTCTTCATTGAGTGAGGAAAAAACGGGCTGTTGATCCCAGTCGGCAAATACCGCCGTTAAc < 1:1216526/88‑1 (MQ=255) aGCCCCCCCAGCCCCCGGCGTTGATCGTCATTGAGTGAGGAAAAAACAGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCt < 1:538189/90‑1 (MQ=255) aGCCCCCCCAGCCCCCGGCGTTGATCGTCATTGAGTGAGGAAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCt < 2:1437770/90‑1 (MQ=255) ccaGCCCCGGGCGTTGATCGTAATTGGGGGAGGAAAAAACAGGCTGTTGACCCCAGTCGGAAAATCCCGCCGTTAACGGTTCTGTGAGaa < 1:176585/90‑1 (MQ=255) ccaGCCCCCGGCGTTGATCGTCTTTGAGTGAGAAAAAAACGGGCTTTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGaa < 1:846117/90‑1 (MQ=255) ccGGGGTTGATCGTCATTGAGGGAGGAAAAAACAGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTg < 1:606689/90‑1 (MQ=255) ccGGCGTTGATCGTCATTGAGTGAGGAAAAAACAGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTg < 1:424504/90‑1 (MQ=255) cGGGGTTTACCGTCATGGGGTGAGGCAAAAAGGGGTTTTTGAAACCGGTGGAAAAACACCGCCGTTAACGGTTCTGTGAGAAAACGCTGc < 2:905662/90‑1 (MQ=255) ggCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCt < 2:1306637/87‑1 (MQ=255) ggCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCt > 1:1306637/1‑87 (MQ=255) gCGTTGATCGTCATTGAGTGAGGAAAAAACGGGCTGTTGATCCCGGCCGCCAAATCCCGCCGTTAACGGTTCTGTGAGAAAACGCTGCAc < 1:1299308/90‑1 (MQ=255) tCGTCATTGGGTGAGGCAAAAACAGGCTGTTGATCCCAGTCGGAAAATCCCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTgg < 1:597201/90‑1 (MQ=255) gggagGAAAAAACAGGCTGTTGACCCCAGGCGCAAAATCCCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAAc < 1:1367431/88‑1 (MQ=255) agGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGc > 1:817079/1‑68 (MQ=255) agGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGc < 2:817079/68‑1 (MQ=255) aGGGTGTTGATCCCAGTCGGAAAATCCCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGc < 1:1024812/90‑1 (MQ=255) | AAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGC > NZ_CP009273/2370544‑2370708 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |