Skip to content

Latest commit

 

History

History
16 lines (9 loc) · 594 Bytes

README.md

File metadata and controls

16 lines (9 loc) · 594 Bytes

Order_VCF_from_Fasta

this small script allows you to order a vcf file based on the order of the reference fasta file

usage is pretty straightforward:

bash order_vcf_ref.sh IN.vcf fastaRef.fa OUT.vcf

IN.vcf -> input vcf file that you want to order

fastaRef.fa -> your reference in fasta format, i.e. hg19.fa

OUT.vcf -> name for the output file

NOTE: the vcf file and the reference file MUST have the same names for chromosomes, i.e. chromosome 1 shoule be either "chr1" in both or "1" in both

this is a work in progress, I need to add errors check for the input arguments