整个图> setwd("F:/XZHMU/huangyue/metagene plot")> bed_file1 <- "Sham_3_summ...

所有bug暂且不表，只记录有用的 环境&包，因为是新电脑，全是新装的if (!require("BiocManager", quietly = ...

library(ggseqlogo) SNI截motif3的4-8位E-value:7.8e-016 prob_matrix <- matrix...

发现当初学的DREME被淘汰了 试了用streme写了一个，先把所有sequences.fasta移动到/home/data/t210424/m...

顺利运行的peakcall3.sh #!/bin/bashinput_dir="/home/data/t210424/m6a/filtered"...

先对单个文件做了peak call， /home/data/t210424/m6a/filtered/singlepeakcall.sh #!/...

/home/data/t210424/m6a/qs.sh做个bam文件过滤 #!/bin/bashinput_dir="/home/data/t...