As part of their submission requirements ENA/ArrayExpress require BAMs to include all reads (e.g. also unmapped ones). TopHat in its wisdom splits unmapped reads from the accepted_hits.bam output.
However, creating a 'complete' BAM is not so straightforward due to problems in the unmapped reads BAM (detailed on seqanswers). I've tried the script mention on SeqAnswers, but it doesn't solve all the problems for me. I still get unmapped reads which don't have a mate in accepted_hits.bam
Given TopHat's popularity this should be a solved problem, right? So, how does everyone deposit their TopHat alignments for publication? Or are fastqs the preferred format.
Any suggestions?
•
written
by
Chris Cole
1 answer
With the help of Christian Brueffer, I managed to resolve this problem. See his github: https://github.com/cbrueffer/tophat-recondition/issues/1#issuecomment-129849801
•
written
by
Chris Cole
Log in to answer this question.
-
Forum: Troubleshooting Tips - bcl2fastq creates duplicate reads
written by DavidStreid •
-
Request to validate my RNAseq workflow
written by umn_bist •
-
HTSeq-count, problem in handling PE data
written by barbara.mariotti81 •
-
Spliced aligners for RNA-seq
written by amirmhzadeh •
-
Extract reads with variant at given position
written by philip.jonsson •
-
QC based on Samtools Flagstat
written by indranipoddar •
-
Tophat : Sam-Flag 115 = Properly-Paired + Read.Reverse + Mate.Reverse ?
written by Pierre Lindenbaum
-
Samtools Count Paired-End Reads
written by dfernan
-
Confused About Solid Rna-Seq Pair End Data
written by Conan •
-
Bam Index Files (Bai) And Picard Tools Buildbamindex
written by Arun
