Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,363,673 | 0 | G | T | 61.5% | ‑2.1 / 15.3 | 13 | L13I (CTT→ATT) | nanK | N‑acetylmannosamine kinase |
Reads supporting (aligned to +/- strand): ref base G (4/1); new base T (0/8); total (4/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.99e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTGGCGTAGGAAGTTCACGACGATCGCGGATCTGCCCGTCAGCGCCAATCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCTGTGTTGTACCACTGACAAATGTGCTC > NZ_CP009273/3363609‑3363756 | gctggcGTAGGAAGTTCACGACGATCGCGGATCTGCCCGTCAGCGCCAATCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGcc < 2:100445/90‑1 (MQ=255) ctggcGTAGGAAGTTCACGACGATCGCGGATCTGCCCGTCAGCGCCAATCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGCCa > 2:370770/1‑90 (MQ=255) gTTCACGACGATCGCGGATCTGCCCGTCAGCGCCAATCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAg > 2:8363/1‑90 (MQ=255) cgaTCGGGGTTCTGCCGGTCAGCGCCAATCGGCGGGGGGGAAATTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCt < 1:107274/90‑1 (MQ=255) gcgGTTCTGCCCGTCAGCGCCAACCGGCGGGGCGGCAATTTTAGTCCCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTc < 1:370770/90‑1 (MQ=255) tCTGCCCGTCAGCCCCAATCAGCGGGGCGGAAATTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGc < 2:270665/90‑1 (MQ=255) gCCCGTCAGCGCCAATCGGCGCGGCGCCATTTTTAGTTCCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCtgt < 1:417307/90‑1 (MQ=255) gCCCGTCAGCGCCAATCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCtgt > 1:232843/1‑90 (MQ=255) ggcgcCATTCTGCGCGGGGGAAATTTTTGCGCCGCCGATATCAATCGCCAGTGGGGTCATAGCTCCGCCTTTTTCATCGCTGTGTTGTAc < 1:82360/89‑1 (MQ=255) cgcAAATCGGCGGGGGGGAAATTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCTGTGTTGTACCa < 2:503113/90‑1 (MQ=255) aaTCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCTGTGTTGTACCACTGa > 2:127053/1‑90 (MQ=255) ggcggGGCGGAAATTTTAGTCCCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCTGTGTTGTACCACTGACaaa < 1:480481/89‑1 (MQ=255) cggcAATTTTAGTCCCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCTGTGTTGTACCACTGACAAATGTGCTc < 1:65315/90‑1 (MQ=255) | GCTGGCGTAGGAAGTTCACGACGATCGCGGATCTGCCCGTCAGCGCCAATCAGCGCGGCGGCAAGTTTAGTACCGCCGATATCAATCGCCAGTGTGGTCATAGCACCGCCTTTTTCATCGCTGTGTTGTACCACTGACAAATGTGCTC > NZ_CP009273/3363609‑3363756 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |