Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • gffread to output sequence, gene_id not output

    HI everyone,

    Many apologies if I'm duplicating, I have searched the forums, google, can't find the specific answer.

    So- I've performed my mRNAseq experiment, used the workflow:

    cufflinks->cuffmerge->cuffquant->cuffdiff

    then used cummeRbund to look at the results.

    From cummeRbund I've generated a list of differentially expressed genes.

    What I'd like to do now is look at the sequence of the genes to see what type of things are differentially expressed (have done a brief GO analysis, would like to search HMM profiles for protein motifs).

    I tried to output the sequences from my merged.gtf file (generated by cuffmerge) using gffread. I can get them to output, but I would really, REALLY, like the gene_id "XLOC_*****" number to be in the fasta header. But it seems that whatever I do, I can't get it out there. I can get almost every single other piece of info from the gtf file there using one or other of the gffread options, but not this.

    Clearly it wouldn't be so hard to write my own script to do this, but I'm under time pressure, and I've leaernt the hard way that duplicating others efficient tools is foolhardy.

    So- am I missing the crucial option here? Or do folks do this (outputting differentially exporessed gene sequences from mRNAseq expts) iin a different way?

    I do have the gene IDs of the annotated genes in the fasta header, but there are some novel/intergenic/anomalous genes which are only really iddentifiable by "XLOC****"

    Many thanks for your help

    Matt

  • #2
    Hi, running into to a similar problem.

    There is a workaround here that may be relevant:



    I can't get my merged.gtf file (output of cuffmerge) to output anything using gffread. I just get empty files.

    I tried to output the sequences from my merged.gtf file (generated by cuffmerge) using gffread. I can get them to output
    Can you detail the code you used to do this?

    Comment


    • #3
      gffread merge_matt_tissue_MSU/merged.gtf -g Oryza_sativa.IRGSP-1.0.21.dna_sm.genome.fa -w test.fa

      for example.

      The fasta genome file is in the same directory- you have to supply the full path if not.

      Depending on which options you use, it seems a little sensitive to the features in your gtf file. For example, using -x gave me noting here as the merged.gtf doesn't have any features labelled as cds, only exons.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Genetic Variation in Immunogenetics and Antibody Diversity
        by seqadmin



        The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
        11-06-2024, 07:24 PM
      • seqadmin
        Choosing Between NGS and qPCR
        by seqadmin



        Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
        10-18-2024, 07:11 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 11-01-2024, 06:09 AM
      0 responses
      30 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 10-30-2024, 05:31 AM
      0 responses
      21 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 10-24-2024, 06:58 AM
      0 responses
      26 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 10-23-2024, 08:43 AM
      0 responses
      57 views
      0 likes
      Last Post seqadmin  
      Working...
      X