Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,276,090 | 0 | T | G | 53.8% | 26.1 / 30.6 | 26 | G252G (GGT→GGG) | agaS | AgaS family sugar isomerase |
Reads supporting (aligned to +/- strand): ref base T (2/10); new base G (14/0); total (16/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.24e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGGAAACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGT > NZ_CP009273/3276003‑3276170 | gTGGAAACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGaa > 1:119036/1‑90 (MQ=255) cGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑ggcg < 1:316995/90‑1 (MQ=255) ggATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑ggcgg < 1:14374/90‑1 (MQ=255) ggATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑ggcgg < 1:4315/90‑1 (MQ=255) ggATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑ggcgg < 1:409378/90‑1 (MQ=255) ggATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑ggggg > 2:400130/1‑90 (MQ=255) ggATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑ggcgg > 1:137899/1‑90 (MQ=255) aTCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGCGGGCTTTTAtt > 1:418348/1‑89 (MQ=255) tctcGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAAATGGGGGGCCtttttta > 2:231465/1‑87 (MQ=255) cGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGCGGCCTTTTATGATtc > 1:167120/1‑90 (MQ=255) ggTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGGGGCCTTTTTTGATTCTCCCAcc > 1:217333/1‑90 (MQ=255) ggCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGCGGCCTTTTTTGATTCTCCCACCCGa > 2:154025/1‑90 (MQ=255) gcgcAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAAAC‑GGCGGCCTTTTTTGATTCTCCTACCGCGTTTCGTCATg > 1:155792/1‑90 (MQ=255) gcgcAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGGGGCCTTTTTTTATTCTCCCACCGGGTTCCGTCCTg > 2:395869/1‑90 (MQ=255) gcgcAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCCTCATg > 1:396142/1‑90 (MQ=255) gcgcAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAAAT‑GGCGGCCTTTTTTTATTTTCCCACCCGGTTCCGTCATg > 2:434655/1‑90 (MQ=255) caCGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCaa < 1:63356/90‑1 (MQ=255) aCGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCaaa < 2:361623/90‑1 (MQ=255) aGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGc < 2:397681/90‑1 (MQ=255) cGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGCGGCCTTTTTTGATTCTCCCACCGGGTTCCGGCATGGACCAAAAACGCTgg > 2:344476/1‑90 (MQ=255) cGGCGCTGAAAGTGCTGGAACTGACGGCGGGGAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCCTCATGGACCAAAATCGCTgg > 2:469474/1‑90 (MQ=255) gAAAGTGCTGGAACTGACGGCGGGTAAAAT‑GGCGGCCTTTTTTGATTTTCCACCCCGGTTCCCGCCTGGGCCCAAAACGCGGGGCGGTGa > 2:254330/1‑90 (MQ=255) gTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGaa < 1:86484/90‑1 (MQ=255) gTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGaa < 1:278402/90‑1 (MQ=255) aCTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCtggt < 1:63946/90‑1 (MQ=255) gACGGCGGGGGAACT‑TGCGGCCTTTTTTGGTTCTCCCACCGGGTTCCGGCTTGGAACAAAAACCCCGGGCGATGACGAAACGCtgggggg > 2:404714/1‑89 (MQ=255) | GTGGAAACGGATCGTTTATCTCGGTAGCGGTGGCTTACAGGGCGCAGCACGCGAGTCGGCGCTGAAAGTGCTGGAACTGACGGCGGGTAAACT‑GGCGGCCTTTTATGATTCTCCAACCGGATTCCGTCATGGACCAAAATCGCTGGTCGATGACGAAACGCTGGTGGT > NZ_CP009273/3276003‑3276170 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |