ViReMa not working after adapter trimming
Entering edit mode
5 days ago
berelly • 0

Hi, I have a viral RNA-seq dataset that I am trying to run through ViReMa to look for deletion junctions. When I input my fastq files directly into ViReMa without trimming adapters first, my results look about how I would expect (multiple junctions returned with several reads each). However, when I trim adapters with cutadapt prior to running ViReMa, I get very few junctions returned, or (usually) no junctions at all.

Is it possible that adapter contamination would be causing artifactual junctions to be called? Or am I doing something super wrong in my trimming step and somehow removing information that ViReMa needs to proceed?

My cutadapt parameters look like this (just trimming Illumina universal adapters):

cutadapt -a AGATCGGAAGAG -o sample_cutadapt.fastq sample.fastq

And my ViReMa parameters looks like this:

python --MicroInDel_Length 20 -DeDup --Defuzz 3 --Seed 25 --N 1 --X 2 --Output_Tag sample -ReadNamesEntry --p 6 bowtie1index/virus sample_cutadapt.fastq sample.results 

I'm new to this type of analysis, so any help is appreciated.

cutadapt trimmomatic ViReMa trimming adapter • 124 views

Login before adding your answer.

Traffic: 1981 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6