SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Illumina Pipeline Version t.wieland Bioinformatics 8 06-09-2011 02:03 AM
Illumina Pipeline on Macintosh hdd102 Illumina/Solexa 1 01-12-2011 03:59 AM
Bowtie to Illumina Pipeline Conversion GERALD Bioinformatics 0 06-27-2010 10:07 AM
Illumina pipeline multicoring joa_ds Bioinformatics 6 03-30-2009 11:42 AM
Help installing Illumina pipeline 1.0 koakeson Bioinformatics 1 12-03-2008 08:39 AM

Reply
 
Thread Tools
Old 11-20-2008, 06:46 AM   #1
bioinfogirl
Junior Member
 
Location: evry

Join Date: Nov 2008
Posts: 3
Question Pipeline Illumina question

Hi,

I'm a new member and I have some questions about Illumina Pipeline GAII.

In Ilumina pipeline documentation there is no documentation about script used by the pipeline. Output of the pipeline show that Goat call different scripts. Does anyone have the details of call script process? What do each script? What is output of each script?...

When I run GERALD as a standalone program with parallelization, I have warning message like 'Warning: did not find any sequences to process (raw)' or 'Warning: none of 25092 data sequences passed filter criteria '(($F[5]>=0.6))' at /env/borneo/pipeline/GAPipeline-1.0/Gerald/qualityFilter.pl line 231.' Does anyone deals with this warning? In this output the lane_tile process are not in order. Does anyone know why? Is it due to the parallelization? This is a problem for me to determine what is the lane_tile concern by this warning.

Last, I would like understand ELAND algorithm. But I didn't find any documentation. Somebody know how can I get this information? saw in documentation that ELAND make several pass alignements. I would like more explanations about the alignment process : pass, error rate...

Thanks
bioinfogirl is offline   Reply With Quote
Old 11-20-2008, 12:22 PM   #2
ScottC
Senior Member
 
Location: Monash University, Melbourne, Australia.

Join Date: Jan 2008
Posts: 246
Default

Hi,

Maybe it would be useful if you gave the commands you are using to initiate the programs.

For a description of the pipeline and how it works, including ELAND, the documents at Genographia.org (mostly written by Nava Whitford) will probably come in handy:

http://www.genographia.org/portal/to...and-algorithms

http://www.genographia.org/portal/to...eline.pdf/view

I thought there was a document on there somewhere about how ELAND works, but I don't seem to be able to find it at the moment. I'm sure someone else will post.

Cheers,

Scott.
ScottC is offline   Reply With Quote
Old 11-21-2008, 02:01 AM   #3
bioinfogirl
Junior Member
 
Location: evry

Join Date: Nov 2008
Posts: 3
Default

Thank you ScottC! genographia give a good description for image analysis but I found no algorithm description about ELAND.

"Maybe it would be useful if you gave the commands you are using to initiate the programs."

I try goat script with default parameter on a run test phiX. Our run is not good we have poor signal intensities. But I'm just starting with GAII and I don't know how deal with this error message and especially how to identify the tile and lane concerned by error in output with the fact that lane and tile are not in order in the output.
If someone deals with this situation, It would be great to share his experience.
bioinfogirl is offline   Reply With Quote
Old 11-21-2008, 09:20 AM   #4
new300
Member
 
Location: northern hemisphere

Join Date: Mar 2008
Posts: 50
Default

Quote:
Warning: none of 25092 data sequences passed filter criteria '(($F[5]>=0.6))' at /env/borneo/pipeline/GAPipeline-1.0/Gerald/qualityFilter.pl line 231
This error means that none of the reads passed "purity filtering". Purity filtering is used to discard reads which are bright in more than one channel (say A and C) and which therefore can not be base called reliably. I believe you can get the purity threshold when running Gerald, but I've not done this for sometime.

So the error is basically tell you "the analysis software doesn't think there are any good reads here so your not going to get any reads from this tile". If this is a GA2 flowcell then 25000 is an abnormally low number of reads anyway and something went wrong. You could post the intensity plots (Linked from Summary.htm, if they were produced) or take a look at the images that came off the device.
new300 is offline   Reply With Quote
Old 11-24-2008, 04:00 AM   #5
bioinfogirl
Junior Member
 
Location: evry

Join Date: Nov 2008
Posts: 3
Default

Quote:
Originally Posted by new300 View Post

So the error is basically tell you "the analysis software doesn't think there are any good reads here so your not going to get any reads from this tile". If this is a GA2 flowcell then 25000 is an abnormally low number of reads anyway and something went wrong. You could post the intensity plots (Linked from Summary.htm, if they were produced) or take a look at the images that came off the device.
Thank you new300!
Effectively the run have low intensities! This explain certainly this error message.
Do you observe desorder of tile and lane in the output pipeline?
bioinfogirl is offline   Reply With Quote
Reply

Tags
algorithms, eland, gaii, pipeline, script

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:12 PM.


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