直接改了/home/data/t210424/m/pcnolambda.sh的output的文件夹 peakcall0423nolambda和p...

peakcallq001 peakcall0422p005 peakcall0422q005 peakcall0423nolambda 保留no...

callpeak for i in {1..9}; do macs2 callpeak \ -t "${align_dir}/RIP...

文件夹/home/data/t210424/m 1.QC是传统QC 2.trim /home/data/t210424/m/trim.sh # ...

不过滤，用测序公司直接给的BAM文件callpeak #MACS2 软件 环境还是 conda activate m6a ~/m6a ./sin...

整个图> 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...