blob: 51630876c354b517914e0094ad5a03cbed7c24bb (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
SelectiveGenomeAmplification
============================
PI: http://brisson.bio.upenn.edu/
To use this you'll need:
- A unix environment
- kmer_total_count, a kmer counter available here: http://github.com/mutantturkey/dna-utils/
- bash or compliant shell.
Setup:
git clone git@github.com:mutantturkey/SelectiveGenomeAmplification.git
cd SelectiveGenomeAmplification
make
sudo make install
Example Usage:
selectiveGenomeAmplification PfalciparumGenome.fasta HumanGenome.fasta;
less PfalciparumGenome_HumanGenome/final_mers
For user customizable variables:
max_mer_distance=5000 max_select=6 min_mer_range=6 max_mer_range=12 \
SelectiveGenomeAmplification.sh PfalciparumGenome.fasta half.fasta
## Customizable variables
range of mers, min and max
C | variable | default | notes
:---- | :---- | :---- | ---- | :----
Y | min\_mer\_range | 6 | minimum mer size to use
Y | max\_mer\_range | 12 | maximum mer size to use
Y | max\_mer\_distance | 5000 | maximum distance between mers in foreground
N | output\_directory | $PWD/$foreground\_$background/ | ex. if fg is Bacillus.fasta and bg is HumanGenome.fasta then folder would be $PWD/Bacillus.fasta\_HumanGenome\_output.fasta/
Y | counts\_directory | $output\_directory/.tmp | directory for counts directory
Y | tmp\_directory=$output\_directory/.tmp | temporary files directory
Y | max\_melting\_temp | 30° | maximum melting temp of mers
Y | min\_melting\_temp | 0° | minimum melting temp of mers
Y | min\_foreground\_binding\_average | 50000 | elminate mers that appear less frequently than the average (length of foreground / # of occurances)
Y | max\_select | 15 | maximum number of mers to pick
Y | max\_check | 35 | maximum number of mers to select (check the top #)
Y | ignore\_mers | Not Enabled | mers to explicitly ignore, space seperated ex. ignore\_mers="ACAGTA ACCATAA ATATATAT"
Y | foreground | Not Enabled | path of foreground file
Y | background | Not Enabled | path of background file
Y | max\_consecutive\_binding | 4 | The maxium number of consecutive binding nucleotides in homodimer and heterodimers
|