Hello,
I want to simulate read-pairs using a read-length greater than 35 (up to 75). If I run MAQ, this works:
maq simulate -N 1000 -1 35 -2 35 out.read.1.fastq out.read.2.fastq human_b37_chr22.fasta calib-36.dat
But this does not:
maq simulate -N 1000 -1 70 -2 70 out.read.1.fastq out.read.2.fastq human_b37_chr22.fasta" calib-36.dat
calib-36.dat was downloaded from http://sourceforge.net/projects/maq/...data/20080929/. The flags "-1" and "-2" do not work if the user wants read lengths greater than those specified in the .dat file.
If anybody knows about other sequencers and other simulators, that is helpful too. I want to sequence a genome subject to random mutations in the sequencing process.
I want to simulate read-pairs using a read-length greater than 35 (up to 75). If I run MAQ, this works:
maq simulate -N 1000 -1 35 -2 35 out.read.1.fastq out.read.2.fastq human_b37_chr22.fasta calib-36.dat
But this does not:
maq simulate -N 1000 -1 70 -2 70 out.read.1.fastq out.read.2.fastq human_b37_chr22.fasta" calib-36.dat
calib-36.dat was downloaded from http://sourceforge.net/projects/maq/...data/20080929/. The flags "-1" and "-2" do not work if the user wants read lengths greater than those specified in the .dat file.
If anybody knows about other sequencers and other simulators, that is helpful too. I want to sequence a genome subject to random mutations in the sequencing process.
Comment