SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Error found while mapping Illumina data using BWA nirav99 Bioinformatics 15 01-14-2013 07:30 PM
Compiling simNGS fails gene coder Illumina/Solexa 1 12-15-2011 05:53 AM
Tophat SOLID PE data not mapped in pairs DerSeb Bioinformatics 4 11-04-2011 11:12 AM
where can i found chip-seq data of Drosophila hannat Epigenetics 3 09-01-2010 11:03 AM
where can i found raw data generated by solexa hannat General 4 01-12-2009 07:56 AM

Reply
 
Thread Tools
Old 07-21-2011, 09:13 AM   #1
oiiio
Senior Member
 
Location: USA

Join Date: Jan 2011
Posts: 105
Default Novoalign: No pairs found with simNGS data

I generated some paired end data using simNGS and simLibrary, as they seem to be very good for illumina HiSeq.

My problem is that novoalign does not recognize the pairs, and returns "No pairs found". It is strange because when mapping my real illumina data, it performs excellently. Any suggestions would be very helpful.

Before suggestions to use different mappers are made, I already tried BWA (which worked fine) and Stampy (same problem as novoalign). Are there known issues with simNGS and mapper compatibility?
oiiio is offline   Reply With Quote
Old 07-21-2011, 07:21 PM   #2
zee
NGS specialist
 
Location: Malaysia

Join Date: Apr 2008
Posts: 249
Default

That is quite odd.

Are you using FASTQ simulated reads. Novoalign needs to have read1 and read2 of a pair in separate files.
Also check the coding of your intensities and use the correct coding for Illumina/Sanger FASTQ.
zee is offline   Reply With Quote
Old 10-03-2011, 06:20 PM   #3
sparks
Senior Member
 
Location: Kuala Lumpur, Malaysia

Join Date: Mar 2008
Posts: 126
Default

Hi,

Novoalign and Stampy both expect reads of a pair to be orientated so that one reads is on +ve strand, the other on -ve strand with 3' ends of reads towards the centre.

5---------------->3 3<-------------------5

This is how paired reads come off the Illumina. BWA doesn't apply this rule and will make proper pairs from any two adjacent mappings and is one of te reasons Novoalign and Stampy are more accurate than BWA.

My guess is that NGSsim isn't producing reads in the correct orientation.
The other problem could be the insert size produced, Novoalign and Stampy both default to around 250bp with std dev of 50bp, if NGSsim produced pairs with mean of much larger than this, say >600, then it's unlikely Stampy or Novoalign would map any pairs.
Colin
sparks is offline   Reply With Quote
Old 11-20-2012, 10:59 AM   #4
dumps
Junior Member
 
Location: Heidelberg

Join Date: Nov 2012
Posts: 1
Default

Hi,

be aware that the --paired option of simLibrary DEACTIVATES generation of paired end reads - Don't ask me why. But maybe your problem is that simple...
dumps is offline   Reply With Quote
Old 11-20-2012, 11:55 AM   #5
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

Please give us an example of the first 12 lines of each fastq pair file.

Please note the " --hdrhd off " command line parameter available with novoalign. I think this turns off checking that names in pairs match.

Do make sure you have matching number of pair elements in fastq1 and fastq2 files.

Last edited by Richard Finney; 11-20-2012 at 11:58 AM.
Richard Finney 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 07:34 AM.


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