Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 832,042 | 0 | C | G | 63.2% | ‑0.6 / 17.0 | 19 | A79G (GCG→GGG) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
Reads supporting (aligned to +/- strand): ref base C (2/5); new base G (12/0); total (14/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC > NZ_CP009273/831962‑832130 | gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:400016/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGTccc > 2:404241/1‑88 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAc > 2:420770/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGATgggg > 2:802558/1‑88 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCg > 2:491683/1‑90 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 1:40870/90‑1 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCAtt < 1:514954/90‑1 (MQ=255) ggAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCGGTGGCGGTTGCGCCtttt > 2:717615/1‑89 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTc < 1:400016/90‑1 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACGCGCGACGGCGAGACCGCCTTTGGtgaa > 2:514954/1‑87 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACGCCCCCTGGCGGGCGCGCCTTTTGTTAg > 1:814999/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGCTGGAGGACGCCCACTTGGttcg > 2:749450/1‑87 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGATGGGGGTCCCGCATTTTGTTAg > 1:9024/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCGATGGGGGTCGCGCATTTTGTCAg > 1:335844/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCGATGGCGAGCGCGCATTTTGTCAg > 2:627382/1‑90 (MQ=255) atGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCCCAGGCGGTCGCGCATTTTTTGAGGGCGCGGGCCATGAAACGGTGgggc > 2:744997/1‑90 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt < 1:532524/90‑1 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCCCGGTGGCGAACGCGCCTTTGGGTAGGGGGCGGGACATGAAGCGATGGCGCTGGGGGGTgag > 1:687435/1‑90 (MQ=255) agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCATGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc < 1:494917/90‑1 (MQ=255) gcggcgGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc < 1:491683/90‑1 (MQ=255) | GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC > NZ_CP009273/831962‑832130 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |