search

nlesc-ave / ave-rest-service

visualize (clustered) single-nucleotide variants across genomes
Apache License 2.0
1 stars 0 forks source link

How to add an accession to an existing bcf file? #29

Open sverhoeven opened 7 years ago

sverhoeven commented 7 years ago

what if you have another acc resequenced that you would like to add to AVE2 ? Do you have feed it another vcf file that has info on all accessions, or just a single vcf

sverhoeven commented 7 years ago

There is bcftools merge -O b -o bigwithnewsample.bcf big.bcf newsample.bcf command to merge bcf files. After the merge the resulting bcf file should be indexed with bcftools index bigwithnewsample.bcf. The newsample.bcf file should created from a vcf, by sorting (vcf-sort -c newsample.vcf > newsample.sorted.vcf) and converting (bcftools view -O b newsample.sorted.vcf > newsample.bcf).

These commands should be tested.