Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,342 | 0 | C | A | 47.5% | 57.1 / 13.3 | 40 | intergenic (+32/‑71) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base C (4/17); new base A (19/0); total (23/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.94e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380264‑2380428 | cTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑at < 2:1561085/90‑1 (MQ=255) cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑attt < 2:1969377/90‑1 (MQ=255) gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGAAAAGAAC‑aatttttttt > 1:983251/1‑90 (MQ=255) aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCt < 2:526732/90‑1 (MQ=255) ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGAAAAGGACAATTTTTTTTTcccca > 2:2025463/1‑88 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑GCCGAAAAAAAGAAAAGCAC‑ATTTTTTTTTTTCTCTTCtttt > 1:2135755/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCCGAAAAAAAGACAAGCAC‑AATTTTTTTTTCCTCCTCtttttttaa > 2:2122431/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCCGAAAAAAAGACAAGCAC‑AATTTTTTTTTCCTCCTCTGTTTTTAg > 2:1361362/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCCAAAAAAAAGACACGCAC‑ATTTTTTTTTTCCCCCTCTTTTTTTAg > 2:1180694/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑TCAGAAAAAAAGAAAAACAC‑AATTTTTTTTTTCTCATCTGTTTTTAg > 1:9747/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAAAAAAGCAC‑AATTTTTTTTTCCTCCTCTTTTTTTAg > 1:1970915/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCCTCtttttttaa > 2:778416/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACAAGCCC‑ATTTTTTTTTTCCTCATCtttttttaa > 2:1786908/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGACATGCAC‑ATTTTTTTTTTCCTCCTCTGtttttaa > 2:1694879/1‑87 (MQ=255) tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCATAAAAA‑AGACAAAAAC‑AATTTATTTTTTCTCCTCTGTTTTTAGAt > 1:2008730/1‑90 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:2025463/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:1078423/90‑1 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAt < 1:1959713/90‑1 (MQ=255) aGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGGAAAAAAGAAAAACAC‑AATTTTTTTTTCCCCCTCTGTTTTTAGATGGGGAACAAGa > 1:946603/1‑90 (MQ=255) gTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAAAGAAAAGCAC‑AATTTTTTTTTCCTCCTCTGTTTTTAGATTGGAAAAATGaa > 2:688916/1‑90 (MQ=255) gTAAATAATAAAAATTTCTCGAGGAGGGT‑CCAGAAAAAGCGACATGCAC‑attttatt > 1:1855203/1‑57 (MQ=255) gTAAATAATAAAAATTTCTCGAGGAGGGT‑CCAGAAAAAGCGACATGCAC‑attttatt < 2:1855203/57‑1 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGG‑TCCCAAAAAAAGACAAGAAC‑ATATTATTTTTCCCCCTCTTTTTTTAAAAGGGAAAAAAGaaaa > 2:1225229/1‑89 (MQ=255) taataaAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGACAAGCAC‑ATTTTTTTTTTCCTCCTCTGTTTTTAAATGCGGAACAAGaaaatta > 2:1757241/1‑86 (MQ=255) taataaAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGCAC‑AATTTATTTTTTCTCCTCTGTTTTTAAATGCGAAAAAAGaaaatta > 1:506243/1‑86 (MQ=255) ataaAAATTTCTCGAGGCGGGT‑CCCGAAAAAAAGACATGCAC‑ATTTTATTTTTCCTCCTCTGTTTTTAGATGCGAAAAAAGAAAATCaaa > 2:983576/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc < 1:1425781/90‑1 (MQ=255) ataaAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGAAAAGCAC‑AATTTTTTTTTTCTCCTCTTTTTTTTAATGAGAAACATAAAAATCaaa > 1:1270005/1‑89 (MQ=255) taaAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg < 2:1286935/90‑1 (MQ=255) aaaTTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggag < 2:1825314/90‑1 (MQ=255) aTTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagga < 2:1970915/90‑1 (MQ=255) ttCTCGAGGCGGGT‑CCAGAAAAAAAGACATGCAC‑ATTTTTTTTTTCCTCATCTGTTTTTAGATGGGAAAAATGAAAATCAAAAGAGGaaa > 2:703032/1‑90 (MQ=255) ttCTCGAGGCGGGG‑CCAGAAAAAAAGACAAGAAC‑ATTTTTTTTTTTCTCCTTTGTTTTTTAATGGGAAAAAAAAAAATTAAAGGAGGaaa > 1:129182/1‑90 (MQ=255) ctcGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa < 2:586815/90‑1 (MQ=255) ctcGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa < 1:688916/90‑1 (MQ=255) ggCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaat < 1:703032/90‑1 (MQ=255) gggT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATggg < 2:744202/90‑1 (MQ=255) ggT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGa < 2:126217/90‑1 (MQ=255) t‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt < 2:1301099/90‑1 (MQ=255) aaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:1856610/1‑90 (MQ=255) | CTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380264‑2380428 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |