SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
htseq-count paolo.kunder Bioinformatics 10 10-22-2014 04:45 AM
HTSeq: mismatched quotes issue? Artur Jaroszewicz RNA Sequencing 11 09-26-2014 11:23 AM
Which ID should be used for HTSeq-count? syintel87 Bioinformatics 11 02-07-2013 12:16 AM
Issue with htseq-count cpleis Bioinformatics 8 10-15-2012 09:31 AM
multiBamCov or htseq-count to count read per feature ? NicoBxl Bioinformatics 1 07-03-2012 02:05 AM

Reply
 
Thread Tools
Old 08-15-2013, 06:12 AM   #1
zuco
Junior Member
 
Location: Switzerland

Join Date: Oct 2011
Posts: 2
Default HTSeq-count issue

First of all,

hello to everyone, I am beginning to do RNAseq analysis in bacteria and this forum has been of immense help so far, unfortunately I have not seen anyone with this issue yet:

My system: Windows 7, running BioLinux 7 in Virtual Machine.

Python 2.7.3
HTSeq 0.5.4p3
Numpy 1.6.1

I was trying to count my mapped reads using a pipe with samtools:

samtools view -h sorted.bam | htseq-count -m intersection-nonempty - reference.gff

The error I got was:

4410 GFF lines processed.
Error: 'itertools.chain' object has no attribute 'get_line_number_string'
[Exception type: AttributeError, raised in count.py:199]




Not sure where to go from here, it feels like I may be missing a module or something.

Any and all help is appreciated!
zuco is offline   Reply With Quote
Old 08-15-2013, 07:02 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

I've had that happen on truncated BAM files before. Have a look at the GFF file after line 4410 (or maybe at the beginning of the BAM file, if that's near the end of your GFF) to see if something is amiss.
dpryan is offline   Reply With Quote
Old 08-15-2013, 11:30 PM   #3
zuco
Junior Member
 
Location: Switzerland

Join Date: Oct 2011
Posts: 2
Default

Hi,

Quote:
I've had that happen on truncated BAM files before.
you were spot on. Turns out the unmapped reads were the issue, I fixed it by changing my pipe so that the samtools command included -F 4:

samtools view -h -F 4 sorted.bam | htseq-count -m intersection-nonempty - reference.gff

Thank you!
zuco is offline   Reply With Quote
Reply

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:59 AM.


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