sortmerna | generation sequence filtering and alignment tool | Genomics library
kandi X-RAY | sortmerna Summary
kandi X-RAY | sortmerna Summary
SortMeRNA is a local sequence alignment tool for filtering, mapping and clustering. The core algorithm is based on approximate seeds and allows for sensitive analysis of NGS reads. The main application of SortMeRNA is filtering rRNA from metatranscriptomic data. SortMeRNA takes as input files of reads (fasta, fastq, fasta.gz, fastq.gz) and one or multiple rRNA database file(s), and sorts apart aligned and rejected reads into two files. Additional applications include clustering and taxonomy assignation available through QIIME v1.9.1. SortMeRNA works with Illumina, Ion Torrent and PacBio data, and can produce SAM and BLAST-like alignments.
Support
Quality
Security
License
Reuse
Top functions reviewed by kandi - BETA
Currently covering the most popular Java, JavaScript and Python libraries. See a Sample of sortmerna
sortmerna Key Features
sortmerna Examples and Code Snippets
Community Discussions
Trending Discussions on sortmerna
QUESTION
I wrote a snakemake pipeline for performing sortmeRNA version 4.2.0. The pipeline is as follows, and works perfectly when I run it for 1 sample:
...ANSWER
Answered 2020-Sep-14 at 09:08The problem here is that {SAMPLES}
in your inputs and outputs is a wildcard. In the shell it is read as the global variable defined above the rule.
You should use {wildcards.SAMPLES}
in the shell section:
QUESTION
I know this is a common error and I already checked other posts but it didn't resolved my issue. I would like to use the name of the database I use for SortMeRNA
rule (rRNAdb=config["rRNA_database"]
) the same way I use version=config["genome_version"]
. But obviously, I can't.
ANSWER
Answered 2020-Jan-21 at 16:26All output
files need to have same wildcards, or else it would cause conflict in resolving job dependencies. Not all files in output:
have {rRNAdb}
wildcard, which is causing this problem. For example, if you have two {rRNAdb}
values, both would write to file "{OUTDIR}/temp/{sample}.fastq"
, which snakemake correctly doesn't allow.
QUESTION
I tried to create a snakemake file to run sortmeRNA pipeline:
...ANSWER
Answered 2020-Jan-16 at 13:46For Snakemake to connect the input of one rule to the output of another, they will need to be identical. The way you describe the output of sortmeRNA
and the input unmerge_paired
do not work together, whether you put temp()
around it or not.
Community Discussions, Code Snippets contain sources that include Stack Exchange Network
Vulnerabilities
No vulnerabilities reported
Install sortmerna
Support
Reuse Trending Solutions
Find, review, and download reusable Libraries, Code Snippets, Cloud APIs from over 650 million Knowledge Items
Find more librariesStay Updated
Subscribe to our newsletter for trending solutions and developer bootcamps
Share this Page