Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,527,926 | 0 | A | C | 61.5% | 2.1 / 12.4 | 13 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/8); total (5/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.77e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATC > NZ_CP009273/4527863‑4528006 | taaaaaTCCGCCACCCTTTCGCAATCCTGTTTT‑GGCTCCCCCATTAATAACCCAGT‑TTTTATCCCT‑TTTATC‑CACAACCGCCGCACCAGc < 1:629098/87‑1 (MQ=255) ccttCCCAAACCTTGTTTGGGCTTC‑CCAATAAAAAACCAGGTTTTTACCCCT‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAgg < 1:576140/88‑1 (MQ=255) aaTCCTGGTTT‑GGCTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg > 2:371458/1‑90 (MQ=255) cccTGTTTT‑GGCTTCCCCATTAAAAACCCAGT‑TTTTATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 1:16796/89‑1 (MQ=255) cTGGTTT‑TGCCTCCCCAATAAAAAACCCGGTTTTTACCCCT‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa < 2:318619/90‑1 (MQ=255) tt‑GGCTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa > 1:146013/1‑90 (MQ=255) ggCTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 2:218693/1‑90 (MQ=255) tcccAATAAAAAACCCGG‑TTTTTCCCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt < 2:52651/88‑1 (MQ=255) ccATTAATAACCCAGG‑TTTTATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 1:293777/90‑1 (MQ=255) caataaaaaGCCCGG‑TTTTTACCCC‑TTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:317110/87‑1 (MQ=255) aTTAAAAACCCA‑G‑TTTTTTCCCC‑TTTTTTCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:66358/90‑1 (MQ=255) gCCAGG‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg > 1:660634/1‑90 (MQ=255) gg‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc > 1:359837/1‑90 (MQ=255) | TTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATC > NZ_CP009273/4527863‑4528006 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |