Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Parashar
    Junior Member
    • Feb 2014
    • 5

    Sorted bamfiles have increased size!!

    BACKGROUND:
    I have RNA-Seq data from Illumina platform.
    I aligned the reads post-QC onto hg19 bowtie index using tophat.
    Now I wished to get the raw count of reads mapped to each gene using Ht-seq-count.

    WHAT I DID:
    However the HT-Seq reports error that it is unable to find the mate pair and asks whether the sam file is properly sorted.
    Hence, after reading posts on the problem. I sorted my bam files using samtoools.
    samtools sort -n accepted_hits.bam accepted_hits_bam_sorted

    CONCERN:
    Now the sorted BAM file has increased in size and I wish to know why his has happened.
    Though it might not be but seems highly unintuitive.
  • dpryan
    Devon Ryan
    • Jul 2011
    • 3478

    #2
    That's not surprising. Remember that BAM files are compressed, with the compression algorithm being more efficient the more similar neighbouring reads are (i.e., coordinate-sorted files should compress to a greater extent).

    Comment

    • Parashar
      Junior Member
      • Feb 2014
      • 5

      #3
      Thanks dpryan.
      Just wanted to confirm if something went horribly wrong.
      I'll now be following your earlier reply on a different post to run HTSeq

      Comment

      • maubp
        Peter (Biopython etc)
        • Jul 2009
        • 1544

        #4
        You sorted by read name (-n), but the default sorting by mapping position puts similar reads together, and therefore their sequence data compresses well, so coordinate sorted BAM files are usually smaller.

        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.
          Today, 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
        • SEQadmin2
          From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
          by SEQadmin2


          Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


          The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
          ...
          06-02-2026, 10:05 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Yesterday, 05:37 AM
        0 responses
        7 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 06-26-2026, 11:10 AM
        0 responses
        17 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 06-17-2026, 06:09 AM
        0 responses
        51 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 06-09-2026, 11:58 AM
        0 responses
        110 views
        0 reactions
        Last Post SEQadmin2  
        Working...