How to Write Loop For Trimming Every FastQ File Individually
I have human Rna-Seq data by reading nanopore. Data has over 1000 .fastq files. I want to trim all fastq files individually. I am using NanoFilt for trimming. All Fastq files are in the directory(not gzip format).
How can I write a loop or command line for trimming all fastq files individually?
I am using Linux Ubuntu 21.04.
• 4,967 views
•
link
updated
by
Carlo Yague
•
written
by
santos48 •
•
written
by
santos48 •
1 answer
Hello, this is a rather basic question. You could start by googling something like "bash loop files directory" to get plenty of examples on how to write loops in bash. The basic structure of the loop is probably going to be something like this:
for file in $directory; do
echo $file
<nanofilt cmd> $file
done
Then if you run into specific issues, feel free to come back and ask for more direction.
• 0 views
•
link
written
by
Carlo Yague
Log in to answer this question.
More posts like this
-
How to Create Loop for all fastq file
written by santos48 •
-
Splitting Big Fastq directory into smaller parts
written by santos48 •
-
FeatureCounts Number of Reads
written by santos48 •
-
Look at Different Locus in the Bam File
written by santos48 •
-
Create Annotation for Long Read Nanopore Rna-Seq Data
written by santos48 •
-
Nanopore Preprocessing Step
written by santos48 •
-
need help with script for hisat2 alignment
written by Antonio.Aubry •
-
Using trimmomatic on multiple single-end read files
written by armel.donkpegan •
-
Linux for loop for 2 inputs and 4 outputs for Trimmomatic fastq quality trimming
written by Rashedul Islam
-
Sra-Toolkit 2.3.1 Fastq-Dump Multiplefiles At Once Windows Command Line
written by deepbiofever •
Can you explain how to run with porechop with all fastq file in the loop ? .
@carlo has explained the basics of a shell loop below. You can try to figure out how to replace
porechopin place ofnanofiltin the loop below. You should take a look atbasenamewhich will allow you to extract just name of sample (here is an example for you to try things out): How to modify my sbatch script to make it be able to be used for other filesWrite a loop. Use
echoto test things out and ask questions with specific examples.I guess this will also work for the same:
I wrote loop like
for file in *.fastq do porechop -i ${file} > /Diroctory/${file} > donebut how can I specify for when all data is trimmed stop the loopIt should not be necessary to stop the loop since the above command will consider each
.fastqfile only once.Thank you for answer. How can I add another
do <something>donein the same bash for example I want to use some filtering tools after the trimming.Please use google to learn shell loops. devarora and Carlo have helped you get started, we cannot feed you the directions for each step.
Just add more lines of code commands between
doanddone. In my answer bellow, the code does two things for each file: it prints the file name and then it runs thenanofitcommand.