| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_017635 | 3,111,916 | 0 | C | G | 30.0% | 50.3 / 26.0 | 30 | A449A (GCC→GCG) | vgrG | type VI secretion system tip protein VgrG |
| Reads supporting (aligned to +/- strand): ref base C (13/8); new base G (6/3); total (19/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTACGGCCTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCGTAAAACACG > NC_017635/3111868‑3111965 | cTACGGCCTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑Ct > 1:81489/1‑50 (MQ=32) tACGGCCTGCATCTGCCGCTGCTGGCGGGAACGGGAAACAGCAACGc‑cct > 2:1844136/1‑48 (MQ=11) cGGCCTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑Cttt > 1:1444977/1‑49 (MQ=32) ggCCTGCATCTGCCGCTGCTGGCGGGGACGGGAGGAAGCATCGC‑Cttttt > 2:2065376/1‑48 (MQ=11) ggCCTGCATCTGCCGCTGCTGGCGGGGAAGGAAGGCAACAACGC‑CCTTgg > 2:1110660/1‑49 (MQ=2) ggCCTGCATCTGCCGCTGCTGGCGGGAAGAGAAGGCAGCACCGCGCTTTg > 2:375032/1‑50 (MQ=2) gCCTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAg > 1:1920446/1‑50 (MQ=32) ccTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑GTTTGAg > 1:352772/1‑49 (MQ=12) ccTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑GTTTGAg > 1:467442/1‑49 (MQ=12) cTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAgg > 1:5693/1‑49 (MQ=32) cTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGa < 2:698818/50‑1 (MQ=32) tGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑GTTTGAgg > 1:1509532/1‑48 (MQ=12) cATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑GTTTGAGGaa < 2:427085/48‑1 (MQ=12) cATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTTAGGAAgg > 1:33515/1‑50 (MQ=12) aTCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGc < 2:83964/50‑1 (MQ=32) gCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAAcc > 1:1279480/1‑50 (MQ=32) ccGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑GTTTGAGGAAGGCAAc < 1:94488/48‑1 (MQ=12) ccGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAAcc < 2:690130/49‑1 (MQ=32) gctgctGGCGGGAACGGAAGTCAGCATCGC‑GTTTGAGGAAGGCAACCCg > 1:193130/1‑49 (MQ=12) gCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAACCCGgaccgacc > 1:1885551/1‑50 (MQ=32) cGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCt < 2:48121/50‑1 (MQ=32) gAACGGAAGTCAGCATCGC‑GTTTGAGGAAGGCAACCCGGACCGACCTta > 1:680932/1‑49 (MQ=12) aaCGGAAGTCAGCATCGC‑GTTTGAGGAAGGCAACCCGGACCGACCTtata > 1:159771/1‑50 (MQ=12) aaCGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTtata > 2:1436752/1‑50 (MQ=32) aGTCAGCATCGC‑GTTTGAGGAAGGCAACCCGGACCGACCTTATATAGCAg < 2:1589720/50‑1 (MQ=11) aGTCAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACATTATATTGCAg < 1:1783487/50‑1 (MQ=255) gTCAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAgg > 2:1984393/1‑50 (MQ=40) cAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCg > 2:1569841/1‑50 (MQ=40) aGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCGt < 1:89327/50‑1 (MQ=40) aGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCGt < 2:1026088/50‑1 (MQ=40) gC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCGTAAAac < 2:1328969/49‑1 (MQ=40) cTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCGTAAAACACg < 2:168203/50‑1 (MQ=38) | CTACGGCCTGCATCTGCCGCTGCTGGCGGGAACGGAAGTCAGCATCGC‑CTTTGAGGAAGGCAACCCGGACCGACCTTATATTGCAGGCGTAAAACACG > NC_017635/3111868‑3111965 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |