Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • jmwhitha
    Senior Member
    • Mar 2013
    • 107

    Add Locus Column to Cummerbund diff_genes.txt

    Good morning SeqAnswers community,

    Today I am interested in adding a "locus" column to my diff_genes.txt file. I'm refering to the code at the top of pg. 575 in the 2012 Nature Protocol by Trapnell et al.

    >gene_diff_data<-diffData(genes(cuff_data))
    >sig_gene_data<-subset(gene_diff_data, (significant=='yes'))
    >write.table(sig_gene_data, 'diff_genes.txt', sep='\t', row.names=F, col.names=T, quote=F)

    The so called "locus" column is found in the genes_exp.diff file in the cuffdiff output. Again, I would like to capture this column in my diff_genes.txt. Furthermore, I would really like to have a column of "/product" ID's from the original GeneBank file, but that seems like it would involve a VBA script that would search the GeneBank file for the "locus tag" and then copy what's in quotations 5 rows down from it. I don't remember how to do that anymore though, so if anyone has any ideas before I manually ctrl find all my locus tags and product ID's, I'd really appreciate the help.

    I found this, but I don't know if it is relevant since there are no examples.

    4.1 Adding additional feature annotation
    Gene- or feature-level annotation can be permanently added to the database tables for future querying. If you have a data.frame where the first column contains the 'tracking_id' (eg. 'gene_id' for genes, 'isoform_id' for isoforms, etc). You can easily add feature level annotation using the addFeatures() function:

    > annot<-read.table("gene_annotation.tab",sep="\t",header=T,na.string="-")
    > addFeatures(cuff,annot,level="genes")


    By default, features added to a CuffSet object are assumed to be gene-level annotations, but the level can selected using the argument level. Features added to a CuffData object are assumed to be of the same type as the 'type' value for that given object (e.g. gene-level features for 'genes', isoform-level features for isoforms, etc.) - See more at: http://compbio.mit.edu/cummeRbund/ma....mojSSK12.dpuf

    Thanks and God bless,
    Jason

Latest Articles

Collapse

  • seqadmin
    Pathogen Surveillance with Advanced Genomic Tools
    by seqadmin




    The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
    Yesterday, 11:48 AM
  • seqadmin
    New Genomics Tools and Methods Shared at AGBT 2025
    by seqadmin


    This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

    The Headliner
    The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
    03-03-2025, 01:39 PM
  • seqadmin
    Investigating the Gut Microbiome Through Diet and Spatial Biology
    by seqadmin




    The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
    02-24-2025, 06:31 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-20-2025, 05:03 AM
0 responses
26 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-19-2025, 07:27 AM
0 responses
33 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-18-2025, 12:50 PM
0 responses
25 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-03-2025, 01:15 PM
0 responses
190 views
0 reactions
Last Post seqadmin  
Working...