Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

FastaAlternateReferenceMaker generate a reference fasta rather than my sample haplotype fasta

$
0
0

Hi, I am using the following command, wanting to generate haplotype fasta for 97 samples:
java -Xmx25g -jar $gatk -T FastaAlternateReferenceMaker -R $ref_genome -o new_output.fasta --variant $vcf

Here are what I did to my .vcf file:
1) I keep all the meta-information
2) I only select a specific region in my .vcf file (around SNP sites within 3Mb in the reference genome)

I attached a .zip file containing the first 120 lines of my .vcf file, I do not know how do deal with this error.
Looking forward to professional guidance, thank you!


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>