Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • cebehind
    Junior Member
    • Oct 2015
    • 5

    #1

    Rare NCBI SRA dataset - how can I use it?

    Hi to all,

    I wanted to check this data and use cufflinks to get the FPKM values, but I haven't been able to. The SRA files are too light (200-300mb approximatly) and I know their usually around 1.5Gb.
    I have little to none experience with processing this kind of data and I am having a hard time coming up with a solution for that. I will really appreciate if someone can explain me how to do that.

    http://www.ncbi.nlm.nih.gov/geo/quer...i?acc=GSE68789 - This is the dataset.

    Thanks a lot.
    Last edited by cebehind; 10-14-2015, 06:09 AM.
  • cebehind
    Junior Member
    • Oct 2015
    • 5

    #2
    Originally posted by mastal
    It looks like the files you are referring to are files with the counts, not the raw data.


    That's the link to SRA data, wich is supposed to be the raw files. But I'm not sure.
    Last edited by cebehind; 10-14-2015, 06:53 AM. Reason: Quote missing.

    Comment

    • mastal
      Senior Member
      • Mar 2009
      • 666

      #3
      The first file listed looked like it was a summary of the counts.

      The other files are .sra files, and I think they might be in a compressed format.

      I think you need the sra toolkit to convert the files to something like fastq.

      Have a look at:


      Comment

      • WhatsOEver
        Senior Member
        • Apr 2012
        • 215

        #4
        Originally posted by cebehind View Post
        http://www.ncbi.nlm.nih.gov/sra?link...rom_uid=283850

        That's the link to SRA data, wich is supposed to be the raw files. But I'm not sure.
        Yes, that is the link to the datafile in the SRA. There you can download the files and convert them to fastq via fastq-dump from the sra toolkit. After that you can go for mapping->cufflinks->whatsoever...

        Alternatively, you can follow your first link to the geo deposit and download the counts table.

        The read numbers are indeed not that high, but for DE it was probably enough - at least it was published with this data amount

        Comment

        • cebehind
          Junior Member
          • Oct 2015
          • 5

          #5
          Originally posted by WhatsOEver View Post
          Yes, that is the link to the datafile in the SRA. There you can download the files and convert them to fastq via fastq-dump from the sra toolkit. After that you can go for mapping->cufflinks->whatsoever...

          Alternatively, you can follow your first link to the geo deposit and download the counts table.

          The read numbers are indeed not that high, but for DE it was probably enough - at least it was published with this data amount
          http://www.ebi.ac.uk/ I used the European Bioinformatics Institute to download the fastq files. Would that work?

          Comment

          • mastal
            Senior Member
            • Mar 2009
            • 666

            #6
            yes, that should be fine.

            Comment

            • cebehind
              Junior Member
              • Oct 2015
              • 5

              #7
              I'm planing on using Galaxy then, I have no idea how to use R. Is there something that I must know beforehand?

              Comment

              • WhatsOEver
                Senior Member
                • Apr 2012
                • 215

                #8
                The EBI hosts those as fastq? That's surprising but should be fine if you are sure the files are the correct ones. From the geo site you get some info on their analysis procedure which will be important if you want to (more or less) exactly reproduce the results (tophat --max-intron-length 10000 --max-multihits 1; Expression table done using ESAT http://garberlab.umassmed.edu/software/esat/; Genome_build: mm9).

                Comment

                • cebehind
                  Junior Member
                  • Oct 2015
                  • 5

                  #9
                  Originally posted by WhatsOEver View Post
                  The EBI hosts those as fastq? That's surprising but should be fine if you are sure the files are the correct ones. From the geo site you get some info on their analysis procedure which will be important if you want to (more or less) exactly reproduce the results (tophat --max-intron-length 10000 --max-multihits 1; Expression table done using ESAT http://garberlab.umassmed.edu/software/esat/; Genome_build: mm9).
                  The thing is I wasn't able to run the ESAT I don't know why. So I'm aiming to get FPKM values for that. I'm having problems with the gene names in cufflinks. I don't know how to do that.

                  Comment

                  Latest Articles

                  Collapse

                  • SEQadmin2
                    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
                    by SEQadmin2




                    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
                    Today, 06:46 AM
                  • SEQadmin2
                    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
                    by SEQadmin2



                    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

                    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
                    ...
                    10-14-2026, 09:55 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by SEQadmin2, 10-28-2026, 09:57 AM
                  0 responses
                  22 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-28-2026, 09:51 AM
                  0 responses
                  15 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-22-2026, 11:13 AM
                  0 responses
                  68 views
                  0 reactions
                  Last Post SEQadmin2  
                  Started by SEQadmin2, 10-22-2026, 10:48 AM
                  0 responses
                  48 views
                  0 reactions
                  Last Post SEQadmin2  
                  Working...