1、軟件安裝
單獨(dú)為leaf cutter創(chuàng)建一個(gè)環(huán)境 R版本為3.4.1(筆者操作系統(tǒng)為linux red hat,嘗試了很多安裝版本和方法诞帐，唯一可行的是：)
conda create -n r3.4.1 -y -c conda-forge r-base=3.4.1
conda activate r3.4.1
conda install -y -c conda-forge -c bioconda -c davidaknowles r-leafcutter
conda install bioconda::regtools
2、Differential Splicing
Step 0. Alignment
STAR --runThreadN 20 --genomeDir /index/mm10_GRCm38.101_index \
--readFilesIn $trim/${name}_clean_1P.fq $trim/${name}_clean_2P.fq \
--outFileNamePrefix $align/${name} --twopassMode Basic --outSAMstrandField intronMotif\
--outSAMtype BAM SortedByCoordinate --outBAMsortingThreadN 20 --quantMode TranscriptomeSAM GeneCounts
Step 1. Converting bams to juncs
##官方code
for bamfile in `ls run/geuvadis/*chr1.bam`; do
echo Converting $bamfile to $bamfile.junc
samtools index $bamfile
regtools junctions extract -a 8 -m 50 -M 500000 $bamfile -o $bamfile.junc
echo $bamfile.junc >> test_juncfiles.txt
done
