piRNA reads数分析
2019-01-02 本文已影响0人
njmujjc
##### -a 输出所有multipy alignment ####
bowtie2 -p 18 -a -x /media/pc/disk2/jjc/piPipes-1.5.0/common/mm9/mm9.piRNAcluster -U pac-zcchc8-control-2-1_trim.fq-S control2-1.sam
samtools sort -@ 18 -o control2-1.bam control2-1.sam; samtools index -@ 18 control2-1.bam
##### 统计reads数目 ####
samtools idxstats control2-1.bam> control2-1.txt
##### normalize 每挖掘百万个miRNA时piRNA的reads数目####
piRNA reads数目 * 1000000 / 总miRNA数目
