



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
