I use STAR for FUSION detect. But some times the fusion partner's sequence is similar to the next exon(in gtf file). Mapping to fusion partner lead to 0 mismatches, and mapping to next exon lead to 5 mismatches in 16bp(this is very low quality).
However STAR will map to next exon(15bp), then softclip remaining bases, and no chimeric alignment output, lead to fusion false negative.
If I increase alignSJDBoverhangMin to 18bp, I will get a correct chimeric alignment. But alignSJDBoverhangMin=18 is too big, lead to a lot of chimeric alignment between adjacent exons.
I think this can be avoid by increase the penalty for a mismatch. But STAR do not have a option to set it. Can I have a way to adjust it(e.g. change the code if it is written as a const value)?
I use STAR for FUSION detect. But some times the fusion partner's sequence is similar to the next exon(in gtf file). Mapping to fusion partner lead to 0 mismatches, and mapping to next exon lead to 5 mismatches in 16bp(this is very low quality).
However STAR will map to next exon(15bp), then softclip remaining bases, and no chimeric alignment output, lead to fusion false negative.
If I increase alignSJDBoverhangMin to 18bp, I will get a correct chimeric alignment. But alignSJDBoverhangMin=18 is too big, lead to a lot of chimeric alignment between adjacent exons.
I think this can be avoid by increase the penalty for a mismatch. But STAR do not have a option to set it. Can I have a way to adjust it(e.g. change the code if it is written as a const value)?