Problem in execuation of filterAndTrim command in DADA2
0
0
Entering edit mode
abhisek001 • 0
@6d5973d2
Last seen 6 weeks ago
India

Hello scientific community... I was running the following command "filterAndTrim" with default parameters for trimming of 16s data in DADA2 pipeline it worked fine. The only problem I have is that it is trimming more than 95% of my reads for all samples why it is happening can anyone tell?

My Fastaq files are taken without any trimming into third-party software like cutadapt, fastp or trimmomatic. If all the reads are got trimmed how can I proceed with so little reads. Any guidance? I also tried to cut the bad sequences from the end of the sequence then most of the reads were filtered out.

out <- filterAndTrim(dataF, filt.dataF, dataR, filt.dataR, truncLen=0,
                     maxN=0, maxEE=Inf, truncQ=2, rm.phix=TRUE,
                     compress=FALSE, multithread=FALSE)

head(out) -

                    reads.in reads.out
AIJ_R1.fastq    63023      1565
AIZ_R1.fastq   104880      2593
AOJ_R1.fastq    76768      2029
AOZ_R1.fastq    75343      1866
FIJ_R1.fastq      68518      1615
Bioconductor dada2 • 406 views
ADD COMMENT
0
Entering edit mode

This problem was solved when I used trimmed sequenced with fastp and used those sequences on DADA2. Raw reads can not give the desired result. Thanks.

ADD REPLY

Login before adding your answer.

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