s4553711
8/29/2016 - 7:29 AM

Note about 23andme script to translate array to VCF

Note about 23andme script to translate array to VCF

rs123 chr1 300 AB
-> no rsid -> rsid = "."
-> AA =~ /-*/  -> next
-> [D,I]B -> next

-> if not find chr1:300 in ref -> log error
-> if find chr1:300 in ref -> stre $ref
-> if $ref == PASS -> next

-> processAltAndGenotype (AB, $ref)
-> if only "A" and A == $ref -> $alt = '.'
-> else
	-> if A != B
		-> if ref == A or B
			-> alt=A or B, GT=0/1
		else
			-> alt=A,B, GT= 1/2
	else
		-> if a == ref
			-> alt=., GT=0/0
		else
			-> alt=A, GT=1/1