This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem with read counts after concatenating FASTQ files

Hello ! I downloaded multiple run accession files (RNAseq-FASTQ) from ENA and tried to concatenate them using cat.

cat *_1.fastq.gz > merged_R1.fastq.gz

The output file is bigger in size, however, when I look at the read counts it is the same as the first of all the files that were concatenated. I am a beginner in NGS data analysis. Can someone please advise what am I doing wrong/missing.

Thank you

rna-seq fastq
• 1,745 views • link
written by Akshaya • Akshaya

What is ${i}? Please show the entire code including how you determined the read number.

• 0 views • link
written by ATpoint ATpoint

I used the readlength.sh script in bbmap

• 0 views • link
written by Akshaya • Akshaya

What is ${i}? Please show the entire code

• 0 views • link
written by ATpoint ATpoint 
for i in $(ls -d */)
do
    cd $i
    cat *_1.fastq.gz > ${i%?}_R1.fastq.gz
    cat *_2.fastq.gz > ${i%?}_R2.fastq.gz
    cd ..
done

I had the files for different samples in individual folders.

• 0 views • link
written by Akshaya • Akshaya

Looks fine, did checking read counts with the suggestion from RamRS help? I guess everything is fine and simply the counting went wrong.

• 0 views • link
written by ATpoint ATpoint

Yes, RamRS's suggestion did help. Thank you

• 0 views • link
written by Akshaya • Akshaya

1 answer

It could be a problem with gzip streaming library. Can you try a simple zcat merged_R1.fastq | wc -l? That should give you 4 * number_of_reads (=number of lines)

• 0 views • link
written by Ram Ram

Yes, this works. Thank you so much!

• 0 views • link
written by Akshaya • Akshaya

Log in to answer this question.

More posts like this