Difference between revisions of "HISAT2 pipeline"
From Crop Genomics Lab.
(Created page with "==mapping == /data1/program/hisat2/hisat2-2.2.1/hisat2 -q --phred33 -x /data1/ref/Gmax_275_v2.0.fa.hisat --summary-file SS2-2_1_hisat_summary -p 6 -1 {forward_input.fq.gz} -2...") |
(→sort/markdup) |
||
(3 intermediate revisions by one user not shown) | |||
Line 1: | Line 1: | ||
==mapping == | ==mapping == | ||
/data1/program/hisat2/hisat2-2.2.1/hisat2 -q --phred33 -x /data1/ref/Gmax_275_v2.0.fa.hisat --summary-file SS2-2_1_hisat_summary -p 6 -1 {forward_input.fq.gz} -2 {reverse_input.fq.gz} | samtools view -bS -q30 > {output}.bam | /data1/program/hisat2/hisat2-2.2.1/hisat2 -q --phred33 -x /data1/ref/Gmax_275_v2.0.fa.hisat --summary-file SS2-2_1_hisat_summary -p 6 -1 {forward_input.fq.gz} -2 {reverse_input.fq.gz} | samtools view -bS -q30 > {output}.bam | ||
+ | |||
+ | ==sort/markdup== | ||
+ | samtools sort {hisat.bam} -o {hisat.sort.bam} | ||
+ | java -jar /data1/program/picard-2.23.4/picard/build/libs/picard.jar MarkDuplicates I={hisat.sort.bam} O={hisat.sort.mkdup.bam} M={hisat.mkdup.metrics} |
Latest revision as of 07:27, 14 September 2020
mapping
/data1/program/hisat2/hisat2-2.2.1/hisat2 -q --phred33 -x /data1/ref/Gmax_275_v2.0.fa.hisat --summary-file SS2-2_1_hisat_summary -p 6 -1 {forward_input.fq.gz} -2 {reverse_input.fq.gz} | samtools view -bS -q30 > {output}.bam
sort/markdup
samtools sort {hisat.bam} -o {hisat.sort.bam} java -jar /data1/program/picard-2.23.4/picard/build/libs/picard.jar MarkDuplicates I={hisat.sort.bam} O={hisat.sort.mkdup.bam} M={hisat.mkdup.metrics}