RNA-seq / Count aligned reads per exons for DEXSeq using own GTF

Description

Given mapped reads in a BAM file and genomic annotation in a GTF file, this tool counts the reads that fall into each non-overlapping exonic part using the script dexseq-count.py.

Parameters

• Does the alignment file contain paired-end data (yes, no) [no]
• Chromosome names in the BAM file look like (chr1, 1) [1]
• Does the BAM file contain paired-end data (yes, no) [no]
• Was the data produced with a strand-specific protocol (yes, no, reverse) [no]
• Mode to handle reads overlapping more than one feature (union, intersection-strict, intersection-nonempty) [union]

Details

Indicate if your data was not produced with a strand-specific RNA-seq protocol, so that a read is considered overlapping with a feature regardless of whether it is mapped to the same or the opposite strand as the feature. Otherwise the read has to be mapped to the same strand as the feature.

Output

Output is a table with counts for each non-overlapping exonic part. In order to use the output in DEXSeq, you need to select all samples and run the tool "Utilities - Define NGS experiment".

References

This tool is based on the HTSeq package by Simon Anders and it is available in the DEXSeq Bioconductor package.