Deleted:featureCounts alias file is not working
0
0
Entering edit mode
vkkodali • 0
@5b1e8323
Last seen 5 weeks ago
United States

I am using v2.0.6 of featureCounts on a Linux machine. I am running the following command:

$ subread-2.0.6-Linux-x86_64/bin/featureCounts -a fc_test.saf -F SAF -o fc_test.counts.out -T 2 -p fc_test.bam

This returns an output where none of the alignments are assigned. That's because I have the following in my SAF vs. BAM:

## SAF
$ head -n2 fc_test.saf
gene-NALCN      NC_086317.1       334712972       334713283       +
gene-NALCN      NC_086317.1       334786484       334786618       +

## BAM
$ samtools view fc_test.bam | head -n 1
SRR5119615.4927416.1    99      NC_086317       334624765       0       73M642546N17M   =       335446699       822024  CCTAAATCCCCTTCCTCAGATGTTGAGGTTCGGACTCTATCAAATATTCTGTACTCTGTCCCTTTTTTTTTTTACTCCCAAATTCAATTT*       NH:i:2  HI:i:1  AS:i:143        nM:i:15 NM:i:11 MD:Z:10T14G4A32G0G0G5A0C1G0G0T13        jM:B:c,1        jI:B:i,334624838,335267383      XS:A:+  MC:Z:90M

My BAM and SAF files use different seq-ids for the chromosome. So, I included the following alias file:

$ cat aliases.txt 
NC_086317.1     NC_086317

And run the following command:

$ subread-2.0.6-Linux-x86_64/bin/featureCounts -a fc_test.saf -F SAF -o fc_test.counts.out -T 2 -p -A aliases.txt fc_test.bam

Still, I see that none of the alignments have been assigned. On the other hand, if I modify my SAF file to remove the version number from the chromosome seq-id, I see ~25% of the alignments assigned.

I think this is a bug in featureCounts and present at least in an older version v2.0.2 and perhaps even earlier. I read in an older post that this could be because of line endings. In my case, all line endings are unix-style.

featureCounts subread Rsubread • 103 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 553 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