FeatureCounts error
Hi,
I'm trying to run FeatureCounts with my nanopore data (cDNA sequencing). I run the command:
featureCounts -g gene_id -Q 3 -a Homo_sapiens.GRCh38.97.gtf -o feature_counts_CLL-0849-02 CLL-0849-02.aligned.bam
and an error occurs:
Assertion failed: (NULL == pairer -> bam_margin_table -> appendix2), function SAM_pairer_finish_margin_table, file input-files.c, line 4612.
What does it mean? What should I do?
• 2,484 views
•
link
updated
by
ATpoint
•
written
by
m.wlodarczyk214 •
•
written
by
m.wlodarczyk214 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNASeq featureCounts Help Needed
written by SomeOne •
-
RNA-seq using rasflow - issues with annotation file
written by cdeantoneo31 •
-
Why does featurecounts give me an output file with only 0s?
written by Shraddha •
-
featureCounts error with 150bp Paired end RNA seq data
written by Ankur B. Sharma •
-
Question: featureCounts error: Assertion `1 != corrected_run' failed.
written by giuseppe0525 •
-
Low fragment assignment RNAseq
written by said3427
-
Large number of Unassigned_NoFeature reads from featureCounts
written by price0416 •
-
FeatureCounts Assertion `l_name < 100' failed.
written by shg018 •
Sorry for my initial comment, I did not notice the input bam file was indeed present. Can you do two things: a) make sure you have the most recent version of the
featureCountsand 2) be sure the BAM file is not corrupted, e.g. by runningsamtools quickcheck -qvvv your.bamand then posting the output?I have samtools 1.11. What command should I run to check my .bam files? I have tried with samtools qiuckcheck -v CLL-0849-02.aligned.bam but there is no output...
Sorry, edited it, it is
quickcheck, I had no coffee today, sorry for my poor-quality comments :)I appreciate your help even without coffee :) So I run the command:
samtools quickcheck -qvvv CLL-0849-02.aligned.bam
And the output looks like that:
verbosity set to 3 checking CLL-0849-02.aligned.bam opened CLL-0849-02.aligned.bam CLL-0849-02.aligned.bam is sequence data CLL-0849-02.aligned.bam has 194 targets in header. CLL-0849-02.aligned.bam has good EOF block.
So the BAM seems ok. If you use the current version of featureCounts and the error persists you can try and drop a message in the subread google group. If possible try to provide a reproducible example, e.g. a small subset of this bam file, plus which operating system you are on and how you installed the tool. Maybe others here have better suggestions on the actual error.