| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,370,630 | 0 | A | C | 60.0% | 2.8 / 11.9 | 20 | Y148D (TAT→GAT) | menH | 2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase |
| Reads supporting (aligned to +/- strand): ref base A (4/4); new base C (0/12); total (4/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.44e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-02 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGT > NZ_CP009273/2370548‑2370712 | gTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCg > 1:951825/1‑90 (MQ=255) tattGTTGCTGCGCAGCCCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATAccg < 2:933643/90‑1 (MQ=255) tgctgcGCAGCGCCCCCAGCCCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACGGGCTGTTGATCCCAGCCGGCAAATACCGCCGTTa < 1:1414412/90‑1 (MQ=255) gcgcCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCtg > 2:59547/1‑90 (MQ=255) ccaGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGaa < 2:546304/90‑1 (MQ=255) ccaGCCCCCGGCGTTGATCGTCATTGAGTGAGGAAAAAACAGGCTTTTGATCCCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGaa < 1:409023/90‑1 (MQ=255) caGCCCCCGGCGTTGATCGTCATTGAGTGAGGAAAAAACAGGTTGTAAATCCCAGGCGACAAATCCCGCCGTTAACGGTTCTGTGAGaaa < 1:1024854/90‑1 (MQ=255) cGGCGTTGATCGTCATTGAGTGAGGAAAAAACAGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGc < 1:291045/90‑1 (MQ=255) cGGCGTTGATCGTCATTGAGGGAGGAAAAAACGGGCTGTTGATCCCGGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGc < 1:202300/90‑1 (MQ=255) gTCATGGAGGGAGGCAAAAACGGGTTTTTGATCCCAGTCGGAAAATCCCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCg < 1:1100377/90‑1 (MQ=255) aTTGAGTGAGGAAAAAACAGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATc < 2:101148/90‑1 (MQ=255) gtgagGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAAc < 2:388947/90‑1 (MQ=255) gagGAAAAAACGGGCTGTTGATCCCGGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGc < 1:894089/90‑1 (MQ=255) agGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCt > 1:627939/1‑90 (MQ=255) gaaaaaaCAGGTTGTTGATCCCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGa < 2:985787/88‑1 (MQ=255) gaaaaaaCAGGCTGTTGATCCCAGTCGGAAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCGCCCATTGGCGATCGGAACGCTGa < 2:252918/88‑1 (MQ=255) aaaaaCAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACg > 2:1027607/1‑90 (MQ=255) gCTGTTGATCCCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCAc < 1:1027607/90‑1 (MQ=255) gCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCAc < 1:59547/90‑1 (MQ=255) tgtgGATCCCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACgt < 1:1577244/90‑1 (MQ=255) | GTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCCGGCGTTGATCGTCATTGAGTGAGGCAAAAACAGGCTGTTGATACCAGTCGGCAAATACCGCCGTTAACGGTTCTGTGAGAAAACGCTGCACCCATTGGCGATCGGAACGCTGACGTTCCGCACGT > NZ_CP009273/2370548‑2370712 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |