| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 832,042 | 0 | C | G | 67.9% | ‑2.8 / 24.7 | 28 | A79G (GCG→GGG) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
| Reads supporting (aligned to +/- strand): ref base C (1/8); new base G (19/0); total (20/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.90e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-05 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831960‑832131 | aTGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGcggcgggc > 2:330462/1‑88 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGggggggggggggag > 1:1486360/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGgggggggggggcgc > 1:306683/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGgggcggcgtc > 2:1019963/1‑86 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGtgat > 2:124497/1‑87 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:1073168/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:1057463/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:15054/1‑90 (MQ=255) ccAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAGTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTc < 1:190871/90‑1 (MQ=255) aTGTACGCTCGTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTC‑GG‑TGGCg > 1:1013513/1‑90 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcggc < 1:1505523/73‑1 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGggcggc > 2:1505523/1‑73 (MQ=255) aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCC‑CA‑Tgggg > 1:1075597/1‑88 (MQ=255) cTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTCGGA‑TGG‑GATCgcgc > 1:90544/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTC‑CA‑TGGGGGTCGCGCGTTTTGTCAg > 1:645068/1‑90 (MQ=255) cAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCC‑CA‑AGGCGATCGCCCATTTTGTGAg > 2:687815/1‑90 (MQ=255) ggggCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCCCCCCC‑GA‑TGGCGATCGCGCATTTTGTCAggg > 1:121147/1‑89 (MQ=255) cACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTcgcg < 2:1372841/90‑1 (MQ=255) aCCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCCC‑GG‑TGGCGATCGGGCATTTTGTCAGGGCGCgg > 1:1237594/1‑90 (MQ=255) ccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGGACCCTC‑CA‑TGGGGGTCGCGCCTTTTGTCAGGGCGGGGc > 2:423145/1‑90 (MQ=255) aCCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACCCCC‑GA‑TGGGGATCGCGCCTTTTGTCAGGGCGCGGCCCATGAAGCGa > 1:1102814/1‑90 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:1326158/90‑1 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:124497/90‑1 (MQ=255) tcatGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGGCACGCTC‑GA‑TGGCGGTCGCCCATTTTGTCAGGTCGCGGCACATGAAGCGGTggg > 1:1559029/1‑89 (MQ=255) catGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCACGCTC‑GA‑TGGCGGTCGCGCATTTTGGCAGGTGGCGGGACAAGAAGCGAAGgcg > 2:287205/1‑90 (MQ=255) atGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGGGATGgggc > 1:934651/1‑90 (MQ=255) tGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCt < 1:1019963/90‑1 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCCCGCCC‑CAGTGGGGA‑CGCGCCCTTTGTGCGGGCGCGGCACAAGAAACGATGGCGCCGGGGGTTggg > 2:524666/1‑88 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACGCCC‑GA‑TGGCGATCCCGCCTTTTGGGAGGGGGCGGGGCCTGAAACGATGGGGGTGGGGGTTggg > 2:976087/1‑88 (MQ=255) gTTAAAGAGGCGGGGGGGGCGGGCACCCCC‑CC‑TGGGGGTCGCGCCTTTTGTCAGGGCGCGGCCCACGAAGGGATGGGGCGGGGGGGTgag > 2:564060/1‑90 (MQ=255) gTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTgag < 2:1013513/90‑1 (MQ=255) ttAAAGAGGCGGGGGGGGCGGGCACGCTC‑GA‑TGGCGATCGCGCATTTTGTCAGGGCGCGGCACATGAAGCGATGGGGCTGGGGGTTggga > 1:396849/1‑90 (MQ=255) ggCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCa < 1:1057463/90‑1 (MQ=255) cGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATc < 1:1317484/90‑1 (MQ=255) cggcggTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACg < 1:423145/90‑1 (MQ=255) | ATGATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831960‑832131 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |