Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 832,048 | 0 | T | G | 60.0% | 41.8 / 61.4 | 50 | V81G (GTC→GGC) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
Reads supporting (aligned to +/- strand): ref base T (3/17); new base G (30/0); total (33/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC > NZ_CP009273/831962‑832130 | gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:991008/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:1354316/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:557538/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCAc > 2:1612956/1‑90 (MQ=255) gtGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCg > 2:693419/2‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGCg > 1:1542621/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCCTGGCg > 1:1870049/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACACGCGATGGCGGTCgc‑gc > 1:1880066/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCCTGGCGATCgc‑gc > 1:1847223/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCgc‑gc > 2:1278055/1‑90 (MQ=255) tCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGGCACGCTCGATGGCGATCgc‑gcc > 2:701830/1‑89 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACGCCCCATGGCGATCGC‑GCCtt > 2:1377636/1‑90 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGC‑GCAtt > 2:728637/1‑90 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCAATCGC‑GCATTTGGTc < 1:1446016/90‑1 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCCCGCGATAGGGGTCGC‑GCCTTTGGttgg > 1:1155685/1‑87 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACGCCCCCTTGCGGGCGC‑GCACTTGGTGAg > 1:1796334/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGC‑GCATTTGGTCAg > 2:763465/1‑90 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGGGGTCGC‑GCATTTTGTCAggg > 2:998183/1‑89 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGGTCGC‑GCATTTGGTCAggg > 1:628333/1‑89 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCGGTGGGGAGCGC‑GCATTTGGTTAggg > 1:951676/1‑89 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCGTGGGGAGCGC‑GCCATTGGTCAGGGGGCggg > 1:1112796/1‑89 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGATCGC‑GCATTTGGGCAGGGCGCGGc > 1:495965/1‑90 (MQ=255) cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACAt < 1:314134/90‑1 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGATGGCGGTCCC‑GCCTTTTGTCAGGGCGCGGGACAAGAAGCg > 1:339350/1‑90 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGGCCGCTCGGAGGTGGTCGCGGGATTTGG‑CAGGGCGCGGCACATGAAACg > 2:1673134/1‑90 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:701830/90‑1 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:236598/90‑1 (MQ=255) atGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCCATGGCGATCCC‑GCCTTTTGGTAGGGCGCGGGCCATgaagggaaggggg > 1:912862/1‑87 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt < 1:65765/90‑1 (MQ=255) gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg < 1:991008/90‑1 (MQ=255) gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg < 2:495965/90‑1 (MQ=255) gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg < 2:339350/90‑1 (MQ=255) aaaaCCGTTAAAGAGGCGGGGGCGGCGGGCACGCTCGATGGCGATCGC‑GCATTTTGTCAGGGCGCGGCACATGAAGCGATGGGGCTgggg > 2:184661/1‑90 (MQ=255) gTTAAAGAGGCGGGGGGGGGGGGCACGCTCGATGGCGATCGC‑GCATTTTGTCAGGGCGCGGCCCATGAAGCGATGGGGCTGGGGGTTggg > 2:2079800/1‑88 (MQ=255) gTTAAAGAGGCGGGGGGGGGGGGCACCCCCGATGGCGGTCGC‑GCATTTTGGCAGGGGGCGGCACATGAAACGAGGGGGGTGGGGGtttgg > 2:1564149/1‑87 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCTCCATGGGGGTCGC‑GCATTTTGTCAGGGCGGGGCACATGAAGCGAGGGGGCTGGGGGTTggg > 2:1594731/1‑88 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTTGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGGTTgag > 2:442037/1‑90 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACAt < 1:9752/66‑1 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACAt > 2:9752/1‑66 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:632521/74‑1 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGAt > 2:632521/1‑74 (MQ=255) agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAgc < 1:184661/90‑1 (MQ=255) ggCGGGGGGGGCGGGCACCCCCGATGGCGATCGC‑GCGTTTTGTCAGGGGGCGGCACAAGAAGCGATGGGGCCGGGGGGTTAGAAAACGCa > 1:582762/1‑90 (MQ=255) ggCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCa < 2:657452/90‑1 (MQ=255) gCGGGGGGGGCGGGCACGCTCGATGGCGATCGC‑GCATTTTGTCAGGGCGCGGCACATGAAGCGGTGGGGGTGGGGGTTGAGAAAGCGCAt > 1:1050743/1‑90 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 2:428210/90‑1 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 1:1354316/90‑1 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc < 2:1112796/90‑1 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc < 1:327764/90‑1 (MQ=255) gcggcgGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAc < 1:998183/90‑1 (MQ=255) | GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGC‑GCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCAC > 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 < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |