-
Notifications
You must be signed in to change notification settings - Fork 401
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
Do you have a recommendation for downstream SV calling, Do it work with Sniffles? #142
Comments
See #137, #118 and fritzsedlazeck/Sniffles#51. Please update to the latest version v2.10 and run minimap2 with "--MD". Let me know if that works. Thanks. |
Hi Heng
We generated two reference genome fragments with 1m.fa containing a 1Mb
deletion from 500K to 1500K. We observed the DEL using ngmlr + sniffles,
but minimap + sniffles doesn't work. Can you please try? Thank you!
Cheers,
Yongfeng
minimap2 --MD -ax asm5 2m.fa 1m.fa >t.sam
samtools view -ubS t.sam >t.bam
sniffles -m t.bam -v t.vcf --skip_parameter_estimation --min_support 1
|
We sent the fasta file to your email: lh3@me.com, Thank you! |
My colleague compared ngmlr and minimap2 and found that the XS tag is missing in minimap2 sam files, which lead to the disconnect between minimap2 and Sniffles. Can you please check? Thank you! |
Could you have a look at fritzsedlazeck/Sniffles#58 and try the latest sniffles? |
Great! The most recent Minimap2 + the newest sniffles worked! Thank you! |
It is difficult to ask minimap2 to output the XS tag. Unlike bwa-mem, minimap2 doesn't compute XS for most of alignments. It is possible to emulate XS, but that would be "hacky". It is good to know that the latest sniffles works. Thanks for the confirmation! |
I believe this can be closed. |
I am interested in calling SVs based on minimap2 mapping. Do you have a recommendation for downstream SV calling? I am using Sniffles (https://github.com/fritzsedlazeck/Sniffles) to call SVs based on ngmlr mapping (https://github.com/philres/ngmlr), but it doesn't work for genome alignment. The curent bam files from minimap2 doesn't work in Sniffles. Do you will have a SV caller for minimap 2 or do you have recommendation what SV caller should I use?
Thank you!
The text was updated successfully, but these errors were encountered: