VCF-PHASED-JOIN(1)

πŸ’– vcf-phased-join(1) πŸ’–

okay so bestie, welcome to vcf-phased-join, a lil User Commands icon that lives in your terminal 🫢


✨ NAME ✨

vcf-phased-join β€” takes a whole bunch of overlapping pre-phased chunks and squishes them together into one single VCF πŸ’…


🌸 SYNOPSIS 🌸

vcf-phased-join [OPTIONS] A.vcf B.vcf C.vcf

πŸ¦‹ DESCRIPTION πŸ¦‹

About: the script grabs multiple overlapping pre-phased chunks and concatenates them into one VCF. no cap, that's the whole vibe πŸŽ€

and the way it does it?? it uses the heterozygous calls from the overlaps to figure out the correct phase. like, the overlaps are literally doing the detective work for you πŸ”βœ¨


πŸ’œ OPTIONS πŸ’œ

  • -j, --min-join-quality <num> β€” the quality threshold for gluing those pre-phased blocks together. default is 10 🩹
  • -l, --list <file> β€” a list of VCFs to join. giving very much "here's the squad, bring them all" πŸ“‹
  • -o, --output <file> β€” the output file name. and when you pass "-", it uses STDOUT and STDERR instead πŸ‘‘
  • -q, --min-PQ <num> β€” breaks pre-phased segments if the PQ value is lower in the input VCFs. default is 0.6 πŸ’”
  • -h, -?, --help β€” this help message. your lil safety net when you forget everything πŸ™ˆ

ngl bestie, remember those defaults: 10 for join quality and 0.6 for min-PQ. they're just sitting there waiting for you to not pass a flag πŸ“

vcf-phased-join 🌈 β€” the version from February 2013, straight out of the VCF-PHASED-JOIN(1) era. that's the sign-off, mwah πŸ’‹βœ¨