pLM-BLAST icon indicating copy to clipboard operation
pLM-BLAST copied to clipboard

Redundant alignment shift

Open KYQiu21 opened this issue 1 year ago • 2 comments

Dear developers,

Thanks very much for this excellent software. I encountered these redundant local alignments after postprocessing (see the figure). These alignments are several-residue shifted with each other, and i wonder how to deal with this.

This case is run using the same code as the "advanced" example, and this sequence is: AVEPKEDTITVTAAPAPQESAWGPAATIAARQSATGTKTDTPIQKVPQSISVVTAEEMALHQPKSVKEALSYTPGVSVGTRGASNTYDHLIIRGFAAEGQSQNNYLNGLKLQGNFYNDAVIDPYMLERAEIMRGPVSVLYGKSSPGGLLNMVSKRPTTEPLKEVQFKAGTDSLFQTGFDFSDSLDDDGVYSYRLTGLARSANAQQKGSEEQRYAIAPAFTWRPDDKTNFTFLSYFQNEPETGYYGWLPKEGTVEPLPNGKRLPTDFNEGAKNNTYSRNEKMVGYSFDHEFNDTFTVRQNLRFAENKTSQNSVYGYGVCSDPANAYSKQCAALAPADKGHYLARKYVVDDEKLQNFSVDTQLQSKFATGDIDHTLLTGVDFMRMRNDINAWFGYDDSVPLLNLYNPSSHHHHHHGSSVNTDFDFNAKDPANSGPYRILNKQKQTGVYVQDQAQWDKVLVTLGGRYDWADQESLNRVAGTTDKRDDKQFTWRGGVNYLFDNGVTPYFSYSESFEPSSQVGKDGNIFAPSKGKQYEVGVKYVPEDRPIVVTGAVYNLTKTNNLMADPEGSFFSVEGGEIRARGVEIEAKAALSASVNVVGSYTYTDAEYTTDTTYKGNTPAQVPKHMASLWADYTFFDGPLSGLTLGTGGRYTGSSYGDPANSFKVGSYTVVDALVRYDLARVGMAGSNVALHVNNLFDREYVASCFNTYGCFWGAERQVVATATFRF

image

KYQiu21 avatar Nov 09 '23 10:11 KYQiu21

Now plmblast.py supports additional argument bfactor X, where X is integer > 0. The bigger bfactor the more sparse search - thus there should be less highly overlapping hits.

Argusmocny avatar Dec 14 '23 10:12 Argusmocny

@Argusmocny I would like to ask if there's an advisable scale for the bfactor? how big should it be (ballpark estimate) to minimize the overlapping hits or speed up the alignment search. Thank you!

papelypluma avatar Dec 24 '23 16:12 papelypluma