| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 832,042 | 0 | C | G | 48.5% | 26.9 / 11.3 | 33 | A79G (GCG→GGG) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
| Reads supporting (aligned to +/- strand): ref base C (2/15); new base G (16/0); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.47e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831956‑832131 | tGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcggc < 2:430341/90‑1 (MQ=255) gCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcgg > 2:626095/1‑90 (MQ=255) gATCCCAAGCTCTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCAc < 2:23725/90‑1 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGgggggggggggccc > 2:1165913/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggc < 1:1374955/81‑1 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:866490/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGGGGGCAc > 2:950149/1‑90 (MQ=255) cccAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCt < 2:390344/90‑1 (MQ=255) cAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCg > 2:26743/1‑90 (MQ=255) tGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACCCCCGTTGGCgg > 2:403888/1‑89 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 2:1401942/90‑1 (MQ=255) cGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATcgc < 2:1398456/90‑1 (MQ=255) ccTGGAGTCAGGGGCACATGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑t < 2:310927/90‑1 (MQ=255) ggAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGGACCCCCGATGGGGGTCGCGCA‑TTTg > 2:360684/1‑90 (MQ=255) aGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGt < 2:1392879/90‑1 (MQ=255) gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGCGATCGCGCAGTTTggg > 1:112602/1‑89 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCCCATGGCGATCGCGCC‑TTTTGTCAg > 2:1364886/1‑90 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCCCCCCCACTGGGGATCGCGCA‑TTTTGTAAggg > 1:318128/1‑89 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGc < 1:360684/90‑1 (MQ=255) tGCAAATTAACCATCATGACAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGGGGTCGCGCA‑TTTTGTCAGGGCGCGGcac > 2:1126279/1‑90 (MQ=255) cAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGGGGTCGCGCC‑TTTTGTCAGGGCGCGGCGCAt > 1:642913/1‑90 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCt > 2:1116086/1‑47 (MQ=255) aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCt < 1:1116086/47‑1 (MQ=255) catcatGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCCCATGGCGAACGCGCC‑TTTTGTCAGGGGGCGGGACATGAAGCGATg > 1:58849/1‑90 (MQ=255) tcatGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTGTCCCCCCATTGGGAACGCGCA‑TTTTTTGAGGGGGCGGCACATGAAGCGGTggg > 1:562755/1‑89 (MQ=255) atGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTCGATGGCGGTCGCGCC‑TTTTGTCAGGGCGCGGGCCATGAAACGAgggggg > 2:1154030/1‑87 (MQ=255) gCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTg < 2:642913/90‑1 (MQ=255) aaCCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGAt < 1:677662/90‑1 (MQ=255) gTTAAAGAGGCGGGGGGGGGGGGCACGCTCGATGGGGATCGCGCA‑TTTGGGCAGGGCGCGGGACATGAAGCGATGGGGCTGGGGGGTgag > 2:434492/1‑90 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGGTCGCGCA‑TTTTGGCAGGGGGCGGGCCATGAAGCGGTGGGGCCGGGGGGTgag > 1:1199152/1‑90 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCCCGATGGCGATCGCGCA‑TTTTGTCAGGGCGCGGGACATGAAGCGGTGGCGCTGGGGGGTgag > 2:54682/1‑90 (MQ=255) tAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGaa < 1:54682/90‑1 (MQ=255) ggCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCa < 1:434492/90‑1 (MQ=255) cGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 2:1209459/90‑1 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc < 1:1154030/90‑1 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc < 2:215060/90‑1 (MQ=255) gggCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGc > 2:676071/1‑90 (MQ=255) ggcggcGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCa < 1:614007/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGc < 2:467064/46‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGc > 1:467064/1‑46 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:676071/90‑1 (MQ=255) cggcggTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:1126279/90‑1 (MQ=255) | TGGCATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCA‑TTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831956‑832131 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |