Tweet
subread-1.5.0-p1-source works
Thank you Wei Shi, the new patched version (subread-1.5.0-p1-source) works for these files now.
Thank you Wei Shi, the new patched version (subread-1.5.0-p1-source) works for these files now.
-
-
That's great. Thanks for letting us know.Comment
-
Error with long CIGARs
I'm trying to map some long reads (PacBio public human data), and use featureCounts to count reads to features (there are other ways of doing this, which I also use, but I'd like to use featureCounts on both the PacBio data and some IBM data to treat both datasets consistently. However, I'm getting the following error:
Can you please update featureCounts to not collapse when a CIGAR is too long?Code:WARNING: cigar string is too long to the buffer. Please only use the compressed BAM format. featureCounts: sambam-file.c:594: PBam_chunk_gets: Assertion `0' failed.
Thanks in advance,
DaryaComment
-
What is the version of featureCounts you use?Comment
-
Segmentation fault (core dumped)
If I run the below command from the commandline it works fine, however it corrupts with "Segmentation fault (core dumped)" if I run within sh script:
Code:featureCounts -T 10 -Q 10 -a library/genes.gtf -o gene_expression.txt 344/RNA-Seq/CD34_2_clean/CD34_2_clean_aligned_sorted.bam
Apparently the patch does not work on my system.Code:========== _____ _ _ ____ _____ ______ _____ ===== / ____| | | | _ \| __ \| ____| /\ | __ \ ===== | (___ | | | | |_) | |__) | |__ / \ | | | | ==== \___ \| | | | _ <| _ /| __| / /\ \ | | | | ==== ____) | |__| | |_) | | \ \| |____ / ____ \| |__| | ========== |_____/ \____/|____/|_| \_\______/_/ \_\_____/ v1.5.0-p3 //========================== featureCounts setting ===========================\\ || || || Input files : 1 BAM file || || P 344/RNA ... || || || || Output file : 344/RNA-S ... || || Summary : 344/RNA-S ... || || Annotation : library/genes.gtf (GTF) || || || || Threads : 10 || || Level : meta-feature level || || Paired-end : yes || || Strand specific : no || || Multimapping reads : not counted || || Multi-overlapping reads : not counted || || || || Chimeric reads : counted || || Both ends mapped : not required || || || \\===================== http://subread.sourceforge.net/ ======================// //================================= Running ==================================\\ || || || Load annotation file library/genes.gtf ... || || Features : 430178 || || Meta-features : 23710 || || Chromosomes/contigs : 49 || || || || Process BAM file 344/RNA-Seq/CD34_2 ... || || Paired-end reads are included. || || Assign fragments (read pairs) to features... || || || || WARNING: reads from the same pair were found not adjacent to each || || other in the input (due to read sorting by location or || || reporting of multi-mapping read pairs). || || || || Read re-ordering is performed. || || || Segmentation fault (core dumped)
Is there any clever fix yet?Last edited by arkanion; 08-04-2016, 07:18 AM.Comment
-
Rsubread_1.20.6 or command line version subread-1.5.0-p1-Linux-x86_64Originally posted by shi View PostComment
-
Can you update to Subread 1.5.1 and try again?Comment
-
I'm also getting a segmentation error.
I use featureCount on 33k files and only 3 of them return this error. I'm getting no stdout or err on these 3 files, just:
44818 Segmentation fault (core dumped)
The query I use:
featureCounts -s 1 -T 8 -t exon -g gene_id -a Homo_sapiens.GRCh37.75.gtf bamfile.bam
Version: 1.5.1
Platform: Linux-x86_64
Link to one of the problem bam files:
Any help would be highly appreciatedLast edited by hbrugge; 11-15-2016, 01:48 AM.Comment
-
It would be really nice if anyone could look into this problemComment
-
I also encounter similar problem:
But after I sort the bam file using PICARD, featureCounts worked fine.Code:SAM_pairer_run: Assertion `1 != corrected_run'
Comment
-
Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.
The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...06-02-2026, 10:05 AM -
With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.
Introduction
Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...05-22-2026, 06:42 AM -
Studying ecosystems means dealing with complex, multi-species communities that are hard to observe at scale. This complexity, however, hides many important questions to be answered, from how biogeochemical cycles work and how climate change can affect species distribution to how conservation strategies can work best.
Genomics, particularly since the expansion of NGS, has transformed ecosystem ecology. By sequencing environmental DNA, we can now assess biodiversity without direct...05-06-2026, 09:04 AM
Comment