Usage: ./slncky.v1.0 bedfile assembly out_prefix
./slncky.v1.0 transcripts.bed mm9 out
./slncky.v1.0 -2 transcripts.bed mm9 out
.
./slncky.v1.0 -w transcripts.bed mm9 out
annotations.config
file in the same directory as itself, but you can set this path if you want:
./slncky.v1.0 -c ./annotations.config transcripts.bed mm9 out
./slncky.v1.0 --bedtools_path /unix/home/jjenny/bedtools-2.17.0/bin transcripts.bed mm9 out
Customize slncky | |
General Options | |
-c / --config | Specify path of annotations.config file. (Default:annotations.config assumed to be in same directory as ./slncky executable.) |
-1 / --no_orth_search | Only filter; don't search for orthologs. |
-2 / --no_filter | Don't filter; only search for orthologs. |
-o / --overwrite | Force overwrite of trancsripts.lncs.bed. |
-n / --threads | Number of threads to use for aligning transcripts. (Default: 5) |
-w / --web | Create browser to visualize aligned lncRNAs. |
--bedtools | Path to bedtools executables. (Default: bedtools assumed to be in user's PATH) |
--liftover | Path to liftover executable. (Default: bedtools assumed to be in user's PATH) |
--lastz | Path to lastz executable. (Default: bedtools assumed to be in user's PATH) |
Filter Options | |
--min_overlap | Minimum overlap for filtering transcripts overlapping CODING or MAPPED_CODING gene. (Default: 0.00) |
--min_cluster | Minimum size of duplication cluster for filtering paralogous transcripts. (Default: 2) |
--min_coding | Minimum identity for filtering transcripts that align to syntenic CODING gene. (Default: slncky learns threshold from input data.) |
--no_collapse | Do not collapse overlapping transcripts into single canonical gene. |
--no_overlap | Do not remove transcripts that overlap CODING or MAPPED_CODING gene. |
--no_self | Do not remove transcripts sharing sequence similarity with each other. |
--no_dup | Do not remove transcripts that share sequence similarity with DUP genes. |
--no_coding | Do not remove transcripts that align to syntenic CODING genes. |
Orthology Search Options | |
--min_noncoding | Minimum overlap for reporting transcripts that align to syntenic NONCODING gene. (Default: slncky learns threshold by aligning to shuffled intergenic regions) |
--no_bg | Do not align to shuffled background; report all alignments. Setting this option will greatly reduce slncky runtime. |
--no_orf | Do not check for conserved ORFs in aligned noncoding transcripts. |
Advanced orthology search options | |
--minMatch | minMatch parameter to pass to liftOver. Decreasing this parameter will increase sensitivity of orthology search. (Default: 0.1) |
--pad | Number of basepairs up- and downstream of lncRNA loci to include in orthology search. Increasing this parameter will increase sensitivity of orthology search. (Default: 50,000) |
--gap_open | Gap open penalty to pass to lastz. Decreasing this parameter will increase sensitivity of orthology search. (Default: 1.0) |
--gap_extend | Gap extend penalty to pass to lastz. Decreasing this parameter will increase sensitivity of orthology search. (Default: 40) |
out.lncs.bed
out.lncs.info.txt
out.filtered_info.txt
/ out.cluster_info.txt
out.canonical_to_lncs.txt
out.orthologs.txt
/ out.orthologs.top.txt
out.orfs.txt
out_browser/
--web
flag was set, slncky will create an evolutionary browser that is viewable from out_browser/browse.html