Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 832,048 | 0 | T | G | 61.9% | 3.9 / 19.3 | 21 | V81G (GTC→GGC) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
Reads supporting (aligned to +/- strand): ref base T (2/6); new base G (13/0); total (15/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.16e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGT > NZ_CP009273/831962‑832133 | 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) gcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGt > 2:814225/1‑90 (MQ=255) | GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGT > NZ_CP009273/831962‑832133 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |