How to obtain unmapped reads from the paired end samples using hisat2
Am having trouble extracting unmapped reads using hisat2. I have 15 samples that I need to keep unmapped reads for a different analysis.
When I use:
hisat2 -x ./HISAT2/Index/Gff-Ind -p 8 -1 ./data/SN1_R1_trimmed.fastq.gz -2 ./data/SN1_R2_trimmed.fastq.gz -S ./data/HISAT2_unmapped_keep/SN1.sam --un-conc ./data/HISAT2_unmapped_keep/unmapp
It works okay but I get a big SN1.sam files and two unmapped reads (unmapped.1 and umapped.2)
When I try to modify the script for obtain the small SN1.bam file from sam file and rename unmapped.1 & .2 file to SN1_unmapped.fastq I get exit code 1 error which have tried to solve in vain.
hisat2 -x ./HISAT2/Index/Gff-Ind -p 8 -1 ./data/SN1_R1_trimmed.fastq.gz -2 ./data/SN1_R2_trimmed.fastq.gz -S ./data/HISAT2_unmapped_keep/SN1.sam -> SN1.bam --un-conc ./data/HISAT2_unmapped_keep/SN1_unmapped.fastq
I'd like some help please.
• 975 views
•
link
updated
by
Ram
•
written
by
otieno43 •
•
written
by
otieno43 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Hisat2 index and alignment question
written by Athena •
-
Hisat2 Error
written by Forough •
-
extract unmapped reads from the paired end samples using hisat2
written by winmoon6 •
-
How to plot alignment statistics; how many reads are mapped to the genome?
written by anamaria
-
bowtie2 --un vs --un-conc
written by blur
-
Unmapped reads to fastq (hisat2 / samtools / bedtools) not matching pair-end reads
written by juan.crescente •
-
denovo assembly of unmapped reads
written by ashish
-
Python scripts to map multiple pair end files to the genome
written by Bioinfonext
-
How to get percentage of reads mapped to each Reference Genome
written by bandanaschapagain •
-
Bowtie2 --un-conc output empy
written by t2
You can't can do this sort of thing on a program command line (unless you were trying to conceptually show what you wanted to do). You will need to pipe the SAM file to
samtoolsto do the conversion. See HISAT2 output direct to bam