Commit 67f0bfd8 authored by Marwan ELKREWI's avatar Marwan ELKREWI
Browse files

Update SNP_calling_in_backcrossing_experiment.md

parent a4c89a28
......@@ -125,15 +125,11 @@ samtools mpileup -B Asex.sort.bam Sex.sort.bam > Asex_Sex.mpileup
```
### 9. Run Popoolation
```
java -ea -Xmx7g -jar /Popoolation2/popoolation2_1201/mpileup2sync.jar --input Asex_Sex.mpileup --output Asex_Sex_java.sync --fastq-type sanger --min-qual 20 --threads 2
perl Popoolation2/popoolation2_1201/snp-frequency-diff.pl --input Asex_Sex_java.sync --output-prefix Asex_Sex --min-count 6 --min-coverage 10 --max-coverage 2000
per SNP
#per SNP
perl Popoolation2/popoolation2_1201/fst-sliding.pl --input Asex_Sex_java.sync --output Asex_Sex.persnp.fst --suppress-noninformative --min-count 2 --min-coverage 10 --max-coverage 2000 --min-covered-fraction 1 --window-size 1 --step-size 1 --pool-size 5:10
per 1kb
#per 1kb
perl Popoolation2/popoolation2_1201/fst-sliding.pl --input Asex_Sex_java.sync --output Asex_Sex.1kb.min10.fst --suppress-noninformative --min-count 2 --min-coverage 10 --max-coverage 2000 --min-covered-fraction 1 --window-size 1000 --step-size 1000 --pool-size 5:10
```
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment