modules/quilt_quilt

Groovy Map containing sample information
e.g. [ id:‘test’, single_end

(Mandatory) BAM/CRAM files

(Mandatory) BAM/CRAM index files

(Mandatory) “Path to file with bam file locations. File is one row per entry, path to bam files. Bam index files should exist in same directory as for each bam, suffixed either .bam.bai or .bai.

(Mandatory) Reference haplotype file in IMPUTE format (file with no header and no rownames, one row per SNP, one column per reference haplotype, space separated, values must be 0 or 1)

(Mandatory) Reference haplotype legend file in IMPUTE format (file with one row per SNP, and a header including position for the physical position in 1 based coordinates, a0 for the reference allele, and a1 for the alternate allele).

(Mandatory) What chromosome to run. Should match BAM headers.

(Mandatory) When running imputation, where to start from. The 1-based position x is kept if regionStart <= x <= regionEnd.

(Mandatory) When running imputation, where to stop.

Buffer of region to perform imputation over. So imputation is run form regionStart-buffer to regionEnd+buffer, and reported for regionStart to regionEnd, including the bases of regionStart and regionEnd.

Number of generations since founding or mixing. Note that the algorithm is relatively robust to this. Use nGen = 4 *.Ne / K if unsure.

(Optional) File with genetic map information, a file with 3 white-space delimited entries giving position (1-based), genetic rate map in cM/Mbp, and genetic map in cM. If no file included, rate is based on physical distance and expected rate (expRate).

Groovy Map containing sample information
e.g. [ id:‘test’, single_end

(Optional) File with positions of where to impute, lining up one-to-one with genfile. File is tab seperated with no header, one row per SNP, with col 1 = chromosome, col 2 = physical position (sorted from smallest to largest), col 3 = reference base, col 4 = alternate base. Bases are capitalized.

(Optional) File with truth phasing results. Supersedes genfile if both options given. File has a header row with a name for each sample, matching what is found in the bam file. Each subject is then a tab seperated column, with 0 = ref and 1 = alt, separated by a vertical bar |, e.g. 0|0 or 0|1. Note therefore this file has one more row than posfile which has no header.

Groovy Map containing sample information
e.g. [ id:‘test’, single_end

(Optional) File with reference genome.