I'm new at this, and I'm trying to do RNASeq. I've done the quality analysis already with fastqc. Now, I'm on the quality control part. I need to use seqtk trimfq or fastx_trimmer to trim the first 5 bases.
The following commands are not working. The seqtk command produces a single output file. The fastx_trimmer command says "writing nucleotides failed: Disk quota exceeded".
I'm doing this on a terminal on OSX to login in to the machine at school.
find *fastq | parallel -j 4 seqtk trimfq -b 5 {} > trimmed_{}
find *fastq | parallel -j 4 --no-notice fastx_trimmer -Q33 -f 6 -i {} -o trimmed_{}
Any suggestions are greatly appreciated! I'd like to do it in parallel since I have approximately 15 .fastq files. After this step I intend to build a hisat index.
0 answers
No answers yet.
Log in to answer this question.
-
How Should I Trim My FASTQ File Now That I Merged Paired-End Reads
written by isu2017 •
-
Extract Sequences From [Gtf File] + [Genome Fasta File], combine matching XLOCs
written by Tom Koch
-
macs2 "Effective genome size" for repetitive genomes
written by Menachem Sklarz •
-
Jupyter and Anaconda not downloading for me
written by chrisclarkson100
-
Running BWA parallel on 1000 samples
written by ebrown1955
-
Converting .narrowPeak file to GRanges for use in Bioconductor ChIPseeker
written by dally
-
Filtering sanger_fastq, illumina_fastq files using Fastx-Toolkit
written by bioinfo
-
Pre-processing And Quality Control of the Raw NGS Data -- > trimming
written by David_emir
-
R1 R2 orientation (RNA-seq)
written by biotech
-
Fast search for many patterns
written by sbliven •
"writing nucleotides failed: Disk quota exceeded" Are you sure you have enough space on disk for the output files? Since you say you are using a computer at school it might well be that your user quota is not that big.