featureCounts -p -T 20 -a ERVL.saf -o ERVL.counts xxx.bam xxx.bam ....
emmm....
有点问题:featureCounts 说丫的不识别,就用bedtools试下吧!
ref:http://www.bio-info-trainee.com/745.html
构建bed:
bedtools multicov -bams T6-DMSO-1_L3_329329.R.uniq.bam T6-DMSO-2_L3_330330.R.uniq.bam T6--1_L3_327327.R.uniq.bam T6-4OHT-2_L3_328328.R.uniq.bam -bed L1.bed > L1.counts
网友评论