Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Milestailsprowe
    Member
    • Jun 2015
    • 31

    #1

    Cuffdiff issue: Sort order of reads in BAMS must be the same

    Hello. I'm having issues with Cuffdiff and can not figure out why. The program says the read order has to be the same but I ran all of then with Samtools sort -n to make sure and still having that issue.

  • Brian Bushnell
    Super Moderator
    • Jan 2014
    • 2709

    #2
    You need to sort by coordinate, not by name (remove the -n flag). The header of the bam should look something like this afterward:

    @HD VN:1.0 SO:coordinate

    Comment

    • Milestailsprowe
      Member
      • Jun 2015
      • 31

      #3
      Originally posted by Brian Bushnell View Post
      You need to sort by coordinate, not by name (remove the -n flag). The header of the bam should look something like this afterward:

      @HD VN:1.0 SO:coordinate
      So do I sort the bam files using @HD VN:1.0 SO:coordinate with samtools?

      Comment

      • GenoMax
        Senior Member
        • Feb 2008
        • 7142

        #4
        No that is how the first header line should look after sorting.

        Just remove -n option when you sort with samtools like Brian said.

        Comment

        • blancha
          Senior Member
          • May 2013
          • 367

          #5
          You can check the headers of the BAM files with samtools view.
          Code:
          samtools view -H fileName.bam

          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...
            Yesterday, 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
          25 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-28-2026, 09:51 AM
          0 responses
          16 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-22-2026, 11:13 AM
          0 responses
          71 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...