Use a command on multiple files

Hi,

I want to give two files as input in the following command using loop and generate two respective files as output

File 1: File1_1.fastq.gz
File 2: File1_2.fastq.gz

bbduk.sh -Xmx1g in1= in2= out1= out2= ref=adapters.fa ktrim=r

Thanks

bash bbtools bbduk.sh

• 1,044 views • link

updated 3.4 years ago by Ram

• written 4.9 years ago by utkarsh.sood •

What did you try? https://www.google.com/search?q=loop+multiple+variables+bash

• 0 views • link

updated 3.4 years ago by Ram

• written 4.9 years ago by Pierre Lindenbaum

0 answers

No answers yet.

More posts like this

Discrepancy in total number of bases in trimmed read1 and read2 files after BBDuk
written 2.7 years ago by CTLong •

Hi all, After performing adapter trimming with `bbduk.sh`, I found that the total number of bases in the read1 file is different compared with the …
How to trim bases with bbduk.sh
written 2.8 years ago by shivangi.agarwal800

Hi, I am using BBMAP following command to trim adapter sequences from the fastq files. ``` bbduk.sh -Xmx1g \ in1=1_R1_001.fastq.gz \ in2=1_R2_001.fastq.gz \ out1=1_R1_001-trimmed.fastq.gz \ …
BBDuk problem: maskMiddle was disabled because useShortKmers=true
written 3.6 years ago by raniai •

Hi, I'm trying to trim adapter from Illumina metagenomics reads using bbduk.sh. The code worked well with some samples, but then it gave me an …
Reading in input file into a Snakemake pipeline
written 4.0 years ago by Timotheus •

Hello, I'm trying to build my first Snakemake worklflow. I've got four sets of short reads for two samples (left/right reads for samples HB001_naive/coevolved) that …
BBDuk returned java.lang.AssertionError
written 4.3 years ago by melissachua90 •

I downloaded all the experiments for `SRP332503` esearch -db sra -query SRP332503 | efetch -format runinfo | cut -d ',' -f 1 | grep SRR …
BBDuk quality filtering not producing expected results
written 4.8 years ago by rebeliscu •

I'm trying to trim/filter low quality reads from paired-end exome-seq data, using BBDuk. I used the command: ``` for ea in $files; do R1="$ea" R2=$(echo …
Bash loop for bbduk
written 5.4 years ago by avelarbio46 •

Hello Everyone! I'm trying to run a loop to trim my rnaseq reads in bbduk. However, I can't seem to input files from a specific …
bbduk.sh exceeding requested threads in an alignment pipeline
written 5.9 years ago by ido.idobar •

Hi all, I'm running a pipeline that includes `bbduk.sh` for reads trimming, then aligning the trimmed reads to a reference genome using `bowtie2`, piping the …
BBDuk does not remove all adapters
written 6.9 years ago by f.dinesh •

Hey all, I'm just getting into rna-seq analysis and have run into a problem that's is very confusing. I ran multiplexed RNAseq samples on a …
solving heap size problem in linux
written 7.4 years ago by Liftedkris •

hi. i am trying to remove contaminants from my fastq files using bbduk but it gives the following error. it says invalid maximum heapn size. …