Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,331 | 0 | T | G | 56.2% | 10.2 / 15.0 | 16 | intergenic (+21/‑82) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base T (2/5); new base G (9/0); total (11/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.81e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTA‑TTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGG > NZ_CP009273/2380251‑2380417 | tATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCGGaaaaa > 1:272948/1‑90 (MQ=255) tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGTCACaaaaaaa > 2:130949/1‑88 (MQ=255) tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGGCCGaaaaaaa > 2:265754/1‑88 (MQ=255) tCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCCGaaaaaaa > 2:481594/1‑88 (MQ=255) tGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCAGAAAAAACGACATGCACATTTTT‑TTTTTCCTc > 2:369395/1‑90 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTA‑TTTTCCCTCATc < 1:240122/90‑1 (MQ=255) ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCTGAAAAAAAGACATGCACATTTTA‑TTTTTCCTCctctt > 2:19402/1‑89 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCAGAAAAAAAGACAAGCACATTTTT‑TTTTTCCTCTTCTGtt > 2:451248/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTATCTGGAAAAAG‑GAAAAACACATTTTT‑TTTTTCCTCTCATGttttgtt > 1:535217/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GGCCAGAAAAAACGACATGCACATTTTT‑TTTTTCCTCCTCtttttttaa > 2:305690/1‑87 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTA‑TTTTCCCTCATCTGTTTTAAGATGCGa < 1:484276/90‑1 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGG‑GTGCAGAAAAAAAGACAAGCACATTTTT‑TTTTTTCTCATCTGTTTTTAGATGGGAAACAAGAAAt > 1:363257/1‑90 (MQ=255) taataaAAATTTCTCGAGGCGG‑GGCCAGAAAAAACGACATGCACAATTTT‑TTTTTCCTCCTCTTTTTTTAAAAGCGAAAAAAGaaaatta > 2:166342/1‑86 (MQ=255) aaaaTTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTATTTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg < 1:90517/90‑1 (MQ=255) aTTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTA‑TTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga < 2:360570/90‑1 (MQ=255) gg‑gTCCAGAAAAAGCGACATGCACATTTTA‑TTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATggg < 2:136296/90‑1 (MQ=255) | TATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGG‑GTCCAGAAAAAGCGACATGCACATTTTA‑TTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGG > NZ_CP009273/2380251‑2380417 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |