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.
•
written
by
prernagoel04 •
1 answer
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.
Log in to answer this question.
-
Specific case of "bedtools intersect"
written by predeus
-
Interpreting Bedtools Intersection Between BAM and BED Files
written by Steven Mick •
-
Bedtools intersect 0- and 1-based coordinate files
written by maxrwjones •
-
Find count of reads in BAM file that fall within BED region
written by KH •
-
Bash script problem
written by valastesolutions •
-
How to calculate allele counts(Ref and alt) using a target bed file and aligned bam sample files
written by prernagoel04 •
-
Why did Bedtools intersect output a zero byte file
written by jaqx008
-
How to apply bedtools intersect function with GNU parallel for bacth BED file comparison?
written by texaswhiskey1969 •
-
Does multiIntersectBed works with -u and -f parameter?
written by ivivek_ngs
-
bedtools intersect - something wrong with chromosome numbers >= 10?
written by flyamer •
