| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 832,048 | 0 | T | G | 61.8% | 22.1 / 38.0 | 35 | V81G (GTC→GGC) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
| Reads supporting (aligned to +/- strand): ref base T (1/12); new base G (21/0); total (23/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.37e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831962‑832131 | gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:885282/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 1:20664/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:1122714/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:254895/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:56515/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGCCAc > 1:510403/1‑90 (MQ=255)gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCCGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:878153/1‑90 (MQ=255) tCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGc < 1:151320/90‑1 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCGCTCGCTGGCg > 2:298001/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGGTgggg > 2:707119/1‑88 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCGCGCCCCTTGGGGATGgcgc > 1:199349/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGGGGgcgcgc > 2:602295/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGCGGGGCACCCCCCCTGGCGATCgcgc > 1:447375/1‑90 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCGCCCGGAGGCGGTCgcgccat > 2:670362/1‑87 (MQ=255) ggAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACGCCCGATGGGGATCGCGCAtttt > 1:4645/1‑89 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCCCCCCATGGCGAGCGCGCCTTTTGGCAg > 1:1107868/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGGTGGGGGGCGCGCATTTGGTccg > 2:979686/1‑88 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGGGGTCGCGCATTTTGTCAg > 2:1122267/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGGGGTCGCGCATTTTGTCAg > 2:529884/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGGATTTGGGCAg > 2:384676/1‑90 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGc < 2:783165/90‑1 (MQ=255) cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACAt < 2:4645/90‑1 (MQ=255) aTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGaa < 2:634007/90‑1 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCGATGGGGGGCGCCCATTTGGGCGGGTCGGGGCCCATTAAgag > 2:710045/1‑88 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:629447/90‑1 (MQ=255) catcatGCCAAAACTGTTAAAGAGGCGGGGGCGGCGGGCACCCTCGATGGCGATCGCGCCTTTGGTCAGGTCGCGGCACCTGAAGCGGTg > 2:438612/1‑90 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt < 1:687282/90‑1 (MQ=255) gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg < 1:1122267/90‑1 (MQ=255) gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg < 1:979686/90‑1 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATTGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGGGTgag > 2:801722/1‑90 (MQ=255) agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc < 2:1035821/90‑1 (MQ=255) gggggCGGCGGGCACGCCCGATGGCGATCGCGCATTTTGGCGGGGCGCGGGCCATGAAGCGGTGGGGGTGGGGGGTGAGGAAAGGCATCa > 1:1100458/1‑90 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc < 1:801722/90‑1 (MQ=255) ggcggcGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCa < 1:424014/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 2:218589/90‑1 (MQ=255) | GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > 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: ‑ |