Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 832,048 | 0 | T | G | 64.3% | 10.6 / 20.6 | 28 | V81G (GTC→GGC) | hcxB | hydroxycarboxylate dehydrogenase HcXB |
Reads supporting (aligned to +/- strand): ref base T (1/9); new base G (18/0); total (19/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831962‑832131 | gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGtgat > 2:124497/1‑87 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:15054/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 2:1057463/1‑90 (MQ=255) gATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGCGGGCAc > 1:1073168/1‑90 (MQ=255) ccAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAGTTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTc < 1:190871/90‑1 (MQ=255) aTGTACGCTCGTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGGGGGGGGCACGCTC‑GG‑TGGCg > 1:1013513/1‑90 (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:124497/90‑1 (MQ=255) ccATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAt < 1:1326158/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) | GATCCCAAGCTATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTC‑GA‑TGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACG > NZ_CP009273/831962‑832131 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |