SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
HTSeq-count on bam dvanic Bioinformatics 18 11-12-2014 05:52 PM
combining htseq-count files into one mbk0asis General 6 06-17-2013 07:07 AM
Using HTSeq with Bowtie2 .sam files all_your_base Bioinformatics 8 10-31-2012 12:39 PM
Error using HTSeq on Mapsplice alignments AnneBiton Bioinformatics 5 08-23-2012 11:34 PM
NEw to Chip-seq and have .bam/.sam/.bam.bai files... then what? NGS newbie Bioinformatics 11 05-25-2011 07:48 AM

Reply
 
Thread Tools
Old 02-07-2014, 06:57 AM   #1
narges
Member
 
Location: Finland

Join Date: Aug 2012
Posts: 29
Default Does HTSeq accepts BAM files from Mapsplice?

Hi,

I have BAM files from Mapsplice and I wanted to quantify the expression levels with HTSeq. But I received following error. Do you think it is regarding the aligner or what is the problem?

samtools view file.bam | htseq-count --stranded=no hg19/Annotation/Genes/genes.gtf > test.txt

Warning: Read file:8:2203:10838:61544/2 claims to have an aligned mate which could not be found. (Is the SAM file properly sorted)
Error: 'itertools.chain' object has no attribute 'get_line_number_string'
[Exception type: AttributeError, raised in count.py:200]

Thanks in advance
narges is offline   Reply With Quote
Old 02-07-2014, 08:21 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,784
Default

I assume these are RNA-seq BAM's. See the note about Mapsplice SAM incompatibility in this note: https://cghub.ucsc.edu/docs/tcga/UNC...eq_summary.pdf
GenoMax is offline   Reply With Quote
Old 02-07-2014, 08:34 AM   #3
narges
Member
 
Location: Finland

Join Date: Aug 2012
Posts: 29
Default

Quote:
Originally Posted by GenoMax View Post
I assume these are RNA-seq BAM's. See the note about Mapsplice SAM incompatibility in this note: https://cghub.ucsc.edu/docs/tcga/UNC...eq_summary.pdf
Yes, they are and I have seen the document. I assumed it is not possible to "reconstruct the fastq files using samtools" but not other purposes like "samtools view ...".
But is there an alternative pipeline rather than samtools for quantifying the expression levels of BAM files from Mapsplice?

Last edited by narges; 02-07-2014 at 08:36 AM.
narges is offline   Reply With Quote
Old 02-07-2014, 09:24 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,784
Default

I am checking with TCGA group for a solution.
GenoMax is offline   Reply With Quote
Old 02-07-2014, 11:16 AM   #5
narges
Member
 
Location: Finland

Join Date: Aug 2012
Posts: 29
Default

Quote:
Originally Posted by GenoMax View Post
I am checking with TCGA group for a solution.
Many thanks for the help.
narges is offline   Reply With Quote
Old 02-10-2014, 08:38 AM   #6
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,784
Default

I have been told that you can use bedtools (coverageBed) with a BED format file. That should work with MapSplice generated BAM's.
GenoMax is offline   Reply With Quote
Reply

Tags
htseq-count, rna-seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:52 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO