Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
This is a sticky topic.
X
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • #16
    Trimming PE reads in Galaxy

    @weijenc

    My suggestion for trimming paired-end reads in Galaxy is:

    (1) Join them using the FASTQ joiner;
    (2) Filter them using the Filter FASTQ tool;
    (3) Split them using the FASTQ splitter.

    Best,
    J.

    Comment


    • #17
      Problem in grooming

      Thanks for the reply to my previous post.

      I have been trying to work with the paired-end dataset (SRR131208, two files). After grooming (solexa to fastq sanger), however, quality values are between 5 and 0. Did I do something wrong?

      Thanks,


      WJ

      Comment


      • #18
        Your data is almost certainly not solexa format; most newer Illumina data is already fastqsanger, in which case the groomer is not needed.

        See the Wikipedia entry for FASTQ for more details:



        You should be able to look at the first few reads of your datasets to determine the FASTQ format.

        Best,
        J.

        Comment


        • #19
          Hi all,
          I am new in the forum and also in the RNA seq analysis field.
          I just a got the results of my RNA sequencing and I am trying to map my reads using tophat on galaxy. I am working with Windows 7. After the FastQC analysis, I converted my reads with FASTQ Groomer and then I run tophat. It has been 2 days and the process hasn't started yet.

          Does tophat (in galaxy) run on windows7?
          Usually how long does it take a mapping analysis (about 20 Mb the size of the genome reference)?

          Thanks

          Comment


          • #20
            @Peppe

            I assume you're using the public Galaxy server at https://main.g2.bx.psu.edu/ , yes?

            If so:

            (a) Galaxy will work fine on Windows, though you'll want to you Firefox or Chrome as your Web browser going forward so that you can use all of Galaxy's functionality.

            (b) The server is very busy right now, so it may take a couple days for your job to start. Do not restart the job or it will go to the end of the wait list. Once your job starts, it should go quickly (4-8 hours is a good estimate) because your genome is small.

            Best,
            J.

            Comment


            • #21
              Thanks jgoecks,
              Yes, I am using the public Galaxy server.
              As you said the job took a while to start, but after that it was quick enough. Now there is an update of the server so I have to wait. After that I will go ahead with Cufflinks,
              Cuffcompare, Cuffmerge and Cuffdiff. Still not sure what it comes out and how to handle the data, but I will figure out. This is the first time that I am going through RNAseq analysis and having a forum like this is very helpful.
              Thanks a lot,
              Peppe

              Comment


              • #22
                Hello,
                I tried to run cufflinks after TopHat mapping and as a genome reference annotation it requires a file in gtf format. I downloaded the genome annotation from the JGI website (I am working with a fungus whose genome was sequenced in 2010). The file annotation that I have is in gtf, but when I extract it it becomes gff, a format not supported by cufflinks in Galaxy. I converted my file with the workflow found at https://main.g2.bx.psu.edu/u/jeremy/...with-cufflinks. The generated genome annotation is in the gtf format (the only difference that I see is in the first column, scaffold become chrscaffold), but when I run cufflinks all the output were 0. Now I am running again cufflinks without using the annotated genome to see if there is any difference in the results. I was wondering if the conversion that I did is correct, if anyone knows if there is another way to do it, or if I should email JGI people to see if they have a GTF format for the genome annotation I need. I don't think I can go ahead with the analysis without a genome annotation in GTF.
                Thanks

                Comment


                • #23
                  As you guys can see in the previous messages, I am new in this field. I am using the public Galaxy server, and for the first time I run sample with TopHat-Cufflink-Cuffcompare-Cuffdiff. Can anyone tell me how to get the data generated by Cuffdiff? Is there any tool in Galaxy for visualizing them? I am interested in gene/transcript differential expression testing. Should I select only the genes with significance YES?
                  Thank you

                  Comment


                  • #24
                    @Peppe

                    Cummerbund is the preferred option for visualizing Cuffdiff output, but it's not yet integrated into Galaxy.

                    However, there is a scatterplot that's handy for visualizing FPKM data for different conditions: click on the visualize icon at the bottom of a Cuffdiff expression dataset and plot column 10 (sample 1 FPKM) vs. column 14 (sample 2 FPKM). And, yes, you may want to filter for only significantly expressed genes/transcripts (using the Filtering tool).

                    Finally, these questions are best asked on the galaxy-user mailing list, where you can communicate with the Galaxy user community: http://user.list.galaxyproject.org/

                    Good luck,
                    J.

                    Comment


                    • #25
                      @Peppe re: JGI gene annotation

                      No need to run that workflow on the annotation in your case. Just make sure that your gene annotation's contig names (column 1) match the names in your genome's FASTA file.

                      J.

                      Comment


                      • #26
                        Thanks Jeremy.
                        I will use the Galaxy user community: http://user.list.galaxyproject.org/

                        Comment


                        • #27
                          I posted a partial solution on iplant's community forum to concerning on how to view Cuffdiff output from Galaxy to cummerbund: http://ask.iplantcollaborative.org/q...-gff3gtf-file/

                          I hope this helps everyone in the future.

                          Comment


                          • #28
                            Hello everyone, my name is andrean, i'm starting to do analysis my RNA sequencing, i want to make a heatmap, and GO analysis, could you tell me how to do using R or another open source program, thank you

                            Comment

                            Latest Articles

                            Collapse

                            • seqadmin
                              Non-Coding RNA Research and Technologies
                              by seqadmin




                              Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

                              Nobel Prize for MicroRNA Discovery
                              This week,...
                              10-07-2024, 08:07 AM
                            • seqadmin
                              Recent Developments in Metagenomics
                              by seqadmin





                              Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                              09-23-2024, 06:35 AM

                            ad_right_rmr

                            Collapse

                            News

                            Collapse

                            Topics Statistics Last Post
                            Started by seqadmin, 10-11-2024, 06:55 AM
                            0 responses
                            11 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 10-02-2024, 04:51 AM
                            0 responses
                            110 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 10-01-2024, 07:10 AM
                            0 responses
                            114 views
                            0 likes
                            Last Post seqadmin  
                            Started by seqadmin, 09-30-2024, 08:33 AM
                            1 response
                            120 views
                            0 likes
                            Last Post EmiTom
                            by EmiTom
                             
                            Working...
                            X