I am new of RNA-seq analysis. I am running HTSeq-count on galaxy with defaut parameters. I use GenePattern Tophat to align my data. The accepted_hit.bam is about 2.55 GiB. And it turns out a fatal error as following:
5700000 SAM alignment record pairs processed. Error occured when processing SAM input (record #14420795 in file /galaxy-repl/main/files/019/167/dataset_19167797.dat): Maximum alignment buffer size exceeded while pairing SAM alignments. [Exception type: ValueError, raised in __init__.py:671]
I don't know what's happened. How to solve this problem? Your advice are appreciated.