Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 832,042 | 0 | C | G | 45.6% | 51.2 / 28.3 | 57 | A79G (GCG→GGG) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
Reads supporting (aligned to +/- strand): ref base C (11/20); new base G (26/0); total (37/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.70e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831954‑832131 | aTTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcg < 2:525269/90‑1 (MQ=255) tGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGgggggc > 1:325934/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGgggggggggggccc > 1:1982767/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:188128/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:397112/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 1:903785/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCAc > 1:866788/1‑90 (MQ=255) aaGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCGa > 2:818774/1‑90 (MQ=255) tATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGGTggg > 1:1276781/1‑89 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACGCCCCATgggg > 2:583261/1‑88 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCCCCCCATGGCg > 2:168786/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGATgggg > 1:330609/1‑88 (MQ=255) aCGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCGCCCCGTGGCggacg > 2:329808/1‑86 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 1:1571325/90‑1 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACGCTCGCTGGCGGTcgggg > 2:1228642/1‑87 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACGCCCCATGGCGGTCgcgc > 2:1571325/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCgcgc > 1:150962/1‑90 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCTCGATGGGGGTCCCGCGtt > 1:2066451/1‑90 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCAtt < 1:97126/90‑1 (MQ=255) aGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCATTTTGt > 1:817244/1‑90 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCATTTTGTc > 1:116396/1‑90 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGTGGTCACGCTCGATGGCGATCGCGCATTTGGTc < 1:397112/90‑1 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATc < 2:783139/78‑1 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATc > 1:783139/1‑78 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCGCCCGATGGGGGTCGCGCATTTTGGCAg > 1:2031998/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCCCGCCGTGGGGAGCGGGCATTTGGTcgg > 2:2210202/1‑88 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACTCTCGGTGGCGGACGCGCCTTTGGTCAg > 2:320361/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCCATGGCGGTCGCGCCTTTGGTCAg > 2:93703/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCCCGCCCCGTGGGGGTCGCGCCTTTGGTCAg > 2:205824/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGGCACCCCCCATGGCGAGCGCGCCTTTTGTCAg > 1:947589/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCAtt < 2:550012/83‑1 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCAtt > 1:550012/1‑83 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGCGCATTTGGTCAg > 1:639401/1‑90 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCCCGCTCGATGGCGATCGCGCATTTTGTCAGGt > 1:295572/1‑90 (MQ=255) gCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCCGTTGCGGTCGCGCATTTTGTCAGGgcgc > 1:2021818/1‑90 (MQ=255) gCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCCCGCCCGATGGCGATCGCGCATTTTGTCAGGgcgc > 1:1463308/1‑90 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGc < 2:434892/90‑1 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGc < 1:19843/90‑1 (MQ=255) cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACAt < 2:2204282/90‑1 (MQ=255) aTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGaa < 2:866788/90‑1 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGggcggc > 1:720651/1‑38 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGggcggc < 2:720651/38‑1 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACCCTCGTGGGGGGTCGCGCCTTTTGGCAGGCCGCGGCACAAGAAGCg > 1:1235226/1‑90 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCa < 2:2186062/43‑1 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCa > 1:2186062/1‑43 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:2194483/90‑1 (MQ=255) catcatGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCCCGATGGCGAACGCGCCTTTGGTCAGGGCGCGGGACAAGAAGGGaag > 1:1764410/1‑88 (MQ=255) tcatGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGAAGGGGATCGCGCAATTTGGCAGGTCGGGGCACATGAAAGGATggg > 1:1189994/1‑89 (MQ=255) tcatGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCCATGGGGGTCGCGCCTTTGGTCAGGGCGCGGGACATGAAGCGATggg > 2:596301/1‑89 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt < 2:1866935/90‑1 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt < 1:93703/90‑1 (MQ=255) ccAAAACCGTTAAAGAGGCGGGGGCGGAGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTgg < 2:1093402/90‑1 (MQ=255) aaaaCCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTgggg < 1:58758/90‑1 (MQ=255) aaCCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGATGGGGCTgggggt > 1:478766/1‑88 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCCCGCTGGCGATCGCGCCTTTGGTCAGGGCGCGGCACATGAAGCGATGGGGCTGGGGGTTgag > 2:18570/1‑90 (MQ=255) agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc < 2:2031998/90‑1 (MQ=255) ggCGGGGGGGGCGGGCACGCCCGATGGCGATCGCGCCTTTTGTTAGGGCGCGGCACATGAAGCGAAGGGGCTGGGGGGTGGGAAAACGCa > 2:885561/1‑90 (MQ=255) cgggggcgggggGCACGCTCCATGGCGGTCGCGCATTTGGTCAGGGCGCGGCACATGAAGCGATGGGGCTGGGGGGTGGGAAAAGGCATc > 2:513951/1‑90 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 1:888227/90‑1 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 1:596301/90‑1 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 2:407123/90‑1 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAGAGCGCATCAGc < 2:1894856/90‑1 (MQ=255) gggCGGCGGGCACGCTCGATGGCGATCGCGCATTTTGTCAGGGCGCGGCCCATGAAGCGATGGGGCTGGGGGTTGAGAAAAGGCATCAGc > 2:888227/1‑90 (MQ=255) ggcggcGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCa < 1:1542892/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:320361/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 2:920581/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:1091552/90‑1 (MQ=255) | ATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831954‑832131 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |