The featurecount function currently only supports SAM and GTF files. The only workaround is to convert your files to these formats.
Using a BioMap object is currently the most straightforward way to convert from SAM to BAM using the Bioinformatics Toolbox.
b = BioMap('test.bam');
write(b,'test.sam','format','sam')
A direct conversion from GFF3 to GTF may not be possible without some simplifying assumptions, and there is currently no way to do this using the Bioinformatics Toolbox. Since GTF is identical to GFF2, you may be able to give the filename of a GFF2 file directly to featurecount.
If you have a third party conversion tool, you can call this tool from within MATLAB by using the ! character, as in this example.
if exist('test.sam') ~= 2
!thirdpartyconversiontool test.bam test.sam
end
if exist('test.gtf') ~= 2
!thirdpartyconversiontool test.gff test.gtf
end
Best Answer