Commit d331b6b1 authored by eboulanger's avatar eboulanger
Browse files

mini SNP positions update 🐧

parent 177af864
......@@ -3,12 +3,16 @@
## Diplodus
# get CDS positions
grep -v "^##" "diplodus/diplodus_CDS.vcf" | cut -f1-2 | sed '1d' > "diplodus/diplodus_CDS_pos.txt"
# get CDS non-synonymous positions
grep -v "^##" "diplodus/diplodus_CDS_nonsynonymous.vcf" | cut -f1-2 | sed '1d' > "diplodus/diplodus_CDS_nonsynonymous_pos.txt"
# get pcadapt outlier positions
grep -v "^##" "../../seaConnect--dataPrep/04-finalPrep/01-Diplodus/dip_adaptive_413.vcf" | cut -f1-2 | sed '1d' > "diplodus/diplodus_pcadapt_pos.txt"
## Mullus
# get CDS positions
grep -v "^##" "mullus/mullus_CDS.vcf" | cut -f1-2 | sed '1d' > "mullus/mullus_CDS_pos.txt"
# get CDS nonsynonymous positions
grep -v "^##" "mullus/mullus_CDS_nonsynonymous.vcf" | cut -f1-2 | sed '1d' > "mullus/mullus_CDS_nonsynonymous_pos.txt"
# get pcadapt outlier positions
grep -v "^##" "../../seaConnect--dataPrep/04-finalPrep/02-Mullus/mul_adaptive_291.vcf" | cut -f1-2 | sed '1d' > "mullus/mullus_pcadapt_pos.txt"
Markdown is supported
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