HISAT2 output direct to bam sorted
1
0
Entering edit mode
22 months ago
sansan96 ▴ 130

Hi all, I am mapping my samples to a reference genome with HISAT2 and I was wondering if it is possible to get the outputs in bam format and also if it is possible to sort by coordinates directly?

Could you help me modify my command?

hisat2 --dta -x genome -1 ../fastq_zea_mays/SRR214880_1.fastq -2 ../fastq_zea_mays/SRR214880_2.fastq -S SRR214880.sam
sorted HISAT2 bam to output direct • 1.5k views
ADD COMMENT
5
Entering edit mode
22 months ago
ATpoint 85k
hisat2 --dta -x genome -1 ../fastq_zea_mays/SRR214880_1.fastq -2 ../fastq_zea_mays/SRR214880_2.fastq \
| samtools sort -o sorted.bam
ADD COMMENT
0
Entering edit mode

Thanks so much, what is the function of "\" ?

ADD REPLY
1
Entering edit mode

It's just a line break so I can write code into two lines for readability and structure. It's not necessary for function.

ADD REPLY
0
Entering edit mode

thanks very much

ADD REPLY

Login before adding your answer.

Traffic: 2435 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6