SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA Soft Clipping Bio.X2Y Bioinformatics 11 03-09-2015 09:08 AM
mapping parameters for SV/CNV discovery maria.b Bioinformatics 2 12-11-2012 03:16 PM
unmapped vs soft-clipping CNVboy Bioinformatics 2 04-08-2012 10:41 PM
Clipping adapters nsl RNA Sequencing 6 05-27-2011 02:08 PM
what's clipping point? jordi 454 Pyrosequencing 1 05-26-2009 01:47 AM

Reply
 
Thread Tools
Old 08-06-2013, 09:05 AM   #1
legaultm
Junior Member
 
Location: Montreal

Join Date: Aug 2013
Posts: 2
Post Clipping and trimming with MrFast (for eventual CNV discovery)

In the mrFAST manual, it says:

Quote:
NOTE: mrFAST will assume that the read lengths are uniform within a file (paired-ends may differ from each other). If the lengths of the reads are not the same in a file, it may crash or output incorrectly formatted SAM files.
Since my sequencing pipelines trims and clips the fastq files by using the fastq-mcf utility, the sequences will not be completely uniform in length. I was wondering if MrFast only worked on raw, untrimmed and unclipped data, or if this particular pre-processing step was acceptable for use with mrFAST.
legaultm is offline   Reply With Quote
Old 08-14-2013, 08:57 AM   #2
legaultm
Junior Member
 
Location: Montreal

Join Date: Aug 2013
Posts: 2
Default

So to anyone who might still be interested. After some testing, I did get a Segmentation fault when using trimmed/clipped reads. Thus, I considered two different options.

1. Use the clip/trim data, but remove additional bases at the end of reads to make the length uniform (and discard the reads that are too short for this process).

2. Use the raw reads.

I chose to use the raw reads, since for CNV calling the low quality bases at the end of the sequences shouldn't be too problematic. Also, I figured the mapping quality shouldn't be affected in a very significant manner.
legaultm is offline   Reply With Quote
Reply

Tags
clipping, mapping, mrfast, trimming, variationhunter

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 04:33 PM.


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