5.3 createMask

Creating a mask BED file

createMask creates a mask in BED format.

5.3.1 Input

Required inputs :

--bam Input_bam_file.bam	Input bam file.
--type type_of_mask	To indicate the type of mask to be created. Options = depth, nonRef, invariant, variant.

Optional inputs :

--pmd Input_PMD.txt	Post-mortem damage parameters (see PMD for generating such a file)
--recal recal.txt	Quality score recalibration file (see recal for further information)
--fasta Input_refrence_genome_file.fasta	Reference genome to calculate percentage of reference bases that are ‘N’ in window. Used when --maxRefN !=1 or type = nonRef

Specific Parameters :

--minDepthForMask integer_value

To create a mask of all sites with depth >= indicated value for which at least one non-ref allele was observed. Used when type = nonRef. Default = 2.

• See Filter parameters to apply specific filters for bases, reads and parsing window setting.

Engine parameters that are common to all tasks can be found here.

5.3.2 Output

*_minDepth2_maxDepth1000000_depthMask.bed.bed	To specify a 0-based bed file with regions to be masked. Default = No region is masked.
*_filterSummary.txt	File containing Filter summary of general filter counts.
*_RGInfo.json	File containing read group info.

5.3.3 Usage Example

#! /bin/bash

# Set atlas path
atlas=$(dirname "$0")/../build/atlas

# Simulate a BAM File
$atlas simulate --logFile simulate.out

# Create mask for depth within [5, 10]
$atlas createMask --type depth --minDepth 5 --maxDepth 10 --bam ATLAS_simulations.bam --logFile createMask.out