mixcr
mixcr copied to clipboard
Align D gene even when V or J is absent
@mikessh
Actually the full list of what is needed (i.e. observed in lymphoma samples):
- Vd-Dd3, Dd2-Jd
- Dd-Ja
- Dh-Jh
- More complex: Dd2-Dd3, Vk-KDE, INTRON-KDE (IGK)
here d,a,h,k
stands for TRD,TRA,IGH,IGK
To speed up it is reasonable to suggest that best D match cannot be more than 30 bases from the available reference point
Requires additional threshold in parameters for minimal score for incomplete alignments. Floating bounds ? Once again, is it possible to add several separate D gene hits into our abstraction? Stas?