Currently, MakeDb.py only takes the first V gene assigned by igblastp when using the igblast-aa command. When there are mutltiple alleles with the top score, all alleles should ideally be reported for completeness and consistency with the igblast nt AIRR output.
Example of one sequence with two identical-scoring alleles:
fasta igblast input: example_igblastp.txt
igblast output: example_input_fasta.txt
current makedb output: example_makedb_airr.tsv
desired makedb output: desired_output_makedb_airr.tsv
Currently, MakeDb.py only takes the first V gene assigned by igblastp when using the igblast-aa command. When there are mutltiple alleles with the top score, all alleles should ideally be reported for completeness and consistency with the igblast nt AIRR output.
Example of one sequence with two identical-scoring alleles:
fasta igblast input: example_igblastp.txt
igblast output: example_input_fasta.txt
current makedb output: example_makedb_airr.tsv
desired makedb output: desired_output_makedb_airr.tsv