Permalink
Cannot retrieve contributors at this time
Name already in use
A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?
MetaGeneView/README.md
Go to fileThis commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
43 lines (33 sloc)
1.23 KB
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# MetaGeneView | |
create a Genome Browser Track of linearised genomic coverage (removing introns). | |
## create coverage from BAM files | |
The script scales the read to represent the coverage as reads per million [RPMs]. | |
``` | |
perl Batch | |
bash BatchCoverage.sh RiboData/bams/MonoVsPoly/ _genome_umi RiboData/gbed/MonoVsPoly/ | |
``` | |
## create a track file | |
``` | |
perl MetaGeneView.pl --bed /path/to/annotation/ncbiRefSeq_rn6_Sep2018.bed --reps 3 --gbed $(ls /path/to/gbed/*.gbed.gz) --query "Camk2a" > myCustomTrack.txt | |
``` | |
* note: use double quotes to escape special characted in query gene symbol | |
## help | |
``` | |
$ perl MetaGeneView.pl -h | |
MetageneView Version 1.0 | |
usage: MetaGeneView.pl --bed annotation.bed --gbed runA.gbed.gz runB.gbed.gz --query GeneOfInterest | |
description: MetageneView creates a BedGraph track file with linear gene representation (no introns). | |
parameters: | |
-q|--query | |
query gene name | |
-r|--reps | |
number of replica, to accumulate coverage in consecutive files | |
-b|--bed | |
valid annotation file in BED12 format | |
-g|--gbed | |
single or a space-separated list of GraphBed files. | |
The GraphBed files need to be block compressed with BGZIP. | |
The script uses Tabix to query the coverage files. | |
-help | |
define usage | |
``` |