ViReMa not working after adapter trimming
0
0
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 ViReMa.py --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
ADD COMMENT

Login before adding your answer.

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