Hello,
I have written a CNV Simulator, which takes fasta(eg. chr22 from hg18, repeat masked, all 'N' removed) file as input and inserts CNV at random locations.
I then generate reads using maq simulate.
To achieve a coverage of around 20x.
I give,
maq simulate -d 200 -N 2466490 -1 75 -2 75 sim_read1.out sim_read2.out chr22_sim_ref_fi.fa simupar.dat
After generation of reads and Alignment by BFAST, to verify I wrote a perl script to bin the sam file into 5kb windows and calculate how many reads map to each window.
I find that most(>95%) of the read are mapping to only the last window.
I read the man page, still not very clear about few issues.
If I specify the parameters for simulate, is it necessary to provide the .dat file?
If yes, then what kind of fastq file should I train simutrain on?
Also the man page mentions that there is a .dat file available for download, but I am not able to find it.
I have also tried with dwgsim now, the same problem occurs, all reads mapping to the last window.
Can anyone point out what I am doing wrong.
Thanking you,
Prakhar Gaur
I have written a CNV Simulator, which takes fasta(eg. chr22 from hg18, repeat masked, all 'N' removed) file as input and inserts CNV at random locations.
I then generate reads using maq simulate.
To achieve a coverage of around 20x.
I give,
maq simulate -d 200 -N 2466490 -1 75 -2 75 sim_read1.out sim_read2.out chr22_sim_ref_fi.fa simupar.dat
After generation of reads and Alignment by BFAST, to verify I wrote a perl script to bin the sam file into 5kb windows and calculate how many reads map to each window.
I find that most(>95%) of the read are mapping to only the last window.
I read the man page, still not very clear about few issues.
If I specify the parameters for simulate, is it necessary to provide the .dat file?
If yes, then what kind of fastq file should I train simutrain on?
Also the man page mentions that there is a .dat file available for download, but I am not able to find it.
I have also tried with dwgsim now, the same problem occurs, all reads mapping to the last window.
Can anyone point out what I am doing wrong.
Thanking you,
Prakhar Gaur
Comment