Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • tophat2 bowtie option help

    Hi,

    I am using Tophat2 on the grid and trying to get single and paired read unmapped bam files by a following command;

    tophat2 -p 16 -o 4560_output -G ~/GRCh37_73.gtf ~/Bowtie2Index2/Ensemble 4650_R1merged.fastq 4650_R2merged.fastq --b2-un-conc ~/4560_output/pefailed.bam --b2-un ~/4560_output/singlefailed.bam

    But I only get unmapped.bam but not pefailed.bam and singlefailed.bam. I really appreciate any help with this.

    Thanks in advance.

  • #2
    "--b2-un-conc" and "--b2-un" aren't options for tophat2. They look like bowtie2 options with "b2-" slapped on the front. That simply won't work as I don't think there's a way to get tophat to do that without hacking the source code a fair bit. However, it would be relatively straight-forward to simply extract that from the unmapped.bam file, which should be name-sorted.

    Comment


    • #3
      Hi Devon,

      I am new to Samtools. Could you send me a pointer of command that do these ?

      Thanks.

      Aki

      Comment


      • #4
        Hi Aki,

        I don't have any paired-end data aligned with tophat2 at my finger-tips, but I imagine if you "samtools sort -n unmapped.bam" to ensure that the file is name-sorted, you would then be able to write a small script (in python/perl/whatever) to then print out reads that either do or do not have a mate (i.e., read with the same name).

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Recent Advances in Sequencing Analysis Tools
          by seqadmin


          The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
          05-06-2024, 07:48 AM
        • 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...
          04-22-2024, 07:01 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 05-14-2024, 07:03 AM
        0 responses
        26 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-10-2024, 06:35 AM
        0 responses
        45 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-09-2024, 02:46 PM
        0 responses
        59 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-07-2024, 06:57 AM
        0 responses
        46 views
        0 likes
        Last Post seqadmin  
        Working...
        X