How to create alternate reference for SNPs?
0
0
Entering edit mode
18 months ago
Parham ★ 1.6k

Hi,

I am looking for ways to create alternate references for metagenomes. I would like to introduce SNPs every other x base pairs in the new reference using the original reference. Ideally if there is a package that does this, I appreciate to know about it, otherwise what is your recommendation to do it?

snps reference • 353 views
ADD COMMENT
1
Entering edit mode

If you have the reference then mutate.sh from BBMap suite can introduce mutations that you control. It will also make a corresponding VCF. If you would like to simulate the meragenome itself first then randomreads.sh from BBMap can do that with metagenome=t flag.

ADD REPLY
0
Entering edit mode

Thanks! Yes, I have the reference and the idea is to make an alternate version of that reference with introducing mutations. If you have experience with mutate.sh I wonder if Substitution rate, 0 to 1 means 0.1 will replace 10 percent of bases?

ADD REPLY

Login before adding your answer.

Traffic: 861 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6