Reference Code backup Executable files
Provided a BAM/SAM file and reference gene model, this module will calculate how mapped reads were distributed over genome feature (like CDS exon, 5’UTR exon, 3’ UTR exon, Intron, Intergenic regions). When genome features are overlapped (e.g. a region could be annotated as both exon and intron by two different transcripts) , they are prioritize as: CDS exons > UTR exons > Introns > Intergenic regions, for example, if a read was mapped to both CDS exon and intron, it will be assigned to CDS exons.
read_distribution.py -i Pairend_StrandSpecific_51mer_Human_hg19.bam -r hg19.refseq.bed12
--version | show program’s version number and exit |
-h, --help | show this help message and exit |
-i INPUT_FILE, --input-file=INPUT_FILE | |
Alignment file in BAM or SAM format. | |
-r REF_GENE_MODEL, --refgene=REF_GENE_MODEL | |
Reference gene model in bed format. |