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
winnowmap PacBio CCS alignment #6
Comments
Although Winnowmap was started using minimap2 as base implementation, over the time I have made significant changes to the algorithm (particularly seeding and anchor chaining). It may be that
I'd be really curious how |
I will report back! |
I can report that winnowmap map_pb parameter is slightly better than minimap2 asm20 parameter for the somatic mutations I am attempting to detect! |
I think pbmm2 from PacBio uses something like this for the minimap2 preset for CCS reads
|
To whom it may concern,
I am writing to inquire whether the authors of the program had compared the results from winnowmap -ax map-pb and -ax asm20 bam file.
Minimap2 recommends using asm20 for mapping CCS reads and when I compared germline mutation results from asm5, asm10 and asm20 bam files, asm20 parameter performed the best for germline mutation calling. I was also wondering if there was reason for recommending the alternative parameter (map-pb) instead of the standard parameter (asm20)
winnowmap recommendation
winnowmap -W repetitive_k15.txt -ax map-pb ref.fa hifi.fq.gz > output.sam
minimap2 recommendation
./minimap2 -ax asm20 ref.fa pacbio-ccs.fq.gz > aln.sam # PacBio CCS genomic reads
Many thanks,
Regards,
Sangjin
The text was updated successfully, but these errors were encountered: