I want to use dwgsim to generate CNVs and InDels.
Searched for documentation, didn't find anything useful except "http://sourceforge.net/apps/mediawiki/dnaa/index.php?title=Whole_Genome_Simulation".
Tried using it myself, its working, am getting an output, hence have some questions!
1) "chr1 763 A C -" what does this line of output mean?
As far as I can understand the program is taking Nucleotides from Ref Chr1, and changing them based on some parameter and then generating reads.
In short can you please explain the algorithm behind this??
In the output format, what is meant by "start end 1" & "start end 2" ?
2) Why does the program produce 3 output files(2 for bwa, 1 for BFast)?
Any help will be highly appreciated!!!
Searched for documentation, didn't find anything useful except "http://sourceforge.net/apps/mediawiki/dnaa/index.php?title=Whole_Genome_Simulation".
Tried using it myself, its working, am getting an output, hence have some questions!
1) "chr1 763 A C -" what does this line of output mean?
As far as I can understand the program is taking Nucleotides from Ref Chr1, and changing them based on some parameter and then generating reads.
In short can you please explain the algorithm behind this??
In the output format, what is meant by "start end 1" & "start end 2" ?
2) Why does the program produce 3 output files(2 for bwa, 1 for BFast)?
Any help will be highly appreciated!!!
Comment