Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,362 | 0 | C | T | 46.3% | 37.4 / 29.9 | 41 | intergenic (+52/‑51) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base C (5/17); new base T (19/0); total (24/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.15e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NZ_CP009273/2380276‑2380451 | 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:1078423/90‑1 (MQ=255) ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGa < 1:2025463/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) 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) gCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAg < 2:1856610/90‑1 (MQ=255) ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt < 2:1403993/90‑1 (MQ=255) tttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTg > 2:1201057/1‑90 (MQ=255) attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt < 1:2122431/90‑1 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 1:1724479/1‑90 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 1:16668/1‑90 (MQ=255) tCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGa > 2:957545/1‑90 (MQ=255) cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAATGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa < 2:1030446/90‑1 (MQ=255) cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa < 2:206141/90‑1 (MQ=255) cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa < 2:1192552/90‑1 (MQ=255) cccTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa < 2:740220/90‑1 (MQ=255) | ATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAA > NZ_CP009273/2380276‑2380451 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |