整個(gè)圖> setwd("F:/XZHMU/huangyue/metagene plot")> bed_file1 <- "Sham_3_summ...

所有bug暫且不表，只記錄有用的 環(huán)境&包，因?yàn)槭切码娔X冕碟，全是新裝的if (!require("BiocManager", quietly = ...

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

發(fā)現(xiàn)當(dāng)初學(xué)的DREME被淘汰了 試了用streme寫了一個(gè)骑晶，先把所有sequences.fasta移動(dòng)到/home/data/t210424/m...

順利運(yùn)行的peakcall3.sh #!/bin/bashinput_dir="/home/data/t210424/m6a/filtered"...

先對(duì)單個(gè)文件做了peak call气破， /home/data/t210424/m6a/filtered/singlepeakcall.sh #!/...

/home/data/t210424/m6a/qs.sh做個(gè)bam文件過濾 #!/bin/bashinput_dir="/home/data/t...