SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Low quality in the middle of the read? MG1655 Illumina/Solexa 0 08-12-2012 12:22 PM
Hello from the middle of the USA smcinturf Introductions 0 07-25-2012 12:01 PM
TopHat doesnt handle 450million reads EBER Bioinformatics 2 05-25-2012 05:09 AM
Bottom Middle Swath and other focus issues csquared Illumina/Solexa 21 03-02-2012 05:19 AM
how does bowtie handle n's in reads? ref sequence? nashp Bioinformatics 0 01-10-2011 04:09 PM

Reply
 
Thread Tools
Old 01-07-2014, 11:07 PM   #1
themysticgeek
Junior Member
 
Location: India

Join Date: Dec 2013
Posts: 4
Question How to handle Ns in the middle of reads

For my illumina data fastqc shows presence of N's at positions 13,14,15 in 101 bp longs reads. If i go for cropping first 15 bases by using trimmomatic, it solves the problem but i lose a lot of data. I wanted to know that if i retain the N's what sort of problems would they cause during alignment(bwa+stampy)/variant calling(unified genotyper) and how can i handle these problems?

If any body faced a similar problem how did you handle it? Similar questions asked on different forums but none answered. Could not find a resourse on how variant calling programs handle N's. Do they ignore them? Or consider them as a variation with low confidence scores?

Last edited by themysticgeek; 01-08-2014 at 01:55 AM.
themysticgeek is offline   Reply With Quote
Old 01-08-2014, 03:04 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,978
Default

Are there N's at those positions in *all* reads? That would almost certainly indicate a technical problem of some kind with this run. In general your sequence provider should not have released this data if that is the case.
GenoMax is offline   Reply With Quote
Old 01-08-2014, 04:42 AM   #3
themysticgeek
Junior Member
 
Location: India

Join Date: Dec 2013
Posts: 4
Default

The N's are in ~50% of the reads. I have attached the Fasqc image for per base n content. This is particular to this sequening run. Did not observe this problem in the other runs

Last edited by themysticgeek; 01-08-2014 at 04:45 AM.
themysticgeek is offline   Reply With Quote
Old 01-20-2014, 04:25 AM   #4
benjaminsb
Junior Member
 
Location: Oxford

Join Date: Aug 2013
Posts: 2
Default Similar problem

I've just run FastQC on a published RNAseq dataset (SRX294957) and I see a very similar pattern:



where on position 21-23, 80% of reads are N's. As I'm only interested in expression, I could accept the low read quality as long as the aligner accepts it. I'm using RSEM+bowtie for the purpose, so I'm wondering if bowtie will match NNN against anything in the reference?
benjaminsb is offline   Reply With Quote
Old 01-20-2014, 04:27 AM   #5
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

It will. You can alter the mismatch score due to an N and the maximum number of allowed Ns if you need to (you'll likely need to tweak the minimum allowable score if you do so).
dpryan is offline   Reply With Quote
Reply

Tags
ambiguous bases, fastqc

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 02:08 AM.


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