一开始感觉要直接做三组平均，写了个peakaverage.sh 暂时没用上 #!/bin/bash# Define input and outp...

环境还是m6a的conda activate m6acd /home/data/t210424/ac4c call了peak结果不是很理想暂时没...

环境还是m6a的 conda activate m6acd /home/data/t210424/ac4c rRNA测了比例竟然是0%所以没去 ...

环境还是m6a的 conda activate m6acd /home/data/t210424/test QC mkdir ./QCclean...

直接改了/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...