Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • NOTEST on Cuffdiff output

    I am new to the Tuxedo Suite and decided to start by going through the workflow in the Nature Protocols publication. All seemed to be going well until I got to the point of differential expression analysis with Cuffdiff. There, the output file gene_exp.diff had NOTEST in the status column for almost all genes. I looked at the Cufflinks manual, which indicated that this output can result from not having enough alignments for testing in any given gene. I then used Samtools to be sure the .bam files that were generated by TopHat had a number of mapped reads that was comparable to the Nature Protocols example. The results, for all accepted_hits.bam files suggested that the number of mapped reads was very similar to what was expected according to the protocol. An example is below:

    $ samtools idxstats C2_R3_thout/accepted_hits.bam

    2L 23011544 4591000 0
    2R 21146708 5029664 0
    3L 24543557 4112826 0
    3R 27905053 5317642 0
    4 1351857 194612 0
    M 19517 0 0
    X 22422827 4157222 0

    I then tried to use the Cuffdiff option -c, to reduce the minimum number of alignments needed at a locus from 10 to 5, but the resulting gene_exp.diff file still had NOTEST in the vast majority of genes.

    In case it helps, the Cuffdiff command that I used was as follows:

    $ cuffdiff -o diff_out -b ./Drosophila_melanogaster/Ensembl/BDGP5.25/Sequence/Bowtie2Index/genome.fa -p 8 -L C1,C2 -u merged_asm/merged.gtf ./C1_R1_thout/accepted_hits.bam,./C1_R2_thout/accepted_hits.bam,./C1_R3_thout/accepted_hits.bam ./C2_R1_thout/accepted_hits.bam,./C2_R2_thout/accepted_hits.bam,./C2_R3_thout/accepted_hits.bam

    If any additional information is necessary for you to help me troubleshoot this issue, please let me know.

  • #2
    I got the similar issue here. Any idea?

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Recent Advances in Sequencing Analysis Tools
      by seqadmin


      The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
      05-06-2024, 07:48 AM
    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin




      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
      04-22-2024, 07:01 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 05-14-2024, 07:03 AM
    0 responses
    24 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-10-2024, 06:35 AM
    0 responses
    44 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-09-2024, 02:46 PM
    0 responses
    58 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 05-07-2024, 06:57 AM
    0 responses
    44 views
    0 likes
    Last Post seqadmin  
    Working...
    X