| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,527,926 | 0 | A | C | 56.2% | 10.6 / 19.5 | 16 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
| Reads supporting (aligned to +/- strand): ref base A (7/0); new base C (0/9); total (7/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.74e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGG > NZ_CP009273/4527853‑4528003 | tGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAAccg > 1:97580/1‑90 (MQ=255) aaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAAccgccg > 1:304909/1‑90 (MQ=255) atccGGCACCCGTTCGAAATCCTGTTTT‑GGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTc < 2:276291/87‑1 (MQ=255) tccGGCATCCGTCCCAAACCCTGTTTT‑GGCTTCCCCATT‑AATAAGCCAGT‑TTTTATCCCCTTTATT‑CACAACCGCCGCACCAGCCGTTcc < 2:183043/88‑1 (MQ=255) cGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGt > 1:762451/1‑90 (MQ=255) ccctgTTT‑GGCCTCCCCATT‑AAAAAGCCCGG‑TTTTTACCCCTTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 1:28370/86‑1 (MQ=255) ccTTGTTTGGGCTTCCCCTTT‑AATACGCCGGT‑TTTTATCCCTTTTTTC‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 2:771608/90‑1 (MQ=255) ttgggCTTC‑CCATTAAAAAAGCCAGT‑TTTTACCCCTTTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 2:304909/87‑1 (MQ=255) ggCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 2:151405/1‑90 (MQ=255) ggCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 2:341579/1‑90 (MQ=255) ggCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 2:421822/1‑90 (MQ=255) cATT‑AAAAAACCCGGTTTTTACCCCTTTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 1:315813/90‑1 (MQ=255) ataaataaCCCGGT‑TTTTACCCCTTTTATT‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:209813/87‑1 (MQ=255) at‑aaaaaGCCGGGTTTTTATCCCTTTTTTC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 1:804186/89‑1 (MQ=255) gCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg > 2:586318/1‑90 (MQ=255) cAGG‑TTTTTCCCCCTTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg < 1:680901/90‑1 (MQ=255) | TGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGG > NZ_CP009273/4527853‑4528003 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |