This is a test version of Biostars. For the public version, visit https://www.biostars.org.
obtain a simulated bam file for variants on existed bam file

Hi all,

I constructed the variant calling pipeline for BAM file. However, when it was applied to positive samples, I did not find any well-known variants in the VCF file. What I am considering is that there might be less-sensitive method in my pipeline.

So I want to create a simulated BAM file based on the existed bam, just change the reads in the bam to create a variant, etc., rs6071.

Is there any convenient tool to do such, and also can I obtain a depth-aimed bam? For example, depth for alleles, A:G, 500:500, OR A:G, 800:200.

Thanks

simulated variant bam
• 4,670 views • link
updated by Ram Ram • written by J.F.Jiang J.F.Jiang

2 answers

BAMSurgeon seems to be able to do what you want (and much more).

• 0 views • link
written by h.mon h.mon

Exactly what I am looking for though bugs are still in this program needing debug

• 0 views • link
written by J.F.Jiang J.F.Jiang

extract the region of interest as fasta, add manually the mutation, run wgsim to create a set of fastq and map those fastqs using your pipeline.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks, wgsim can works fine to simulate a set of fastq with mutations, however, the program will randomly select mutation site in the given intervals. What I want to is to mutate the specific site, and can change the alleles depth seperately.

• 0 views • link
written by J.F.Jiang J.F.Jiang

Hi,

Where you able to perform this task. if so can you please let us know

• 0 views • link
written by bassanio • bassanio

Are you aware you are resurrecting a four-years old thread?

Did you read the manual and tried to use BAMSurgeon? In this case, did you run into any problems in running BAMSurgeon?

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this