| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 832,048 | 0 | T | G | 52.5% | 43.5 / 56.4 | 40 | V81G (GTC→GGC) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
| Reads supporting (aligned to +/- strand): ref base T (3/16); new base G (21/0); total (24/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831962‑832131 | gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:843719/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:809576/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:76258/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:369105/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:704717/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAc > 1:309690/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAc > 1:732038/1‑90 (MQ=255) aTCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACg < 2:424685/90‑1 (MQ=255) cAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCg > 1:345634/1‑90 (MQ=255) aaGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACGCCCGa > 1:776504/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCCCGATGGCg > 1:177475/1‑90 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 2:627708/90‑1 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 2:732038/90‑1 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCATGGGGGTCgcg‑c > 1:626751/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCTCGATGGCGATCgcg‑c > 1:15477/1‑90 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTc < 1:588108/90‑1 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTc < 2:809576/90‑1 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCCGTGGGGGGCGCG‑CATTTTGGCAg > 1:826663/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCTCGATGGCGATCGCG‑CATTTTGTCAg > 1:322930/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCCCCCCCGATGGGGATCGCG‑CCTTTGGTCAg > 1:941859/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGATCGCG‑CATTTGGGCAg > 1:293170/1‑90 (MQ=255) gggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTc < 1:358736/90‑1 (MQ=255) cACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCGATGGGGATCCCG‑CCTTTGGGCAGGGcgcg > 1:552617/1‑90 (MQ=255) cACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTcgcg < 1:295392/90‑1 (MQ=255) aCCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCCCGATGGCGATCGCG‑CCTTTTGTCAGGGCGCgg > 2:204727/1‑90 (MQ=255) cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACAt < 2:909526/90‑1 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGCCGCCCCCTGGCGGACGCG‑CCTTTTGTTAGGGCGCGGGACAAGAAGCg > 2:700985/1‑90 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCCCCCGGGCGCACGCGCCATTTGTTC‑GGGCGCGGCACCAGAAggg > 2:761245/1‑88 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:700985/90‑1 (MQ=255) gtGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGgcgc < 1:810137/89‑1 (MQ=255) gTTAAAGAGGCGGGGGGGGGGGGCCCGCCCCATGGCGATCGCG‑CCTTTTGTCAGGGGGCGGGCCATGAAACGGTGGGGGTGGGGGTTgac > 2:810137/1‑89 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGGCGCCCCCCGTGGCGAACGCG‑CCTTTTGTCCGGGCGCGGGCCATGAAGCGATGGGGCCGGGGGTTgag > 2:548239/1‑90 (MQ=255) ttAAAGAGGCGGGGGCGGCGGGCACGCCCCATGGGGAGCGCG‑CCTTTTGGGAGGGCGCGGGGCACGGAGCGAGGGGGCTGGGGGGTgaga > 2:810865/1‑90 (MQ=255) ggggCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAg < 2:776504/90‑1 (MQ=255) gcggcgGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc < 1:612588/90‑1 (MQ=255) gcggcgGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc < 2:193025/90‑1 (MQ=255) cggcggTCACGTTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 2:365506/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 2:293170/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:76258/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg > 2:586797/1‑90 (MQ=255) | GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCG‑CATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831962‑832131 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |