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 06:52 PM
combining htseq-count files into one mbk0asis General 6 06-17-2013 08:07 AM
Using HTSeq with Bowtie2 .sam files all_your_base Bioinformatics 8 10-31-2012 01:39 PM
Error using HTSeq on Mapsplice alignments AnneBiton Bioinformatics 5 08-24-2012 12:34 AM
NEw to Chip-seq and have .bam/.sam/.bam.bai files... then what? NGS newbie Bioinformatics 11 05-25-2011 08:48 AM

Reply
 
Thread Tools
Old 02-07-2014, 07: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, 09:21 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,800
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, 09: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 09:36 AM.
narges is offline   Reply With Quote
Old 02-07-2014, 10:24 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,800
Default

I am checking with TCGA group for a solution.
GenoMax is offline   Reply With Quote
Old 02-07-2014, 12:16 PM   #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, 09:38 AM   #6
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,800
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 06:01 PM.


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