Bug dans should-vdj-to-tap.py losqu'il n'y a pas de N décrit et une alternative dans le V
Bug rencontré par @boreec lorsqu'on filtre les germlines sur les 5 meilleurs gènes :
Cette description (dans gosh_igh.should-vdj.fa
) (IGHV3-30*18, IGHV3-30-5*01) IGHD3-22*01 IGHJ2*01
est transformée en cette expression régulière : [[:space:]]*(IGHV3-30[*]18.*|IGHV3-30-5[*]01) [:space:]]*IGHD3-22[*]01.*IGHJ2[*]01
Or il y a un .*
avec le 3-30*18
mais pas après le 3-30-5*01
. On n'ignore donc pas la partie N dans le cas où on a le gène 3-30-5*01. Ce .*
devrait être sorti de l'alternative.