Closed stevekm closed 7 years ago
In at least one case, a sample variant with no coding change somehow passed this step and generated a summary table, while the variant itself still lacked the Amino Acid change field in the annotation (e.g. p.R123H
item). Need to investigate why the lack of this field results in errors for some cases but not others.
This error also occurs if an intronic variant gets through filtering criteria and is processed by the function
In the script
merge_vcf_annotations.py
, functionsplit_df_col2cols
gives an error if the number of fields output bymerge_df['AAChange'].str.split(':')
is fewer than the number of replacement columns (5 given by default).