modules/quilt_quilt
QUILT is an R and C++ program for rapid genotype imputation from low-coverage sequence using a large reference panel.
Description
QUILT is an R and C++ program for rapid genotype imputation from low-coverage sequence using a large reference panel.
Input
(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.
*.{txt}
(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)
*.{hap.gz}
(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).
*.{legend.gz}
(Mandatory) When running imputation, where to start from. The 1-based position x is kept if regionStart <= x <= regionEnd.
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).
*.{txt.gz}
(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.
*.{txt}
(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.
*.{txt}