What Is The Best Way To Get Rid Of Low Quality Bases In A Sequence
3
0
Entering edit mode
7.9 years ago
kmkdesilva ▴ 90

Hi all,

I would like to know what is the best way to get rid of low quality bases in a sequence using fastx tool kit. Is quality filtering the sequence using a given phred score or trimming the sequence.

By quality filtering we might loose bases in the middle of the sequence. From trimming we can chop of the ends. So I would like to know what is the best way.

Regards, kmkdesilva

fastx trimming • 3.5k views
ADD COMMENT
0
Entering edit mode
7.9 years ago
IV ★ 1.3k

Quality filtering by use of fastq_quality_filter does not clip bases within the sequence but rather rejects reads that do not meet the user defined criteria.

For instance you would want to keep reads that have more than q quality for at least p percentage of the read length.

fastq_quality_filter -q min_quality -p min_percent -i file -o outfile

You could select to trim when low quality base calls are at the last few nucleotides of your reads and to quality filter when there are more generalized quality issues.

In general, fastx toolkit tools are considered quite strict and you end up losing quite a lot of your sequences.

Other tools for trimming, adapter clipping and/or quality filtering are more gentle (such as scythe, cutadapt, etc).

Quality trimming and filtering should be performed with care, in order to actually increase the quality and quanity of what you get in the end.

Adapter removal to my opinion should be performed more vigorously.

Cheers,

IV

ADD COMMENT
0
Entering edit mode
7.9 years ago
Haluk ▴ 190

I used the PoPoolation before, and it works well.

ADD COMMENT
0
Entering edit mode
7.9 years ago
rob234king ▴ 600

Primers and adaptors need to be watched out for, especially 454 data but the best trimming tends to be done by the software in question, atleast for mapping. For instance I dont trim using third party software with novoalign because it has a clipping function built in which I can use and when tested with my trimming it performs significantly better. Or sometimes even with all the junk will map better rather than trimming and just filter after. For assembly I found Quake good at automatically cleaning up reads and did show benefit. If I have to trim then trimmomatic is my favourite.

ADD COMMENT

Login before adding your answer.

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