Hello!
Could anyone please suggest a tool which would simulate short reads with only mutations supplied by the user. I am trying to evaluate variant calling methods, and have a list of variants in specific gene locations and want to find out which methods would be best for finding the variants. So, basically I am looking for a tool which would take in reference fasta and vcf file and generate Fastq reads, but in addition to that I have 2 conditions - a) There should not any mutations apart from the intended location, b) the intended mutation should be in the intended location only in some of the reads and not all the reads (ie, if the read coverage in a particular region is 60x, the mutations should be only in say 20 or 30 reads and not in all the 60 reads)
I looked through some of the previous Biostar posts for read simulation, but couldn't exactly find the best tool.
Any suggestion would be greatly appreciated.
Thanks, Raman
•
written
by
Raman2 •
1 answer
convert a VCF to bam : vcf to bam
and then convert the bam to fastq using samtools fastq
mix this ouput with the output of htslib/wgsim
Log in to answer this question.
-
Creating a personalized genome assembly for read simulation from illumina WGS data
written by scsc185 •
-
Simulate reads for positions in a vcf file
written by berndmann •
-
Generate with specific converage
written by Tim •
-
Insert mutations in a reference genome from a vcf
written by IP
-
Are there any tools to simulate mutations?
written by MAPK
-
Example RNA-seq datasets that require specific normalizations
written by GreenDiamond
-
Simulation Tools To Generate Next-Generation Sequencing Reads along with their VCF
written by manojkumarbioinfo
-
Simulated microarray expression data
written by Jan-Niklas •
-
The problem with sensitivity analysis of somatic variant caller
written by phongphak.06 •
-
Tool: Art: Simulation Tools To Generate Synthetic Next-Generation Sequencing Reads.
written by Istvan Albert
