Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,380,341 | 0 | G | A | 54.5% | 3.5 / 29.2 | 22 | intergenic (+31/‑72) | yfbL/yfbM | M28 family metallopeptidase/YfbM family protein |
Reads supporting (aligned to +/- strand): ref base G (1/9); new base A (12/0); total (13/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCT > NZ_CP009273/2380264‑2380428 | gacAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑at < 1:71098/88‑1 (MQ=255) tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑att < 2:737127/90‑1 (MQ=255) aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑atttt < 2:418623/90‑1 (MQ=255) gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTccc < 1:577389/90‑1 (MQ=255) gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGCAC‑ATTTTTTTTTTCCtttc > 2:71098/1‑87 (MQ=255) tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATc < 1:950891/90‑1 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCCTCtttt > 2:548255/1‑90 (MQ=255) ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAAAGAAAAGAAA‑AATTTTTTTTTTCTCATCTGtt > 1:687241/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAACG‑CAAGAAC‑CCTTTTTTTTTCCCCATCtttttttaa > 2:857756/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGT‑CCCGAAAAAACGACAAGCAC‑ATTTTTTTTTTCCTCATCCGTTTTTAg > 2:104557/1‑90 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGAAAAGCAC‑ATTTTTTTTTTTCTCCTCTGtttttat > 1:283523/1‑87 (MQ=255) tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCAGAAAAAACGACAAGGAC‑ATTTTTTTTTTCCTCCTCttttttttg > 2:890676/1‑87 (MQ=255) caacaGTAAATAATAAAAATTTCTCGAGGCGGGG‑GCCGAAAAAAAAACAAGCAC‑AATTTTTTTTTCCTCCTCTGTTTTTAAATGGGaaaa > 1:725209/1‑89 (MQ=255) acaGTAAATAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAAAGACAAGCACAATTTTTTTTTCCCCCTTTTTTTTTAA‑ATGGGaaaaaa > 2:649050/1‑87 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGAAAt > 2:1235508/1‑90 (MQ=255) aaaTAATAAAAATTTCTCGAGGCGGGG‑CCCGAAAAAACGACAAGCAC‑ATTTTTTTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGaaaa > 2:138038/1‑89 (MQ=255) ataaAAATTTCTCGAGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc < 1:28095/90‑1 (MQ=255) ttCTCGAGGCGGGG‑CCAGAAAAAACGACATGCAC‑ATTTTTTTTTTCCTCCTCTTTTTTTAAATGGGAAAAAAAAAAATCAAAGGGGGgaa > 2:1302250/1‑90 (MQ=255) aGGCGGGT‑CCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataaa < 2:1240825/90‑1 (MQ=255) aaaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 1:104557/90‑1 (MQ=255) aaaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc < 2:1100269/90‑1 (MQ=255) aaGCGACATGCAC‑ATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt > 1:637002/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 < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |