| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 832,048 | 0 | T | G | 63.2% | 28.2 / 41.9 | 38 | V81G (GTC→GGC) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
| Reads supporting (aligned to +/- strand): ref base T (4/10); new base G (24/0); total (28/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.12e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831961‑832131 | tGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCa > 2:1403880/1‑90 (MQ=255)tGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCa < 2:783195/90‑1 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAAGCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:1244931/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:755160/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:546491/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 1:284604/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:197860/1‑90 (MQ=255) cAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCTAAAACCGTTAAAGAGGCGGGGGCGGGGGTCACGCCCg > 2:269784/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGATGGCg > 2:767583/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGATGGCg > 2:714205/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTccc > 1:1084163/1‑77 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCCATGGCg > 1:1059022/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCAc < 2:1084163/79‑1 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 1:767583/90‑1 (MQ=255) gCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCCATGGCGGGcgcg > 1:1101168/1‑90 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCGCGCCCGATGGCGGTCgcgc > 1:1136624/1‑90 (MQ=255) ccTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑t < 1:632775/90‑1 (MQ=255) cTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGGTCGCGCA‑tt > 2:632775/1‑90 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACAATCGGTGGGGAGCGCGCA‑ttttggt > 2:1263509/1‑86 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCA‑TTTTGTCAg > 2:30463/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGGTCGCGCA‑TTTTGTCAg > 1:160079/1‑90 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCCCGCCCGATGGGGATCGCGCC‑TTTTGTCAGGt > 1:426199/1‑90 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCa > 2:147256/1‑59 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCa < 1:147256/59‑1 (MQ=255) gCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTcgc < 1:302342/90‑1 (MQ=255) ttAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAg > 2:301899/1‑90 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCCATGGCGGTCGCGCC‑TTTTGTCAGGGCGCGGCACATGAAGCg > 2:1062284/1‑90 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGAt < 2:40292/90‑1 (MQ=255) tcatGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACCCCCGATGGCGGTCGCGCC‑TTTTGTCAGGGCGCGGGGCCTGAAGCGATggg > 2:1507082/1‑89 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCTCGATGGGGATCGCGCA‑TTTTGTCAGGTCGCGGGACATGAAGCGAGGGGGCt > 2:861920/1‑90 (MQ=255) ccAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTgg < 1:550517/90‑1 (MQ=255) aaaaCCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGGTCGCGCA‑TTTGGTCAGGGCGCGGCACATGAAGCGATGGGGCTgggg > 2:715983/1‑90 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGGCGCGGCACATGAAGCGATGGGGCTGGGGGTTgag > 2:182816/1‑90 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGGCACGCCCGATGGGGATCGCGCA‑TTTTGTCAGGGCGGGGCACCAGAAACGAGGGGGCGGGGGGGTgag > 1:704208/1‑90 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 1:1403880/90‑1 (MQ=255) gggggCGGCGGGCACGCTCGGTGGCGATCGCGCATTTTGGTTAGG‑CGCGGCACACGAAGCGGTGGGGCTGGGGGGTGGGAAAACGCGTCa > 1:794112/1‑90 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:861920/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCACTGGGGATTGAGAAAGCGCATCAGCACg > 1:391274/1‑90 (MQ=255) | TGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831961‑832131 |
| 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: ‑ |