0
$\begingroup$

I want to benchmark my DNA sequencing pipeline. In order to do that I want some gold-standard files in every step(.vcf,.bam,.fastq). I want to generate/simulate a bam file of reads for a given set of vcf. How to do that?

The GIAB project has gold standard vcf, bam, and fastq. But they are not perfect for pipeline testing. Their size is quite big. Also, they are real data. My question was how to generate simulated Bam from vcf.

$\endgroup$
4
  • 1
    $\begingroup$ Have you tried making a modified fasta containing the variants and then simulated reads using that. Normally then you'd mix those together with reads simulated against the reference sequence at different ratios. $\endgroup$
    – Devon Ryan
    Nov 28, 2023 at 14:51
  • $\begingroup$ Actually, I am planning to do that now. The thing is, I couldn't determine which one is the most convenient way. I am planning to do that using Python. Can you write it in the answer and suggest a convenient way to do that? $\endgroup$ Nov 28, 2023 at 15:20
  • 1
    $\begingroup$ Hint, you want to make a fasta from a VCF (e.g., biostars.org/p/360900) and then use any random read simulator. No need to program anything, that'd be a waste of time. $\endgroup$
    – Devon Ryan
    Nov 28, 2023 at 21:47
  • $\begingroup$ Sorry, I didn't read your comment earlier properly. Actually, I used many simulators like NEAT, ART, MASON, etc. But, unfortunately, they are too much buggy. If you go to their issues, you will see their problems. That's why I wanted to write my own code in the first place. However, among the simulators, NEAT probably is the least buggy one. $\endgroup$ Nov 29, 2023 at 18:23

0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.