Skip to content

Commit 0afc477

Browse files
author
Wenzheng Li
authored
Update README.md
1 parent 01d3554 commit 0afc477

File tree

1 file changed

+12
-1
lines changed

1 file changed

+12
-1
lines changed

README.md

Lines changed: 12 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -46,8 +46,19 @@ determined based on what is required in the next script.
4646
In the next step we use an adaptation of a script originally due to
4747
ASDF.
4848
```
49-
./cnvAnalysis.R bin_counts.bed SampelName bins_5k_hg19_gc.txt bins_5k_hg19_exclude.txt
49+
./cnvAnalysis.R bin_counts.bed SampleName bins_5k_hg19_gc.txt bins_5k_hg19_exclude.txt
5050
```
51+
The input file `bin_counts.bed` is the same as described above. The input file
52+
`bins_5k_hg19_gc.txt` is the GC content of each bin. The input `bins_5k_hg19_exclude.txt`
53+
is used to exclude certain parts of the genome that attract an unusual amount of reads.
54+
The format is simply the line numbers, in the corresponding bed file, of the bins
55+
to exclude from the CNV analysis. The first output is a PDF file
56+
`{SampleName}.5k.wg.nobad.pdf` for the CNV profile. In addition,
57+
two tables are saved: one table
58+
`{SampleName}.hg19.5k.nobad.varbin.data.txt` with the information
59+
(chromosome, genome position, GC content, bin count, segmented value) for each bin, and
60+
the other table `{SampleName}.hg19.5k.nobad.varbin.short.txt`
61+
summerizing the breakpoints in the CNV profile.
5162

5263
## Simulating SMURF-seq reads for evaluating mappers
5364

0 commit comments

Comments
 (0)