Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
seaconnect
ddocent_radseq
Commits
32dd907f
Commit
32dd907f
authored
Aug 13, 2019
by
peguerin
Browse files
Final.recode.vcf correction
parent
358df771
Changes
2
Hide whitespace changes
Inline
Side-by-side
00-scripts/rm_indels.sh
View file @
32dd907f
...
...
@@ -3,4 +3,4 @@
SPECIES
=
$1
CONTAINER
=
$2
singularity
exec
"
${
CONTAINER
}
"
vcftools
--vcf
04-ddocent/
"
${
SPECIES
}
"
/TotalRawSNPs.vcf
--remove-indels
--recode
--recode-INFO-all
--out
05-vcf/
"
${
SPECIES
}
"
\ No newline at end of file
singularity
exec
"
${
CONTAINER
}
"
vcftools
--vcf
04-ddocent/
"
${
SPECIES
}
"
/Final.recode.vcf
--remove-indels
--recode
--recode-INFO-all
--out
05-vcf/
"
${
SPECIES
}
"
\ No newline at end of file
README.md
View file @
32dd907f
...
...
@@ -206,7 +206,7 @@ Genotypes are stored as `VCF` files. We keep only SNPs (Single Nucleotide Polymo
SPECIES=mullus
CONTAINER=/entrepot/working/seaconnect/seaConnect--dDocent/seaconnect.simg
## filter VCF remove INDELs
singularity exec "${CONTAINER}" vcftools --vcf 04-ddocent/"${SPECIES}"/
TotalRawSNPs
.vcf --remove-indels --recode --recode-INFO-all --out 05-vcf/"${SPECIES}"
singularity exec "${CONTAINER}" vcftools --vcf 04-ddocent/"${SPECIES}"/
Final.recode
.vcf --remove-indels --recode --recode-INFO-all --out 05-vcf/"${SPECIES}"
```
Final result is generated as
`VCF`
files into
[
05-vcf
](
05-vcf
)
folder and contains only SNPs.
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment