Trimmomatic command
1
0
Entering edit mode
9 weeks ago

The command i am using for trimming my adaptors from the sample data is

java -jar /data/sata_data/home/Rajdeep/Rajdeep/Trimmomatic-0.39/trimmomatic-0.39.jar PE -phred33 CC1_R1.fastq.gz CC1_R2.gz CC1_R1_paired.fastq.gz CC1_R1_unpaired.fastq.gz CC1_R2_paired.fastq.gz CC1_R2_unpaired.fastq.gz ILLUMINACLIP:/data/sata_data/home/rah=jdeep/Rajdeep/Trimmomatic-0.39/adaptors/TruSeq3-PE.fa:2:30:10 MINLEN:25

Is it correct to use this command before aligning with STAR or Hisat2?

STAR RNASeq Trimmomatic • 621 views
ADD COMMENT
0
Entering edit mode

Please provide more information. You've shown a command to trim PE data, but haven't said what you're trying to achieve. But yes, trimming raw reads before alignment is generally a good idea, but it depends on what the QC looks like, which you've also not provided.

ADD REPLY
1
Entering edit mode
9 weeks ago

Keep in mind for RNASeq

1) if the library was done properly, it's very unlikely you have many fragments shorter than your read length. You probably won't see much adapter at all. 2) splice-aware aligners don't mind soft trimming to make an alignment. Even if you had adapters on a read, it will probably still align to the right place, and be counted as belonging to the correct gene.

Now, if you were doing DNASeq, or do novo assembly, then removing adapters would be more important. But your typical RNASeq experiment, where all you want are gene/transcript counts? Don't agonize. It's not going to be much a problem, even if you did nothing.

ADD COMMENT

Login before adding your answer.

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