Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,991,529 | 0 | A | T | 52.6% | 12.7 / 25.5 | 19 | T63S (ACC→TCC) | uvrD | DNA helicase II |
Reads supporting (aligned to +/- strand): ref base A (0/9); new base T (10/0); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGC > NZ_CP009273/3991440‑3991606 | gTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTa < 2:612377/90‑1 (MQ=255) ggTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGACGTTTTcc > 1:117233/1‑90 (MQ=255) gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGCGGTTTTCAACaaaac > 1:951154/1‑88 (MQ=255) gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGAGGGTTTTCAACaaaaa > 1:461070/1‑88 (MQ=255) gcgcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGACGGTTTTCAACaaaaa > 1:96261/1‑88 (MQ=255) gcgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGCGGTTTCCaaaaaaaagg > 2:230886/1‑86 (MQ=255) tGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAAcagcgggggggg > 1:192690/1‑87 (MQ=255) aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc < 2:699530/90‑1 (MQ=255) aTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc < 2:297253/90‑1 (MQ=255) aTCGCCTGGTTGCTGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGGGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCAt < 2:679042/90‑1 (MQ=255) gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAAGTTTTCCAACAAAACGGCGGCGGGGGGGCGTCATCGt > 1:299799/1‑90 (MQ=255) ggTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTTCCAACAAAACGGGGGGGGGGGGGGGGCCTCGTATCg > 1:423589/1‑90 (MQ=255) gatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCa < 2:210374/90‑1 (MQ=255) tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAAc < 2:252445/90‑1 (MQ=255) aaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc < 1:175020/90‑1 (MQ=255) aaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCAc < 2:459040/90‑1 (MQ=255) tCGCCATACTCGATTATGGCGGTGACGTTTTCCAACAAAACGGCGGCGGGGGTGGGGCATCGTATCGGGGGACAGATGGGGGCGAGcccg > 1:943300/1‑88 (MQ=255) ccATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGc < 2:555894/90‑1 (MQ=255) ttATGGCGGTGACGTTTTCCAACAAAACGGGGGGGGGGGTGGGGCCTCGTATCGGGGGACCGATGGGGACGAGCCCGGGCGGGGTGTggg > 1:888293/1‑90 (MQ=255) gCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGc < 2:173101/90‑1 (MQ=255) | GTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGC > NZ_CP009273/3991440‑3991606 |
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: ‑ |