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

所有bug暫且不表取具，只記錄有用的 環(huán)境&包，因為是新電腦腺毫，全是新裝的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寫了一個骂澄，先把所有sequences.fasta移動到/home/data/t210424/m...

順利運(yùn)行的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...