diff --git a/R/io.R b/R/io.R
index 77c1cf2..ab7c658 100644
--- a/R/io.R
+++ b/R/io.R
@@ -1203,7 +1203,7 @@ fn_load_genotype = function(fname_geno, ploidy=NULL, force_biallelic=TRUE, retai
                             "The first 3 columns do not correspond to 'chr', 'pos', and 'allele'."))
                     return(error)
                 }
-                vec_loci_names = paste(df[,1], df[,2], df[,3], df[,4], sep="\t")
+                vec_loci_names = paste(df[,1], df[,2], df[,4], sep="\t")
                 # TODO: enable below as the default so time in the future
                 # vec_loci_names = paste(df[,1], df[,2], df[,3], df[,4], sep="\t")
                 vec_entries = colnames(df)[c(-1:-4)]
diff --git a/man/fn_load_genotype.Rd b/man/fn_load_genotype.Rd
index 8eb4b7d..1eb0507 100644
--- a/man/fn_load_genotype.Rd
+++ b/man/fn_load_genotype.Rd
@@ -12,6 +12,7 @@ fn_load_genotype(
   retain_minus_one_alleles_per_locus = TRUE,
   min_depth = 0,
   max_depth = .Machine$integer.max,
+  new_format = FALSE,
   verbose = FALSE
 )
 }
@@ -44,6 +45,8 @@ followed by \code{bcftools -m - ...} to  split mult-iallelic sites into bialleli
 
 \item{max_depth}{if input is a VCF file: maximum depth per locus beyond which will be set to missing data (Default=.Machine$integer.max)}
 
+\item{new_format}{use the new allele frequency table format? (i.e. includes all the alleles between \code{pos} and \code{allele}, e.g. header: \code{all_alleles} with rows: \code{A|T}, and \code{C|G|T}) (Default=FALSE)}
+
 \item{verbose}{show genotype loading messages? (Default=FALSE)}
 }
 \value{