Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bam2fastq discarded reads

    Hi all,

    I've been using bam2fastq on my tophat output and it's been great, runs really quickly, except for the number of reads being discarded. For example this was for one of my output files from tophat

    This looks like paired data from lane 239.
    Output will be in x_1.fastq and x_2.fastq
    60465861 sequences in the BAM file
    60465861 sequences exported
    WARNING: 6585459 reads could not be matched to a mate and were not exported

    That's 10% of the reads being discarded, and in other files it's even more (I ran it on another file just now and 17% of the reads were discarded). What I don't understand is that the PE files which were put into tophat were quality filtered with software to directly handle PEs and so both files have the same number of reads and all the reads have mates and both PE files are the same order (tophat freaks out otherwise) so why is bam2fastq discarding these reads? If any reads didn't have a mate then tophat would have returned an error.
    Last edited by Derek-C; 01-14-2013, 04:47 AM.

  • #2
    I remember Bam2fastq only export read pair. So, all singleton mapped will be discarded. In your case, it may be better to write a small script to parse the singleton out
    Marco

    Comment


    • #3
      Originally posted by marcowanger View Post
      I remember Bam2fastq only export read pair. So, all singleton mapped will be discarded. In your case, it may be better to write a small script to parse the singleton out
      Are you certain about the singletons thing? It would make sense, butl I just did some checking there with samtools flagstat and according to the numbers I'm looking at for this one accepted_hits.bam file from tophat there were more singleton reads then there were reads discarded by bam2fastq (~3 million more). So if bam2fastq doesn't export singletons, what happened to those 3 million reads.

      Actually there are options in tophat for setting it so that there should be no singletons, hopefully that'll fix this.
      Last edited by Derek-C; 01-14-2013, 07:29 AM.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Exploring the Dynamics of the Tumor Microenvironment
        by seqadmin




        The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
        07-08-2024, 03:19 PM
      • seqadmin
        Exploring Human Diversity Through Large-Scale Omics
        by seqadmin


        In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
        06-25-2024, 06:43 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 07-10-2024, 07:30 AM
      0 responses
      19 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-03-2024, 09:45 AM
      0 responses
      197 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-03-2024, 08:54 AM
      0 responses
      207 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 07-02-2024, 03:00 PM
      0 responses
      191 views
      0 likes
      Last Post seqadmin  
      Working...
      X