DWGSIM 0.1.4: whole genome NGS simulator

[nilshomer]

We are pleased to release DWGSIM version 0.1.4.

New features include:
- support for Ion Torrent data
- more metrics to focus on within dwgsim_eval
- set the random seed to generate deterministic variants
- can input a mutations.txt file or a bed-like file to specify the variants to simulate (the -m or -b option)
- new mutations.txt format to specify mutation strand, not just ploidy

This release also includes bug fixes:
- some insertions were not always left-justified
- command line option checking
- better usage

Please see the Documentation.

[TiborNagy]

Hi,

Could You show me an example how could I generate IonTorrent reads? I have try the following commands (Version: 0.1.8):

dwgsim -c 2 -B -f auto reference.fasta iontorrent/testdata
dwgsim with these settings wait forever.

dwgsim -c 2 -B -f force brca1.fasta iontorrent/masodik
dwgsim with these settings give me an error message:
[dwgsim_core] Updating error rate for end 1
[dwgsim_core] 0dwgsim: src/dwgsim.c:271: generate_errors_flows: Assertion `opt->flow_order_len != i' failed.
Aborted

Thanks

[nilshomer]

Hi,

Could You show me an example how could I generate IonTorrent reads? I have try the following commands (Version: 0.1.8):

dwgsim -c 2 -B -f auto reference.fasta iontorrent/testdata
dwgsim with these settings wait forever.

dwgsim -c 2 -B -f force brca1.fasta iontorrent/masodik
dwgsim with these settings give me an error message:
[dwgsim_core] Updating error rate for end 1
[dwgsim_core] 0dwgsim: src/dwgsim.c:271: generate_errors_flows: Assertion `opt->flow_order_len != i' failed.
Aborted

Thanks

I think you forgot to include the flow order (-f) option. I am going to add some better command line checking. Also, it looks like you are generating paired reads (use -2 0 to turn that off).

So try (assuming a TACG flow order):
dwgsim -2 0 -c 2 -B -f TACG force brca1.fasta iontorrent/masodi

[TiborNagy]

Thank You. The flow order string was the main problem.

[chintanspy]

no output for ion torrent data

I am trying to generate data for ion torrent.

I tried following command:
dwgsim -2 0 -c 2 -B -f ATGC genome/e_coli_K12_DH10B.fasta sim_SE_ion

Output:
[dwgsim_core] Updating error rate for end 1
[dwgsim_core] 1000000
[dwgsim_core] Updated with scaling factor 0.45297!
[dwgsim_core] Escherichia_coli_K-12_DH10B length: 4686137
[dwgsim_core] 1 sequences, total length: 4686137
[dwgsim_core] Currently on:
[dwgsim_core] 7046823
[dwgsim_core] Complete!

Files:
-rw-r--r-- 1 root root 0 2012-04-17 21:43 sim_SE_ion.bfast.fastq
-rw-r--r-- 1 root root 0 2012-04-17 21:43 sim_SE_ion.bwa.read1.fastq
-rw-r--r-- 1 root root 0 2012-04-17 21:41 sim_SE_ion.bwa.read2.fastq
-rw-r--r-- 1 root root 0 2012-04-17 21:43 sim_SE_ion.mutations.txt

my ouput files are empty can you explain why is it happening so?

[chintanspy]

version 1.11 vs 1.10

I forgot to mention the version of dwgsim which i used...it was 1.11..i had downloaded from git

Now, i downloaded version 1.10. from sourceforge.net, to my surprise I am getting data for ion torrent and not for illumina and SOLiD.

Can someone explain what is the problem???

[chintanspy]

question solved

Finally, got the answer...it was the 'hard disk space' issue.

But i had doubt..as in....how is illumina data diffferent from ion torrent data??
Both the outputs are in fastq and basespace

How is the flow order(for ion torrent data) implemented to generate ion torrent data?

[nilshomer]

Check out the differences in the technology to understand why Ion Torrent is different from Illumina. You can even download public datasets that each company provides. Understanding the technology will help you understand flow order and the like.

[chintanspy]

I guess i was not clear about my question.......

I have known the difference between the 2 technologies

I wanted to know how the 'flow order' information used to generate simulated data?

Sorry if I am still unclear.

[nilshomer]

I guess i was not clear about my question.......

I have known the difference between the 2 technologies

I wanted to know how the 'flow order' information used to generate simulated data?

Sorry if I am still unclear.

454 and Ion Torrent produce estimates of hompolymers, one for each flow. The errors occur when these homopolymers are misestimated. So for Ion Torrent data, dwgsim introduces errors by misestimating the hompolymer length.

