Error with HTSeq

Running HTSeq to get raw read counts associated with chromosomal positions/annotated genes from sorted .bam files:

Error occurred when processing SAM input (record #30675291 in file {PATH}/tophat/Exp-1-Basal/accepted_hits.sorted.bam):
 Maximum alignment buffer size exceeded while pairing SAM alignments.
 [Exception type: ValueError, raised in]

Simon Anders, the author of the package from EMBL Heidelberg, suggested in this forum post to remove the “-r pos” or “–order=pos” argument. I did this and it seems to work fine now. Interesting that I had run HTSeq before with a different data set of sorted bam files and did not get this error.

HTSeq manual page:


About Lisa Johnson

PhD candidate at UC Davis in Molecular, Cellular, and Integrative Physiology
This entry was posted in Bioinformatics, software. Bookmark the permalink.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s