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

How to deal with the same reads’ name? #19

Closed
wzhangahas opened this issue Jun 17, 2021 · 3 comments
Closed

How to deal with the same reads’ name? #19

wzhangahas opened this issue Jun 17, 2021 · 3 comments

Comments

@wzhangahas
Copy link

When I check the different chromosomes reads’ names, there are some reads’ names belonged to different chromosomes, should we remove the same reads’ names in different chromosomes, or just keep them in different chromosomes? The size of chromosome 1 of Arabidopsis is 29.9M. When I assembled the fq of chromosome 1 of Arabidopsis with the same reads’ names by canu 2.0, the size of contig is 36.2M, and when I assembled the fq of chromosome 1 of Arabidopsis without the same reads’ names by canu 2.0, the size of contig is 19.6M.

@mawad89
Copy link
Collaborator

mawad89 commented Jun 17, 2021

You don't need to delete anything. Use the fq file as it is. You can try delete the reads with low mapping quality (Ex: Q<10) but my suggestion is to keep it.
If the Chr1 size in Arabidopsis is 36.2M. usually you will find some contigs belongs to another chromosomes or duplicated contigs, ignore them.

@wzhangahas
Copy link
Author

wzhangahas commented Jun 18, 2021 via email

@mawad89
Copy link
Collaborator

mawad89 commented Jun 18, 2021

Hi Wei
If you mean the assembly you can find it at https://doi.org/10.5281/zenodo.3840274. if you want something else let me know.

All the best;
Mohamed

@mawad89 mawad89 closed this as completed Aug 26, 2021
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

2 participants