[chintanspy]

454 and Ion Torrent produce estimates of hompolymers, one for each flow. The errors occur when these homopolymers are misestimated. So for Ion Torrent data, dwgsim introduces errors by misestimating the hompolymer length.

You mean to say the Carry Forward and Incomplete Extension errors are introduced, right?
it is also depended on the length of the read, right?

Also, does -f option take pattern for flow order ?

[nilshomer]

You mean to say the Carry Forward and Incomplete Extension errors are introduced, right?
it is also depended on the length of the read, right?

Also, does -f option take pattern for flow order ?

CAFIE is not explicitly modeled, no, but a more simple error model (overcall/undercall) is used. The error rate can be adjusted to increase/decrease across the read. Finally, any flow order can be used as long as it has at least one instance of each nucleotide.

[chintanspy]

Thank you for quick replies. Things are pretty much clear.

There is an option to generate strand specific reads:

-S INT generate reads [0]:
0: default (opposite strand for Illumina, same strand for SOLiD/Ion Torrent)
1: same strand (mate pair)
2: opposite strand (paired end)

The default is '0' which means it should generate reads from same strand for Ion torrent. I am getting the reads from opposite strand as well.

I tried the following command
dwgsim -1 30 -e 0 -E 0 -r 0 -X 0 -y 0 -R 0 -2 0 -c 2 -B -f ATGC test.fa ionData/test

my test.fa had following sequence
AAAATGCAAAATCTGAAAAAACGTTTTGGGAAAAAAAAAA

####Output reads
count Reads
8 AAAATCTGAAAAAACGTTTTGGGAAAAAAA
9 AAAATGCAAAATCTGAAAAAACGTTTTGGG
2 AAATCTGAAAAAACGTTTTGGGAAAAAAAA
7 AAATGCAAAATCTGAAAAAACGTTTTGGGA
8 AATCTGAAAAAACGTTTTGGGAAAAAAAAA
8 AATGCAAAATCTGAAAAAACGTTTTGGGAA
7 ATCTGAAAAAACGTTTTGGGAAAAAAAAAA
4 ATGCAAAATCTGAAAAAACGTTTTGGGAAA
6 CAAAATCTGAAAAAACGTTTTGGGAAAAAA
2 CCCAAAACGTTTTTTCAGATTTTGCATTTT
9 GCAAAATCTGAAAAAACGTTTTGGGAAAAA
7 TCCCAAAACGTTTTTTCAGATTTTGCATTT
9 TGCAAAATCTGAAAAAACGTTTTGGGAAAA
5 TTCCCAAAACGTTTTTTCAGATTTTGCATT
5 TTTCCCAAAACGTTTTTTCAGATTTTGCAT
7 TTTTCCCAAAACGTTTTTTCAGATTTTGCA
7 TTTTTCCCAAAACGTTTTTTCAGATTTTGC
3 TTTTTTCCCAAAACGTTTTTTCAGATTTTG
2 TTTTTTTCCCAAAACGTTTTTTCAGATTTT
5 TTTTTTTTCCCAAAACGTTTTTTCAGATTT
8 TTTTTTTTTCCCAAAACGTTTTTTCAGATT
5 TTTTTTTTTTCCCAAAACGTTTTTTCAGAT

can you please explain why is it happening so?

[nilshomer]

Two things:
1. Could you try the latest code here: https://github.com/nh13/dwgsim. See commit: 6c33e0e5c64c816a5b95c7eac155eef2b4b8155c
2. Can you post two reads in your FASTQ that map on opposite strands?

[chintanspy]

1. Could you try the latest code here: https://github.com/nh13/dwgsim. See commit: 6c33e0e5c64c816a5b95c7eac155eef2b4b8155c

I will try and let you know the output.

2. Can you post two reads in your FASTQ that map on opposite strands?

The following reads map on to opposite strand
CCCAAAACGTTTTTTCAGATTTTGCATTTT
TCCCAAAACGTTTTTTCAGATTTTGCATTT
TTCCCAAAACGTTTTTTCAGATTTTGCATT
TTTCCCAAAACGTTTTTTCAGATTTTGCAT
TTTTCCCAAAACGTTTTTTCAGATTTTGCA
TTTTTCCCAAAACGTTTTTTCAGATTTTGC
TTTTTTCCCAAAACGTTTTTTCAGATTTTG
TTTTTTTCCCAAAACGTTTTTTCAGATTTT
TTTTTTTTCCCAAAACGTTTTTTCAGATTT
TTTTTTTTTCCCAAAACGTTTTTTCAGATT
TTTTTTTTTTCCCAAAACGTTTTTTCAGAT