Question: wgsim outputs different files
0
gravatar for signauncreal
3.9 years ago by
signauncreal0 wrote:

Hey guys,

im simulating SVs on the human genome using svsim. The output files im using with vgsim to simulate sequencing reads. I have observed something i cannot make sense of:

depending on whether i use svsim in contig or whole genome mode, the result i get using wgsim on either files differs from each other. Shouldn't both output files from wgsim be completely identical?

 

Can you guys help me out?

 

sequencing • 1.3k views
ADD COMMENTlink modified 2.9 years ago by Biostar ♦♦ 20 • written 3.9 years ago by signauncreal0
0
gravatar for iraun
3.9 years ago by
iraun3.6k
Norway
iraun3.6k wrote:

From the SVsim manual:
"The default mode is called contig mode. In this mode, the fasta file will contain only small regions of the genome surrounding the target location or breakpoints of the SVs. The other major moade is called Whole Genome Mode (WGM). In this mode, the entire mutated genome is output to the fasta file."

So, according to the above paragraph, the fasta output file that SVsim generates depends on the mode. And as the fasta file changes, the reads that you are generating using wgsim will change too.

ADD COMMENTlink written 3.9 years ago by iraun3.6k
0
gravatar for signauncreal
3.9 years ago by
signauncreal0 wrote:

that is true, but the files give out the same information, just a marginal portion (roughly 5%) of the information given in the file originating from svSimContig is missing in the file given from svSimWG. That is not expected given the information you just stated. Anyways, what im trying to emphasize is that both output files, given they have the same format, should not differ from each other (if not stated otherwise in the documentation of wgSim).

ADD COMMENTlink modified 3.9 years ago • written 3.9 years ago by signauncreal0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1996 users visited in the last hour