Question: STAR Alignment Score reduced after trimming
0
gravatar for dec986
3 months ago by
dec98690
United States
dec98690 wrote:

Hello,

I'm comparing three different QC tools (trimgalore, trimmomatic, & bbduk) and then seeing how STAR alignments are affected by all of them, comparing against a noQC run. Specifically, I am comparing the "AS" and "nM" fields in the BAM file. I figured that the alignment scores should increase, but the QC file without any QC gets the best alignment scores, albeit with higher mismatch rates.

I am concerned that the alignments after QC may be worse because of this, specifically the means of the alignment scores can be seen:

trimgalore  - noQC  = -5.57
trimmomatic - noQC  = -2.99
bbduk       - noQC  = -1.15

Why are the QC tools causing worse alignment scores in STAR?

rna-seq star • 199 views
ADD COMMENTlink modified 3 months ago by Devon Ryan79k • written 3 months ago by dec98690
4
gravatar for Devon Ryan
3 months ago by
Devon Ryan79k
Freiburg, Germany
Devon Ryan79k wrote:

There's little point in trimming reads with STAR (or other tools that do local alignment), at most trim adapters and very low quality bases (phred scores up to ~3 or so). Presuming you did quality trimming, it's likely that the score you used was much too high, which decreased the number of bases that STAR would have otherwise aligned and therefore decreased the alignment scores. Even adapter trimming isn't perfect, since inevitably you'll end up trimming off a bit of sequence. For both our command line and Galaxy RNAseq workflows, we skip trimming entirely and use STAR with the raw reads. That tends to produce quite nice results.

ADD COMMENTlink modified 3 months ago • written 3 months ago by Devon Ryan79k

thank you! who is "our"?

ADD REPLYlink written 3 months ago by dec98690
1

who is "our"?

@Devon is currently a bioinformatician/data manager at the Max Planck Institute for Immunobiology and Epigenetics in Freiburg, Germany.

ADD REPLYlink written 3 months ago by genomax48k

great, thanks, I have to let my colleagues know. We're comparing several different RNA-Seq pipelines.

ADD REPLYlink written 3 months ago by dec98690

Evidently a great manager too

ADD REPLYlink written 3 months ago by Kevin Blighe19k
1

Thankfully I manage the data and not the people :)

ADD REPLYlink written 3 months ago by Devon Ryan79k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1144 users visited in the last hour