This app takes an array of fastq files and runs amplivar.
Amplivar creates a list of unique reads and then performs alignment using Blat.
Analysis of cancer samples
- A indexed reference genome (
.tar.gz
format) - A list of fastq files.
- A list of usual suspects
- A list of amplicon flanking sequences (by default this is the file saved in 001_ToolsReferenceData)
- Value for minimum reported variant frequency as percentage of amplicon reads (optional - default = 1 (1%). If left blank uses amplivar default (5%))
- Value for minimum coverage for variant calling (optional - leaving blank uses amplivar default of 10 (default))
- Value for minimum number reads containing the variant allele (optional - leaving blank uses amplivar default of 5 (default))
- The branch of amplivar_blat to use. The master branch should be used for all production runs, and this is set as the default value. However this can be overwritten by settings in an Nexus workflow therefore this should always be specified in the dx run command.
- A directory of bam files (in a folder called bam)
- A directory of bam files indexes
- A directory containing the read depth for each amplicon (flanked file)
- One folder per sample including all other files created during analysis
- Clone the amplivar github repo and checks out the specified branch
- Downloads and unzips reference genome
- If required makes 2 bit reference file.
- Run BLAT server
- Test the BLAT server is working using a example fasta file
- converts the encoding of the flanking file using dos2unix, as this can cause the app to fail
- Run amplivar using -m Variant_Calling and arguments defined in inputs
- To speed up the upload of the many files produced by amplivar these are zipped on a sample by sample basis.
- This app uses a copy of Amplivar modified to exclude variant calling (https://github.com/moka-guys/amplivar_blat)
- The app was made by Viapath Genome Informatics