How to specify samples for GATK FastaAlternateReferenceMaker in the option --use-iupac-sample
2
1
Entering edit mode
6 months ago
b03601015 ▴ 10

I'm now trying to create a fasta file that contains consensus sequences of multiple samples (for multiple sequence alignment) from a multiple sample vcf file. I tried the following, but it failed to recognize my samples.

gatk FastaAlternateReferenceMaker --output 20210124quinoa_Cq1A.fa -reference /pfs/work7/workspace/scratch/ho_chuang-data-0/ref/quinoa_ref_V2_Cq1A.fa --variant quinoas_Cq1A.vcf.gz --use-iupac-sample S2H1_NKD181100721_H7VJJDSXX_L2,SRR4300210,SRR4300211,SRR4300212,SRR4300213,SRR4300214,SRR4300215 and --use-iupac-sample [S2H1_NKD181100721_H7VJJDSXX_L2,SRR4300210,SRR4300211,SRR4300212,SRR4300213,SRR4300214,SRR4300215] and --use-iupac-sample {S2H1_NKD181100721_H7VJJDSXX_L2,SRR4300210,SRR4300211,SRR4300212,SRR4300213,SRR4300214,SRR4300215} and --use-iupac-sample "S2H1_NKD181100721_H7VJJDSXX_L2,SRR4300210,SRR4300211,SRR4300212,SRR4300213,SRR4300214,SRR4300215"

But none of them worked. How should I list my samples so that the software could recognize it?

SNP genome GATK • 357 views
ADD COMMENT
1
Entering edit mode
6 months ago
Buxus ▴ 10

An alternative is to maybe use bcftools consensus.

ADD COMMENT

Login before adding your answer.

Traffic: 1135 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6