I have raw data fastq.gz files for some samples. Using Hisat2 for aligning reads to genome and used samtools converted sam to bam format. Also used samtools sort for sorting the bam file. For all the samples I used following command:
hisat2 -p 8 --dta --rna-strandness RF --trim3 3 -x genome_snp_tran -1 /A.1.fastq.gz -2 /A.2.fastq.gz | samtools view -Sb - > A.bam
samtools sort -T /tmp/A.sorted -o A.sorted.bam A.bam
A.1.fastq.gz 8.9G
A.2.fastq.gz 9.7G
For most of the samples this worked. For one sample I got an error at sorting step. I got the A.bam file (40G) but got an error for sorting.
[E::bgzf_flush] hwrite error (wrong size)
[E::bgzf_close] file write error
[bam_sort_core] failed to create temporary file "/tmp/A.sorted.0468.bam": No space left on device
What I have to do now?
•
written
by
Vasu
1 answer
The error you see is quite obvious: there is no disk space left on the device for the intermediate files written while sorting.
Note that you can shorten your command considerably and avoid an intermediate file by piping hisat2 directly into samtools sort:
hisat2 -p 8 --dta --rna-strandness RF --trim3 3 -x genome_snp_tran -1 /A.1.fastq.gz -2 /A.2.fastq.gz | samtools sort -o A.bam -
This probably requires a recent samtools version.
You'll additionally want to set a temporary directory with the -T option, like -T /some/place/with/space/FileName.
As you can see I added the sorting step to the end of the hisat command.
Log in to answer this question.
-
samtools index: failed to create index. ([E::hts_idx_push] NO_COOR reads not in a single block )
written by DareDevil
-
htslib/c what is the correct way to use bgzf_thread_pool ?
written by Pierre Lindenbaum
-
Adding @RG header line to bam and sam files in hisat2
written by Nemo •
-
samtools markdup
written by PROFESOR •
-
Vg gamsort temporary folder
written by arunkumarramrsh •
-
HIsat2 bam and summary output - multiple output file type
written by Payal
-
Error in extracting Counts from the bam file using htseq-count
written by Biologist
-
how can I align several fastq files with hisat2 and get htseq-count
written by Learner
-
bgzf_blocked error when using samtools
written by milanlove •
-
Error with samtools - indexing alignments: Writing to standard output failed: broken pipe - python
written by coyot001 •
