The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!
No common samples in VCF and BAM headers, so nothing could possibly be phased!
I want to phase some DNA-seq data.
java -jar GenomeAnalysisTK.jar -T ReadBackedPhasing -R ref.fasta -I readnames.bam --variant test.vcf -L Chr.list -o phased_SNPs.vcf --phaseQualityThresh 20.0
My vcf file looks like this and only contains information for 1 sample
INFO=<ID=SAF,Number=.,Type=Float,Description="Specific Allele Frequency">
FILTER=<ID=DP,Description="Minimum depth of 10">
FILTER=<ID=SAF,Description="Allele frequency of at least 0.3 with base quality minimum 0">
CHROM POS ID REF ALT QUAL FILTER INFO
NC_024331.1 131 . G GA . PASS SAF=0.655738;DP=61
NC_024331.1 147 . C G . PASS SAF=0.320000;DP=25
NC_024331.1 422 . C A . PASS SAF=0.414545;DP=275
I previously had an error message saying my bam file did not have read names. I ran
java -jar AddOrReplaceReadGroups.jar I=sorted.bam O=readnames.bam RGLB=LaneX RGPU=NONE RGSM=AnySampleName RGPL=illumina
Now I am getting an error
ERROR MESSAGE: No common samples in VCF and BAM headers, so nothing could possibly be phased!
Is there somewhere in the header of the vcf I can add AnySampleName?