Hi everyone,
is there a way to introduce a specific artificial deletion in a bam file? I'm searching something that take a particular interval of coordinates and remove it from the bam file e.g. chrx:3122135-3425222
•
written
by
enee •
2 answers
Forgive me if I wasn't very clear, I mean that I would need to remove the reads that map a certain chromosomal interval from the bam file.
samtools view -L deletions.bed -O BAM --unoutput output.bam in.bam > /dev/null
e I have to find a way that partially removes the reads and not all of them
split your bam with
samtools view --subsample 0.5 -O BAM --unoutput output1.bam -o output2.bam in.bam
remove reads with deletion from output2.bam
and then samtools merge with output1.bam and output2.with_del.bam
If you just want to remove the reads, it's relatively easy.
If you want to simulate a truncation, this is difficult, because you will need to modify instead of remove the reads so they can cross both breakpoints of a truncation.
Log in to answer this question.
-
Jellyfish docker image error with bam files
written by enee •
-
Introduce SNPs on FASTA
written by Riccardo •
-
Fine-tuning NCBI taxid / taxidlist
written by theclubstyle •
-
Extracting length of deletions from a specific region from a BAM file
written by a.beggs •
-
Remove reads from bam file that partially covers a region
written by arbarharw •
-
How to introduce artificial mutation in bam
written by bassanio •
-
Identify deletion in multiple bam files compared to the same bacterial reference
written by David
-
Tool to calculate extreme most positions in a bed file for a given window
written by Aishwarya Kulkarni
-
correct notation of indels/duplication in .vcf file
written by User6891
-
Intersect Gene Annotation With Specific Position Or Genomic Interval
written by Nicolas Rosewick
A bam file is a representation of (typically) sequencing reads mapped against a reference genome. Is is not clear what you want. Do you want for the reference genome to have the deletion? Or do you want to remove the reads mapping to that region?
Forgive me if I wasn't very clear, I mean that I would need to remove the reads that map a certain chromosomal interval from the bam file. I need these samples with the artificial deletions to test a tool I'm developing.