Skip to content

Instantly share code, notes, and snippets.

What would you like to do?
# variant_demo_simplify.tsv hold basic info of variant, looks like:
# chr1 69511 . A G
# chr1 865694 . C T
# 1. split input tsv file, cause if records are too many, a server error will be returned,
## although offical site stat that up to 50000 VCF rows may be posted one time
# 2. use NCBI variation services API to post all rows which are stored in a file
split -l 10000 /path/to/variant_demo_simplify.tsv
for s in $(echo "xa?"); do
curl -X POST \
--header "Content-Type: text/plain; charset=utf-8" --header "Accept: text/plain; charset=utf-8" \
--data-binary "@${s}" \
"" \
> variant_demo_${s}_recalibrated.tsv
# Here I'm using hg19, if you are using hg38, you should use another assembly version in API url
# If add "-i" a header will appear in the returned result
# For more info, please refer to offical website at:
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.