Paired end alignment issue

Hi,

I have paired end data which im trying to align using tophat:

tophat -r 200 -p 4 -o /.../Tophat_output/ -G /pathtogtf /pathtobowtieindexfiles Sample10_untrim1.fastq Sample10_untrim2.fastq

Im getting only one bam alignment file. Im not sure if I should be getting 2 bam files.

paired end sequencing bowtie tophat alignment

• 2,394 views • link

updated 10.2 years ago by Antonio R. Franco

• written 10.2 years ago by aj123

1 answer

You should get the accepted bam file, and another one of unmapped sequences

But don't expect two bam accepted files, one for each of the paired sequences. You only get one

• 0 views • link

written 10.2 years ago by Antonio R. Franco

Does it matter if the unmapped sequences files are larger than the accepted hits.bam file?

• 0 views • link

written 10.2 years ago by aj123

Possibly. Look at the alignment_summary.txt file that should have been produced. That will tell you what happened during the alignments.

• 0 views • link

written 10.2 years ago by GenoMax

yes thank you that answered my question.

• 0 views • link

written 10.2 years ago by aj123

More posts like this

hisat2 error-Paired end reads not equal or lack of RAM/disk space?
written 10 months ago by aj123

41467051 SRR8571937_1.fastq.gz 41030872 SRR8571937_2.fastq.gz My paired end reads are not equal. Please let me know if there is something I can do about this because …
Demultiplex based on sequence, not headers
written 5.9 years ago by damt0320 •

Hi. Im trying to demultiplex a dataset but im having a couple of problems with my files. I have reallized that each pair of barcodes …
Removing fused transcripts from RNA seq data
written 9.3 years ago by Max Ivon

Hello! I have performed alignment of paired-end data with tophat with the following options: tophat --no-coverage-search --max-multihits 1 -p 4 -o <...> <...> <...> -G …
low percentage mapping... tophat2.1.1
written 9.7 years ago by mra8187 •

hi im working on RNA-seq data of honey bee apis.mellifera.mellifera for trim the data i am using trimmomatic and it is the resaults : single …
RNA Seq pipeline
written 10.2 years ago by aj123

Hi, Im trying to create a simple RNA Seq pipeline from fastqc to DE in python (or perl). I want to read in a configuration …
length of quality sequence and length of read do not match in fastq
written 10.4 years ago by ifudontmind_plzz

I am trying to use cutadapt to filter low quality read from the paired end fastq files (paired end human exome sequence) . But im …
paired-end alignment with tophat2
written 10.5 years ago by zizigolu

Hi, I have 12 paired-end fq files. In tophat2 protocol this is the syntax tophat -p 8 -G genes.gtf genome file-1.fq file-2.fq but I read …
Both Paired And Unpaired Read Data Input For Subread Rna Sequencing Read Alignment
written 12.6 years ago by alpha2zee •

Does anyone know if the Subread RNA seq. aligner software accepts both paired and unpaired read data-files for the same sample, the way Tophat does? …
Aligning Paired-End Reads In Single-End Mode
written 13.5 years ago by GPR

Hello, I have a question on how to align paired-end reads. In cases of very large fastq files, which make aligners like TopHat crash in …
How To Generate A Paired End File Suitable For Bwa
written 13.8 years ago by GPR

A question about BWA and paired-end reads. In TopHat, one inputs the paired-end fastq file twice and the tool understands you have paired-end reads. I …