philippdre / omniCLIP

omniCLIP is a CLIP-Seq peak caller
GNU General Public License v3.0
15 stars 9 forks source link

One issue that require your attention. #29

Open Ruitulyu opened 1 year ago

Ruitulyu commented 1 year ago

Please use STAR to map your par-clip data, bowtie2 generated bam files will report the MD/NM tag error. nohup STAR --runThreadN 10 --genomeDir hg19_STAR/ --readFilesCommand zcat --readFilesIn PAR-CLIP.rep1.R1_val_1.fq.gz PAR-CLIP.rep1.R2_val_2.fq.gz --outFileNamePrefix PAR-CLIP.rep1 --outSAMtype BAM SortedByCoordinate--outSAMattributes MD NM &