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

re-root to wh1? #13

Open
roblanf opened this issue Aug 14, 2020 · 4 comments
Open

re-root to wh1? #13

roblanf opened this issue Aug 14, 2020 · 4 comments

Comments

@roblanf
Copy link
Owner

roblanf commented Aug 14, 2020

One suggestion from Russ CD is to re-root on WH1, since i use that as a reference for the alignments. This would make mapping mutations more sensible looking.

@AngieHinrichs
Copy link

Actually, Wuhan-Hu-1/EPI_ISL_402125/MN908947/NC_045512 (29903 bases) and WIV04/EPI_ISL_402124/MN996528 (29891 bases) are identical except that Wuhan-Hu-1's polyA at the end is 12 bases longer than WIV04's. There is no WIV04 base that differs from Wuhan-Hu-1 (thank goodness). Since the beginning and end of the sequence alignment are masked, there should be no difference between the two sequences as far as phylogeny is concerned, right? So I think this is moot.

@AngieHinrichs
Copy link

Sorry, I assumed you were using WIV04|EPI_ISL_402124 because that is GISAID's preferred reference (https://www.gisaid.org/epiflu-applications/hcov-19-reference-sequence/). But I see in global_tree_gisaid.sh it's actually hCoV-19/Wuhan/WH04/2020|EPI_ISL_406801|2020-01-05. That is definitely different (relative to WIV04 or NC_045512 it has 8782 C>T and 28144 T>C).

@AngieHinrichs
Copy link

We've been rerooting to Wuhan-Hu-1/NC_045512 for our purposes, so this doesn't directly affect us. I was hoping to give you a reason to simply close this issue, but I was just mixing up WIV04 and WH04, oops.

@roblanf
Copy link
Owner Author

roblanf commented Nov 17, 2020 via email

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