Skip to content
Snippets Groups Projects
Commit e9b3ba36 authored by alice.boizet_cirad.fr's avatar alice.boizet_cirad.fr
Browse files

return 0/1 instead of AB in vcf

parent a79543cc
No related branches found
No related tags found
No related merge requests found
...@@ -109,7 +109,7 @@ def generate_chr_for_one_ind (chromosome_size, marker_density, err_rate, mean_de ...@@ -109,7 +109,7 @@ def generate_chr_for_one_ind (chromosome_size, marker_density, err_rate, mean_de
#print("reads",genotype1, genotype2) #print("reads",genotype1, genotype2)
#print(genotype) #print(genotype)
# Append the genotype and the depth of the current marker to the segment list # Append the genotype and the depth of the current marker to the segment list
finalGenotype = str(genotype1) + str(genotype2) + ":" + str(site_depth) + ":" + str(x) + "," + str(y) + ":.:.:.:.:.:." finalGenotype = str(genotype) + ":" + str(site_depth) + ":" + str(x) + "," + str(y) + ":.:.:.:.:.:."
segment.append(finalGenotype) segment.append(finalGenotype)
return segment return segment
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment