Question: Error: htseq-count exceeds memory buffer, Solution: Set sort option on tool form to Yes
0
gravatar for chaudhry
11 months ago by
chaudhry0
chaudhry0 wrote:

Hello, I am dealing with RNA-seq data and I was successfully able to align my raw FASTQ data to the built-in input hg19 genome using HISAT2. I am now having problems with measuring my reads using htseq-counts. I am using the hg19 GTF file for the genome and I keep getting either an error saying:

An error occurred with this dataset: Fatal error: Unknown error occured 100000 GFF lines processed. 200000 GFF lines processed. 300000 GFF lines processed. 400000 GFF lines processed. 500000 GFF lines processed. 600000 GFF lines processed. 700000 GFF lines processed. 800000 GFF lines processed.

Can anyone point me in the right direction? Thank you!

ADD COMMENTlink modified 11 months ago by Jennifer Hillman Jackson25k • written 11 months ago by chaudhry0
0
gravatar for Jennifer Hillman Jackson
11 months ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

The full error ends with:

Maximum alignment buffer size exceeded while pairing SAM alignments.
[Exception type: ValueError, raised in __init__.py:671]

Try setting the option Force sorting of SAM/BAM file by NAME to Yes.

Help from the tool form: This option can be used for for paired-end data that has many unmapped mates. Use this if you get the warning about paired end data missing or not being properly sorted.

Hope that helps! Thanks, Jen, Galaxy team

ADD COMMENTlink written 11 months ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 171 users visited in the last hour