Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • assembly unmapped reads from tophat

    Hi,
    I have used Tophat/Cufflink to map and assembly my RNA-seq pair-end data(the insert size is 200bp while each read has 90bp length) with a reference genome. Now I am thinking how to de novo assembly the unmapped reads so that i can find some noval transcripts out of the reference genome. My tophat version is 2.0.8, so i got a output file named unmapped.bam. I am running bam2fastx with -o and -P options but running into errors. It generates a fastq without -P option, but i am not sure whether i can use this single fastq file to do a de novo assembly. Could anyone help me?
    Last edited by leifive; 03-07-2013, 06:17 AM. Reason: edit the title

  • #2
    I used bam2fastq tools, it can detect whether the reads in the bam file are pair end, and generate reads1 and reads2 file automatically if it is true.

    Comment


    • #3
      Originally posted by yzzhang View Post
      I used bam2fastq tools, it can detect whether the reads in the bam file are pair end, and generate reads1 and reads2 file automatically if it is true.
      Thanks, yzzhang. I have installed bam2fastq and running on my 500MB size unmapped.bam file. But converting it to paired end fastq files is taking much longer than the original tophat run which was only last 2 hours. One night has passed, and now the pair-end file size only reach 15MB. Is there someting wrong? or all i need is just waiting?

      Comment


      • #4
        Oh, It's my fault. I forget to run bam2fastq with -readname option. It's OK now. Thanks again.

        Comment

        Latest Articles

        Collapse

        • 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 on Modified Bases...
          Yesterday, 07:01 AM
        • seqadmin
          Current Approaches to Protein Sequencing
          by seqadmin


          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
          04-04-2024, 04:25 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        39 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        41 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        35 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-04-2024, 09:00 AM
        0 responses
        55 views
        0 likes
        Last Post seqadmin  
        Working...
        X