[go: nahoru, domu]

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

Align short reads to an uncorrected miniasm assembly #90

Closed
sjackman opened this issue Jan 5, 2018 · 8 comments
Closed

Align short reads to an uncorrected miniasm assembly #90

sjackman opened this issue Jan 5, 2018 · 8 comments
Labels

Comments

@sjackman
Copy link
sjackman commented Jan 5, 2018

Hi, Heng. Do you have a suggestion for parameters to align Illumina reads to an uncorrected miniasm assembly with minimap2?

@lh3 lh3 added the question label Jan 5, 2018
@lh3
Copy link
Owner
lh3 commented Jan 5, 2018

I don't know. It will be hard. Minimap2 won't work well at 15% error rate on ~150bp reads. Probably stampy or bwa-mem will work better.

@sjackman
Copy link
Author
sjackman commented Jan 5, 2018

Thanks for the info, Heng. I'll try bwa mem -xont2d. I'll also try first correcting with Racon.

@sjackman
Copy link
Author
sjackman commented Jan 5, 2018

To map nanopore reads to an uncorrected miniasm assembly for Racon, do you suggest minimap2 -xmap-ont?

@lh3
Copy link
Owner
lh3 commented Jan 5, 2018

Perhaps also add -xmap-ont -w5 to increase sensitivity.

@sjackman
Copy link
Author
sjackman commented Jan 5, 2018

I'll try that. Thanks, Heng.

@lh3
Copy link
Owner
lh3 commented Jan 11, 2018

Any luck with either bwa-mem or minimap2?

@sjackman
Copy link
Author

Rather than align to the uncorrected Miniasm assembly, I used two rounds of Racon correction with the Nanopore reads alone, and then aligned the short reads to the assembly using BWA-MEM. If I try aligning to the uncorrected assembly, I'll post back here. Thanks for your help, Heng.

@jdmontenegro
Copy link

Hi guys,

I wonder if this was ever sorted. what parameters were used to map short reads to a one-round polished miniasm assembly? Currently racon developper recommends using minimap2 with -x sr parameter, but I can see here that bwa mem should work better, which one would you finally recommend and what parameters should be fine-tuned?
Cheers,

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

No branches or pull requests

3 participants