diff --git a/README.md b/README.md index cbb7789..6d345d3 100755 --- a/README.md +++ b/README.md @@ -79,7 +79,7 @@ By default, DNA PromethION reads (R9.4.1) will be simulated. Specify the `-x STR - `dna-r10-min`: genomic DNA on MinION R10.4.1 flowcells - `dna-r10-prom`: genomic DNA on PromethION R10.4.1 flowcells - `rna004-min`: direct RNA on MinION RNA004 flowcells -- `rna004-prom`: direct RNA on MinION RNA004 flowcells +- `rna004-prom`: direct RNA on promethION RNA004 flowcells If a genomic DNA profile is selected, the input reference must be the **reference genome in *FASTA* format**. *squigulator* will randomly sample the genome from a uniform distribution and generate reads whose lengths are from a gamma distribution (based on `-r`). If a direct RNA profile is selected, the input reference must be the **transcriptome in *FASTA* format**. For RNA, *squigulator* will randomly pick transcripts from a uniform distribution and the whole transcript length is simulated. diff --git a/src/model.c b/src/model.c index 17c14d3..9d6f6e2 100644 --- a/src/model.c +++ b/src/model.c @@ -172,7 +172,7 @@ uint32_t set_model(model_t* model, uint32_t model_id) { } else if(model_id==MODEL_ID_RNA_RNA004_NUCLEOTIDE){ kmer_size=9; num_kmer=262144; - inbuilt_model=rna004_70bps_u_to_t_rna_9mer_template_model_builtin_data; + inbuilt_model=rna004_130bps_u_to_t_rna_9mer_template_model_builtin_data; assert(num_kmer == (uint32_t)(1 << 2*kmer_size)); //num_kmer should be 4^kmer_size } else{ assert(0); diff --git a/src/model.h b/src/model.h index 8a9eb7d..dd5e8c7 100755 --- a/src/model.h +++ b/src/model.h @@ -267284,7 +267284,7 @@ static float r10_4_nucleotide_9mer_template_model_builtin_data[] = { 111.8712688 , 3 // TTTTTTTTT }; -static float rna004_70bps_u_to_t_rna_9mer_template_model_builtin_data[] = { +static float rna004_130bps_u_to_t_rna_9mer_template_model_builtin_data[] = { 100.9501553 , 3 , // AAAAAAAAA 100.8738303 , 3 , // AAAAAAAAC 100.4230422 , 3 , // AAAAAAAAG