Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,634,359 | 0 | A | G | 42.1% | 52.0 / 13.4 | 38 | G46G (GGA→GGG) | dtpB | dipeptide/tripeptide permease DtpB |
Reads supporting (aligned to +/- strand): ref base A (5/17); new base G (16/0); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGC > NZ_CP009273/3634223‑3634496 | tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTGGGtt > 2:342394/1‑139 (MQ=255) tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGAtt < 1:622/139‑1 (MQ=255) caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAaga < 1:419369/139‑1 (MQ=255) caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGGTGAAAAAGTTGGGTTTTCGCaa > 1:180474/1‑137 (MQ=255) cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCTTTTAACAGCTTGGATTTTCGCAAGGGCAgg > 2:411982/1‑139 (MQ=255) cccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAgg < 1:194451/139‑1 (MQ=255) aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTTAACAACTTGGGTTCTCGCACGAGGTGGC‑tc > 2:339460/1‑138 (MQ=255) aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑tt > 1:498342/1‑139 (MQ=255) ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTAAAAAGCTTGGGTTTTTGCCAGAACAGGC‑tttt > 2:500357/1‑139 (MQ=255) ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAACTTGGGTTTTCGCAAGAGCAGGC‑tttt > 2:325125/1‑139 (MQ=255) gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAACTTGGGTTTTCGCAAGAGCAGGC‑TTTTGTTACt > 2:467629/1‑139 (MQ=255) gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGGTTTTTGCAAGAGCACGG‑TTTTGTCACt > 1:308618/1‑139 (MQ=255) gctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAACAGCTTGGGTTTTCGCAAGAGCCGGG‑TTTTTTCACt > 2:131348/1‑139 (MQ=255) tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTAAAACGCCTCGGTTTTCGGAAGAGGAGGC‑TTTTTTTACttt > 1:337322/1‑139 (MQ=255) agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTgg < 2:520244/139‑1 (MQ=255) aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGc < 1:591917/139‑1 (MQ=255) cTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAACTTGGGTTTTTGCAAAAGCAGGCGTTTTTTAATTTTTggggttt > 2:591917/1‑134 (MQ=255) tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCCTTGAGTCTCGCAAGAAGCGGG‑TTTTTTTACTTTTGGTGGtttt > 1:374994/1‑139 (MQ=255) tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAACAGCTTGGATTTTCGGCAGAGCAGGG‑TTTTTTCACTTTTGGTGCtttt > 1:96412/1‑139 (MQ=255) gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGGTTTTCGCAAGAGCAGGG‑TTTTTTTACTTTTTGTGCTTTtgt > 2:408514/1‑138 (MQ=255) gCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTtg < 1:16119/138‑1 (MQ=255) aTTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTgcg < 2:498342/139‑1 (MQ=255) tttttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTgcgc < 2:64085/139‑1 (MQ=255) tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAGCTTTGGTTCTCCCAAAAGCAGGC‑TTTTGTCACTTTTTGTGGTTTTTCTGCGCTGGTCTATGGCCt > 1:201092/1‑139 (MQ=255) ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGc < 1:408514/136‑1 (MQ=255) gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt < 1:74985/139‑1 (MQ=255) gCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTAt < 1:325125/139‑1 (MQ=255) cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTAAACAGCTTGGGTTCTCGCCAGAGCAGGC‑TTTTTTCACTTTTGGTGCTTTTTCTGCGCTGGTCTATGGGCTCCTTTTCATTtgcgggcttg > 1:440007/1‑134 (MQ=255) cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTAAAAAGCCTGGGTTCTTGCAAGAGCGGGC‑TTTTTTTACTTTTTGTGCTTTTGCTGCGCTGGTCTATGGGCTCCTTTTCATTggcgggtctt > 1:145288/1‑134 (MQ=255) cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGGTTTTTTTTTGAACACCCTGGGATTTCCCCGGAGGATGC‑TTTTGTTTCTTTTTGTGATTTTTCGTCGCTGGTGTTCGGCCCCGTTTTCATTgggggggttt > 2:304538/1‑134 (MQ=255) cGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGCTTTCTTCGGTAAACCGCTTGGGTTCTCCCAAGAGAAGGG‑TTTTGACCCCTTTTGTGGGTTTGCTGCGCTGGGCTATGGCCTCATTTTCATTtgcgggcctt > 1:33982/1‑134 (MQ=255) ctacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGacc < 1:283861/137‑1 (MQ=255) ggCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACctgc < 1:487090/139‑1 (MQ=255) tACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACg < 1:88985/139‑1 (MQ=255) cTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCa < 1:131348/139‑1 (MQ=255) gTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAg < 1:467629/139‑1 (MQ=255) gCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtgtg < 1:258773/139‑1 (MQ=255) cTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGc > 2:97126/1‑139 (MQ=255) aTTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGc > 2:38394/1‑138 (MQ=255) | TGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGC > NZ_CP009273/3634223‑3634496 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |