Mouse VDJ Data

[zktuong]

Discussed in #310

^{Originally posted by hungms September 5, 2023}
Hi Kelvin,

Thanks for developing this great package.

I am currently analysing mouse VDJ data with your package in the singularity container but I see all my IGHA calls from igblast were assigned as IGHA2 after correction.

2MpOG_1A.pdf

I see that in Step 4: Assigning constant region calls, you mentioned:

Further, to deal with incorrect constant gene calls due to insufficient length, a pairwise alignment will be run against curated sequences that were deemed to be highly specific in distinguishing IGHA1 vs IGHA2, and IGHG1 to IGHG4.

Can I check if the IGHA2 calls coming from this step (the curated sequences are from human I presume?). If this is the case, can I presume that it won't be appropriate to run this step for mouse VDJ to correct the other constant region calls? Or is there an alternative approach to this?

Thank you again in advance.

need to create a step in the singularity container as well to switch the correction_dict, or disable correction.