From 0a4cc80f9d74db3593e1279185945bcf1d153481 Mon Sep 17 00:00:00 2001 From: Simon Waldherr Date: Sat, 13 Oct 2018 00:05:21 +0200 Subject: [PATCH] link to new, recommended DOI resolver in README.md --- README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/README.md b/README.md index 9e3e1ae..a5a84cb 100644 --- a/README.md +++ b/README.md @@ -107,7 +107,7 @@ Aligning is a computationally intensive process and will take far more CPU time Many of the Hi-C reads can be determined to be noise, rather than signal. It is not strictly necessary to remove these reads from your BAM files before running LACHESIS. But it will reduce the files' size, reduce the I/O time required by LACHESIS to read them in, and may reduce LACHESIS' error rate. -We have included the scripts `PreprocessSAMs.pl` and `PreprocessSAMs.sh`, which perform the filtering methods we used. `PreprocessSAMs.pl` removes all reads that do not map within 500bp of a restriction enzyme site, following the suggestion of [Yaffe & Tanay (Nature Genetics, 2011)](dx.doi.org/10.1038/ng.947). To do this, it uses the script `make_BED_around_RE_site.pl` (also included) and also the bedtools library. `PreprocessSAMs.pl` also removes unpaired reads, which LACHESIS cannot use. Make sure to set the variable `$RE_site` in `PreprocessSAMs.pl`. You can also use the script `PreprocessSAMs.sh` as a batch script to run `PreprocessSAMs.pl` on a set of SAM/BAM files (set the variables `SAMs` and `ASSEMBLY`). The final set of SAM/BAM files is specified in the INI file as the parameter `SAM_FILES`, which names the files, and `SAM_DIR`, which names the directory they are in. +We have included the scripts `PreprocessSAMs.pl` and `PreprocessSAMs.sh`, which perform the filtering methods we used. `PreprocessSAMs.pl` removes all reads that do not map within 500bp of a restriction enzyme site, following the suggestion of [Yaffe & Tanay (Nature Genetics, 2011)](https://doi.org/10.1038/ng.947). To do this, it uses the script `make_BED_around_RE_site.pl` (also included) and also the bedtools library. `PreprocessSAMs.pl` also removes unpaired reads, which LACHESIS cannot use. Make sure to set the variable `$RE_site` in `PreprocessSAMs.pl`. You can also use the script `PreprocessSAMs.sh` as a batch script to run `PreprocessSAMs.pl` on a set of SAM/BAM files (set the variables `SAMs` and `ASSEMBLY`). The final set of SAM/BAM files is specified in the INI file as the parameter `SAM_FILES`, which names the files, and `SAM_DIR`, which names the directory they are in. #### 5. Aligning the draft assembly to the reference assembly, if there is one