Fastqc loop
1
0
Entering edit mode
15 months ago
camillab. ▴ 160

Hi,

I want to do bulk-RNAseq so I sent to sequence 12 samples 9pair end strategy) and I need to do the fastqc, alignment etc. I am trying at the moment do to a for loop to get the fastqc since it takes ages but I am not able to do it so I am not sure what I am doing wrong.

For one sample I would do:

/Users/camilla/Desktop/FASTQ/FastQC/fastqc /Users/camilla/Desktop/FASTQ/00_fastq/P3D40_R1_001.fastq.gz —extract -o /Users/camilla/Desktop/FASTQ/QC_before

but I am not able to make it for a loop and save the output file in the folder (QC_before).

I found this code here script

for file in *.fastq
doo
echo " fastq $file"
done

but it's not really working. Can someone help me?

Thank you and apologies if it is a stupid question!

Camilla

fastQC r RNA-seq • 838 views
ADD COMMENT
0
Entering edit mode

Avoid that "lessons" repository - it seems to be written by a bioinformatics beginner - and not maintained at all. It also has quite a few issues that I do not wish to waste my time on. Look for answers with upvotes here or on Stack Overflow.

ADD REPLY
1
Entering edit mode
ADD COMMENT
0
Entering edit mode

thank you! I don't know how did I not find it earlier!

ADD REPLY

Login before adding your answer.

Traffic: 1815 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6