Hi Elsayed

At the risk of overwhelming you with examples:

The Genomics Virtual Lab also provides basic and advanced RNA-seq and Variant Detection tutorials; with a (small) public Galaxy instance specifically aimed at running these two tutorials.

Ron Horst

Project Manager, Genomics Virtual Lab

Level 6, QBP, The University of Queensland, QLD 4072

t. 07 3346 2276 | m. 0417 538 723 | e. r.horst@uq.edu.au | w. www.genome.edu.au

From: galaxy-dev-bounces@lists.bx.psu.edu [mailto:galaxy-dev-bounces@lists.bx.psu.edu] On Behalf Of Jennifer Jackson
Sent: Saturday, 4 January 2014 6:59 AM
To: Elsayed Hejazy
Cc: Galaxy Dev
Subject: Re: [galaxy-dev] Help Needed

Hello Elsayed,

Your protocol below seems to be a mix of a variant detection and an RNA-seq workflow. To build a workflow for RNA-seq, you will want to compare your steps with the protocols in the link that I sent you.

If you want more examples, many more can be found here (for both RNA-seq and variant analysis, these are distinct analysis). Some have workflows that you can import and use directly or modify.

https://wiki.galaxyproject.org/Learn#Other_Tutorials

Shared Data -> Published Pages

In particular, there is a sample workflow in this tutorial that will help you to understand what the different steps are for and what order they go in.
https://usegalaxy.org/u/jeremy/p/galaxy-rna-seq-analysis-exercise

Best,

Jen
Galaxy team

On 1/3/14 10:57 AM, Elsayed Hejazy wrote:

Dear Dr. Jennifer,

A lot of thanks for your reply its really mean alot for me.

i am still NGS junior i trying to do the following steps kindly give me the right order for these procedures

Data Description: i have to samples each sample consists of 14 FASTQ file (7 forward and 7 reverse ) i think this mean its paired end from Illumina then i will try the following workflow to got best results

1- Drag tow input dataset into workflow one for forward sequences file and one for reverse to use paired end option in TOPHAT tool later and when i run this workflow i will select multiple selection for the 7 forward files to analyse all of them at the same time

2- Drag FASTQC and link with last step for each to got if these file may be illmina 1.8 version or older.

3- Drag FASTQ Groomer and link with last step if files older than 1.8 version to prepare as .FASTQSANGER format.

4- Drag Filter FASTQ and link with last step to remove redundancy of sequences.

5- Drag FASTQ trimmer to remove unwanted ends of sequenced may occur

6- Drag Manipulate FASTQ and link with last step (i dont know why).

the above six steps done twice to generate to files as output to make as input for the following steps.

7- Drag TOPHAT for illumina and make it accept paired end files and link each file generated from QC to TOPHAT this step used to align and map with reference genome.

8- Drag Cufflinks and link with aligned BAM file generated from TOPHAT to create an assembly

9- Drag Cuffmerge and link with GTF file from Cufflinks this step to merge all assemblies generated in Cufflinks

10- Drag Cuffcompare and link with last step to got detailed reports for accuracy of all generated assemblies.

11- Drag MPileup and link with TOPHAT BAM file to generate file containing SNPs sites.

12- Drag Pileup-to-Interval and link with MPileup step to filter the number of output SNPs to successive one or the most accurate. (i dont know what is the difference between this tool and Filter Pileup).

- i dont know what is the tools used to know the copy number variation CNV

i need to know how to separate human sequences from sample may infected with any other sequences (is this at alignment stage)

i need to know the perfect order of steps if this order is not completely right.

Is this what i should do to make a good NGS workflow got all possible information form dataset

Really i am so so so sorry for disturbance - waiting your reply

Best Regards,

elsayed
On Thu, Jan 2, 2014 at 6:05 PM, Jennifer Jackson <jen@bx.psu.edu> wrote:
Hello Elsayed,

Protocol help for RNA-seq analysis can be found here:
https://wiki.galaxyproject.org/Support#Tools_on_the_Main_server:_RNA-seq

The QA/QC steps should be done before mapping, on individual datasets (such as replicates) or on partial or merged datasets as needed (if that is how the data was sequenced, just be consistent). Only keep in mind that the larger the dataset, the more compute some of these steps can require.

Hopefully this helps,

Jen
Galaxy team

On 1/1/14 1:32 AM, Elsayed Hejazy wrote:
i need to know more about the order of steps for RNA Seq data analysis

Is alignment and assembly should done first to combine all FASTQ reads into one file and start analysis or i should do Quality Control ,filtering , trimming and manipulation first and do alignment and assembly at the end of analysis to use tophat for example or any other analysis tools in Galaxy

Thank you

elsayed
___________________________________________________________
Please keep all replies on the list by using "reply all"
in your mail client.  To manage your subscriptions to this
and other Galaxy lists, please use the interface at:
  http://lists.bx.psu.edu/
 
To search Galaxy mailing lists use the unified search at:
  http://galaxyproject.org/search/mailinglists/
-- 
Jennifer Hillman-Jackson
http://galaxyproject.org

--

Jennifer Hillman-Jackson

http://galaxyproject.org