Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Soft clip disappeared after CallMolecularConsensusReads #1016

Open
jonghwanyoon opened this issue Dec 2, 2024 · 0 comments
Open

Soft clip disappeared after CallMolecularConsensusReads #1016

jonghwanyoon opened this issue Dec 2, 2024 · 0 comments

Comments

@jonghwanyoon
Copy link

hello

When calling Molecular consensus reads using fgbio's CallMolecularConsensusReads, the soft clip disappears.
I went through CorrectUmis, GroupReadsByUmi, CallMolecularConsensusReads, alignment on a bam file mapped with raw fastq, and noticed that the soft clips disappeared on important singleton reads.

Given the presence of soft clips in the bam file resulting from GroupReadsByUmi, it seems that the issue is in CallMolecularConsensusReads.

I only gave CallMolecularConsensusReads the java option and the --min-reads=1 parameter.
There doesn't seem to be any parameter associated with CallMolecularConsensusReads (I could be wrong)
This is the same in 2.1.1, when collapsing between chromosomes became possible, and in 2.4.0, the current latest version.

Although it is a singleton read, it is an important read that represents a clinically meaningful fusion, so I don't want the soft clips to disappear.
May I request your opinion on this issue?

Thank you!

image

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant