Question: Matching mutation in bed files to their corresponding bam files
gravatar for prernagoel04
20 days ago by
prernagoel040 wrote:

Hi all. I am doing an RNA Seq analysis for which I have created a few BAM files after alignment of the sequences to reference genome. I already have a .bed file which has some specific mutations from the same data(i.e a single file having mutations from all samples). I want to match those mutations to the particular sample (BAM) files to which they belong since the name assigned to mutations in bed files are different. I am trying to use bedtools intersect for this but it doesn't seem to be working properly. Could you suggest me any other way to do so or how to do it using bedtools itself. Thank you.

rna-seq bedtools • 76 views
ADD COMMENTlink modified 18 days ago by bari.ballew120 • written 20 days ago by prernagoel040
gravatar for bari.ballew
18 days ago by
bari.ballew120 wrote:

BAM files contain reads aligned to a reference sequence. Bed files contain regions defined by chromosome, start, and stop positions. Neither of these files contain called mutations. It sounds like you may need to use a variant caller to first call variants in your BAM files. Then, you can intersect the resulting VCF with your bed file to find mutations in the regions of interest.

ADD COMMENTlink written 18 days ago by bari.ballew120
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1964 users visited in the last hour