Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 832,042 | 0 | C | G | 61.1% | 7.3 / 34.7 | 36 | A79G (GCG→GGG) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
Reads supporting (aligned to +/- strand): ref base C (5/9); new base G (22/0); total (27/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.13e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831954‑832131 | aTTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcg > 1:783195/1‑90 (MQ=255) ttGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCgggggggg > 1:1598454/1‑90 (MQ=255) ttGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcgg < 2:692512/90‑1 (MQ=255) tGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCa > 2:1403880/1‑90 (MQ=255) tGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCa < 2:783195/90‑1 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAAGCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:1244931/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGcgggggcgggggggcc > 1:1519910/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGggggggggggggac > 2:1080772/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGgggggggggggccc > 2:302342/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcggggac > 2:894494/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcgggccc > 2:569443/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcggccga > 1:788092/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:546491/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 1:284604/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:197860/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCAc > 2:755160/1‑90 (MQ=255) cAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCTAAAACCGTTAAAGAGGCGGGGGCGGGGGTCACGCCCg > 2:269784/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGATGGCg > 2:714205/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCGATGGCg > 2:767583/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) | ATTGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > 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: ‑ |