Biostar Beta. Not for public use.
Addition of Y chromosome to the ref sequence
0
Entering edit mode
16 months ago
micro32uvas • 10
Pakistan

Hello,

I have aligned my WGS samples to a the latest bovine reference(B.tau5.0), available at NCBI. It is just yesterday, that in a pipeline used by 1000bull genome project, they have used the same reference merging the Y chromosome in it (B.taur5.0Y).

Now, I am not sure about it. Should I align all of my samples again? It would take huge time and resources. Is there any way to compensate and add the Y chromsome alignment that could save the current bams?

Any help is greatly appreciated.

I used bwa-0.7.17 and bwa mem for the alignment.

Ainnie

ADD COMMENTlink
0
Entering edit mode
12 months ago
France/Nantes/Institut du Thorax - INSE…

Should I align all of my samples again?

yes, adding new contigs removes false positives (think of a read matching 100% with 'Y' -> no variant, but matching with mismatches on chromosome '1' -> variant)

ADD COMMENTlink
0
Entering edit mode

What I dont get, is that the providers of the ref seq, (NCBI and Bovine database) have not given any newer version for the said purpose (The latest came by April 2018, and its done by PacBio). Now one of the pipelines, shared by the 1000 Bull genome Project, includes the Ref Seq (NCBI's Ref seq+ Y) and its provided at their site.

If the main providers have not updated a new/ modified reference, should I be even following the 1000Bull genome project's pipeline? Its more of a nightmare to start over, but is it worth the confidence?

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1