-
Notifications
You must be signed in to change notification settings - Fork 1
/
charr_dosaf
46 lines (43 loc) · 960 Bytes
/
charr_dosaf
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
cd /home/tkess/scratch/CharrGDL_WG
module load htslib
module load samtools
#this is with 2pops - PALE
/home/tkess/angsd/angsd \
-anc /home/tkess/scratch/CharrGDL_WG/Charr_Genome/GCA_002910315.2_ASM291031v2_genomic.fasta \
-bam PALEbams.txt \
-out PALE \
-dosaf 1 \
-GL 1 \
-doCounts 1 \
-nThreads 20 \
-minMapQ 30 \
-minQ 20 \
-minInd 10 \
-doMajorMinor 5 \
-doMaf 2 \
-setMinDepth 10 \
-uniqueOnly 1 \
-remove_bads 1 \
-only_proper_pairs 1 \
-maxDepth 400 \
-rf charr_chroms.txt
#this is with 2pops - DARK
/home/tkess/angsd/angsd \
-anc /home/tkess/scratch/CharrGDL_WG/Charr_Genome/GCA_002910315.2_ASM291031v2_genomic.fasta \
-bam DARKbams.txt \
-out DARK \
-dosaf 1 \
-GL 1 \
-doCounts 1 \
-nThreads 20 \
-minMapQ 30 \
-minQ 20 \
-minInd 10 \
-doMajorMinor 5 \
-doMaf 2 \
-setMinDepth 10 \
-uniqueOnly 1 \
-remove_bads 1 \
-only_proper_pairs 1 \
-maxDepth 400 \
-rf charr_chroms.txt