FeatureCounts Number of Reads
0
0
Entering edit mode
4 months ago
santos48 ▴ 30

Hi everyone, I have over 4000 human Rna-seq data reading by nanopore Minion, I merged all files with the 'cat' command. I am using featureCounts for counting features of aligned Rna-Seq. After the featureCounts when I control with multiQC report I saw the number of reads 4k but it should be over the million.

Should I merge after aligned against the ref genome before start analysis?

RnaSeq FeatureCounts • 387 views
ADD COMMENT
0
Entering edit mode

I have over 4000 human Rna-seq data reading by nanopore Minion

What does that mean? You have 4000 files, or samples? Please explain.

ADD REPLY
0
Entering edit mode

4000 fastq files

ADD REPLY
0
Entering edit mode

How did you merge them with cat? Can you show your command?

ADD REPLY
0
Entering edit mode
cat /*.fastq >all.fastq
ADD REPLY
0
Entering edit mode

If this is the precise command then this has two problems:

a) / means the root folder where most likely these files are not stored, you probably meant to use ./ which means current directory, and

b) it might be that the newly created file gets concatenated to itself, better:

find . -maxdepth 1 -name "*.fastq" | xargs cat > new.fastq
ADD REPLY

Login before adding your answer.

Traffic: 2091 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