Find Resistance with ShortBRED
This tool allows you to detect and quantify the presence of antibiotic resistance (AR) marker genes by running DIAMOND with a database of peptide marker sequences which represent the genes of interest. The Find Resistance with ShortBRED tool works similarly to the quantify step of ShortBRED, a public bioinformatics pipeline and resource. To learn more about the ShortBRED-Quantify tool [Kaminski et al., 2015], see https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1004557.
Find Resistance with ShortBRED quantifies the presence of Antibiotic Resistance (AR) marker genes in a sample of NGS short reads. It is possible to output a sequence list containing all the input reads which contained a marker (each read in the output is annotated with metadata describing the properties of the marker detected in the read).
Antibiotic Resistance marker databases for use with the tool can be downloaded using Download Resistance Database (Download Resistance Database).
To start the tool, go to:
Toolbox | Microbial Genomics Module (
) | Drug Resistance Analysis (
) | Find Resistance with ShortBRED (
).
The tool accepts a nucleotide sequence or sequence list as input.
In the next dialog (figure 12.5), several parameters are available:
Figure 12.5: References and search parameters.
- Reference marker database: select the antimicrobial resistance marker database.
- Genetic code: select the genetic code to use when DIAMOND translates the nucleotide sequences to proteins.
- E-value: specify an expectation value to use as threshold for qualifying hits with DIAMOND.
- Sensitivity: select DIAMOND sensitivity setting.
- Percent identity: defines a minimum threshold for the percent identity of an alignment. This value is used by Find Resistance with ShortBRED to determine whether a hit found by DIAMOND is sufficiently good to be validated as a true hit (equivalent to the parameter "-id" of the ShortBRED-Quantify tool).
- Minimum alignment length: the minimum length of the DIAMOND alignment. This value is used by Find Resistance with ShortBRED to determine whether a hit found by DIAMOND is sufficiently good to be validated as a true hit (equivalent to the parameter "-pctlength" of the ShortBRED-Quantify tool - see citation at the beginning of this section).
- Minimum read length: the minimum read length. This value is used by Find Resistance with ShortBRED to determine whether a read is long enough to be processed by Find Resistance with ShortBRED (equivalent to the parameter "-minreadBP" of the ShortBRED-Quantify tool).
The Find Resistance with ShortBRED tool will output a resistance abundance table, a result summary table, an optional report, and an optional sequence list.
The result summary table provides an overview of all the resistance phenotypes in the applied database and reports the number of identified resistance genes and markers, and number of reads assigned to each phenotype.
The optional report output contains general information about the input sample, the marker database used and a short result summary.
The optional sequence list output contains all reads from the input sample which were found to contain one of the AR marker sequences. Each read in the sequence list is annotated with metadata describing the properties of the marker detected in the read.
Subsections