Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • krish91
    Junior Member
    • Apr 2014
    • 3

    checking the differential expression from .bed files but using cufflinks/cuffdiff, ??

    Hi,
    I have been trying to see the differential expression of mrna seq by using there .bed files from epigenome data.

    I have converted .bed files to .bam using bedtools and i am using cufflinks for obtaining FPKM values and using cuffdiff for checking the differential expression, but both are resulting 0 FPKM values.
    command used:

    cuffdiff -o Homo_sapiens.GRCh37.75.gtf ips_cell_line/gsm706050_s_sorted.bam UW.Fetal_Spinal_Cord/gsm1101679_sorted.bam

    cufflinks -G Homo_sapiens.GRCh37.75.gtf UW.Fetal_Spinal_Cord/gsm1101679_sorted.bam

    version: cufflinks 1.3.0-2ubuntu1
    Last edited by krish91; 04-10-2014, 10:53 PM.
  • krish91
    Junior Member
    • Apr 2014
    • 3

    #2
    I have actually sorted my .sam file obtained from .bam by

    sort -k 3,3 -k 4,4n gsm1101679.sam > gsm1101679.sam.sorted
    even after sorting it i am getting the same error.
    ERROR:
    Error: this SAM file doesn't appear to be correctly sorted!
    current hit is at chr10:60106, last one was at chr1:249235968
    Cufflinks requires that if your file has SQ records in
    the SAM header that they appear in the same order as the chromosomes names
    in the alignments.
    If there are no SQ records in the header, or if the header is missing,
    the alignments must be sorted lexicographically by chromsome
    name and by position.

    can someone give me a clue how to proceed further..

    Comment

    Latest Articles

    Collapse

    • GATTACAT
      Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing
      by GATTACAT
      Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
      07-01-2026, 11:43 AM
    • SEQadmin2
      Nine Things a Sample Prep Scientist Thinks About Before Sequencing
      by SEQadmin2


      I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

      Here are nine questions we think about, in roughly the order they matter, before...
      06-18-2026, 07:11 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 07-02-2026, 11:08 AM
    0 responses
    25 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-30-2026, 05:37 AM
    0 responses
    23 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-26-2026, 11:10 AM
    0 responses
    23 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-17-2026, 06:09 AM
    0 responses
    55 views
    0 reactions
    Last Post SEQadmin2  
    Working...