Genome 540 Homework Assignment 5

Due Sunday Feb. 14, 11:59pm

1. Write a program to identify regions of elevated copy-number using the D-segment algorithm described in lecture, using the number of read-starts at every base to determine the base's score. This program should:
   • read in two input files
     • a mapped read count file (described below)
     • a 'scoring scheme' file that indicates the score for each read count (0, 1, 2, or >=3 reads, scoring scheme described below)
   • output the following information (see template file for example):
     • numbers of background and elevated copy-number segments
     • a list of elevated copy-number segments including their start position, end position, and score (rounded to two decimal points)
     • annotations for the three segments with the highest scores (use the UCSC genome browser and put down anything interesting about the segment)
     • histograms of read-start counts (i.e. number of positions with 0, 1, 2, and >=3 read-starts) for a) background segments, and b) elevated copy-number segments
2. Run your program on this file using the following scoring scheme:
   • score for 0 reads: -0.1077
   • score for 1 read: 0.4772
   • score for 2 reads: 1.0622
   • score for >=3 reads: 1.6748
   • D = -20
   • S = -D = 20
3. The input file has three columns: chromosome, position, and a read start count. The file was created based on the start positions of all reads mapping to chromosome 16 for the individual CHM13. Sequencing was performed on the Illumina platform, and the reads were mapped to the human reference hg38 using BWA. To create the likelihood scores, the number of read starts mapping to chromsome 16 were used to estimate the mean number, m, of read starts per base. This was then used to create two Poisson distributions: one with mean m, and the other with mean 1.5 m. The first distribution should apply to read starts in 'background' (= homozygous reference) while the second distribution should reflect a region heterozygously duplicated in the sample relative to the reference. The scoring scheme is the log likelihood ratio (with log base = 2) of these two distributions given the number of read starts observed. As discussed in lecture, 2^S should correspond (very approximately -- & ignoring the constant K) to the expected average spacing between segments of score >= S in 'background'.
4. Using this example input and this example scoring scheme your output file should look like this template. Use the same template structure for your output on the actual file. Please put everything into ONE plain text file - do not send an archive of files or a tar file, or a word processing document file. Compress it (using either Unix compress, or gzip -- if you don't have access to either of these programs let us know), and send it as an attachment to both Phil (phg (at) uw.edu) and Dani (dfaivre (at) uw.edu).