FeatureCounts Assertion `l_name < 100' failed.
0
0
Entering edit mode
shg018 • 0
@shg018-11542
Last seen 7.6 years ago

Hi I have been trying to run featureCounts on my mapped bam extension files with the command: featureCounts -s 2 -t exon -g gene_id -a /oasis/.../genes_l2.gtf -o /oasis/.../CompleteRNA.txt file.bam

and have been getting:

    ==========     _____ _    _ ____  _____  ______          _____  
    =====         / ____| |  | |  _ \|  __ \|  ____|   /\   |  __ \ 
      =====      | (___ | |  | | |_) | |__) | |__     /  \  | |  | |
        ====      \___ \| |  | |  _ <|  _  /|  __|   / /\ \ | |  | |
          ====    ____) | |__| | |_) | | \ \| |____ / ____ \| |__| |
    ==========   |_____/ \____/|____/|_|  \_\______/_/    \_\_____/
  v1.5.1

//========================== featureCounts setting ===========================\ || || || Input files : 1 BAM file || || P FF1_1_R.bam || || || || Output file : /oasis/ ... || || Summary : / ... || || Annotation : /oasis/ ... || || Dir for temp files : /oasis/ ... || || || || Threads : 1 || || Level : meta-feature level || || Paired-end : no || || Strand specific : reversely stranded || || Multimapping reads : not counted || || Multi-overlapping reads : not counted || || Overlapping bases : 0.0% || || || \===================== http://subread.sourceforge.net/ ======================//

//================================= Running ==================================\ || || || Load annotation file /oasis/.../genes_l2. ... || || Features : 690739 || || Meta-features : 46170 || || Chromosomes/contigs : 165 || || || || Process BAM file FF1_1_R.bam... || || Paired-end reads are included. || || Assign reads to features... || featureCounts: readSummary.c:1578: process_pairer_header: Assertion `l_name < 100' failed. /var/spool/torque/mom_priv/jobs/6660772.tscc-mgr.local.SC: line 18: 668 Aborted featureCounts -s 2 -t exon -g gene_id -a /oasis/.../genes_l2.gtf -o /oasis/.../CompleteRun1.txt ${file}


I was wondering about the possible reason for said error. Can it be due to low mapping percentage in bam file.

featurecounts rnaseq • 1.4k views
ADD COMMENT
0
Entering edit mode

Could you please provide your bam file and gtf anotation file offline so that we can take a close look?

ADD REPLY

Login before adding your answer.

Traffic: 891 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6