SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
filtering out low quality sequences from Sanger Sequencing MAmineHassani Bioinformatics 1 05-03-2018 04:47 AM
All sequence bases have the same quality score. xiangwulu Bioinformatics 6 06-24-2016 08:45 AM
Illumina RNA-seq: should I trim low quality bases prior to mapping? JonB RNA Sequencing 1 02-14-2013 04:49 AM
Periodical illumina read length distribution after trimming of low-quality bases luxmare General 4 12-20-2010 04:18 PM
Quality trimmming / Mask low quality bases? bbimber Bioinformatics 9 03-25-2010 02:40 PM

Reply
 
Thread Tools
Old 10-17-2018, 11:30 PM   #1
MNMoller
Junior Member
 
Location: Denmark

Join Date: Oct 2018
Posts: 3
Default Filtering low quality bases

Hi,

I would like to filter out or mask bases in my SAM file that has a quality value below a defined number. Is there a way to do this?

Alternatively I want to filter out the bases or mask them in my FASTQ file, but this has to be done without disturbing the BWA-MEM alignment, which I intend to do afterwards. This could possibly be done by replacing the law quality base call with a "N".

I am aware of quality trimming and removing reads by MAPQ, but this will not remove low quality bases within the read.

I would appreciate any input!
MNMoller is offline   Reply With Quote
Old 10-18-2018, 10:09 PM   #2
finswimmer
Member
 
Location: Europe

Join Date: Oct 2016
Posts: 55
Default

Hello MNMoller,

could you please describe why you want to do this?

fin swimmer
finswimmer is offline   Reply With Quote
Old 10-18-2018, 10:24 PM   #3
MNMoller
Junior Member
 
Location: Denmark

Join Date: Oct 2018
Posts: 3
Default

Yes, and this may require that I explain a little bit about my setup.

I am working on determining the fidelity of different enzymes. Therefore I have amplified a 200 bp target in a PCR and then sequenced it and counted the mistakes. I want to count the mistakes made by the enzyme and not the method. If I can get rid of the low quality base calls then I am also more sure that the basecall made by the instrument is correct. I have tried to overcome this by measuring background of the sequencing, which in theory should be sufficient.
However when quality trimming my read (from both ends), my results start to make more sense. So I guess the bottom line is that I want to test if a "quality trimming" within the sequence, will help the results make even more sense.
MNMoller is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:09 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO