Comparative Snp Analysis
Hello, I am trying to compare the degree of A-to-G editing in a near-to-isogenic pair of cell lines. I have two biological replicates and have mapped with Bowtie and BWA, followed by a samtools mpileup...
View ArticleWhat Is The Fastest Method To Determine The Number Of Positions In A Bam File...
I have two very large BAM files (high depth, human, whole genome). I have a seemingly simple question. I want to know how many positions in each are covered by at least N reads (say 20). For now I am...
View ArticleBedtools: Top N Most Similar Regions When Comparing Two Bed/Wig/Bam Files?
Is there an easy way of finding, probably with bedtools, given a window size, the top N most correlated regions when comparing two bed/wig files? For example, in comparing two bed/wig/bam files that...
View ArticleGeneral Considerations For Genomic Overlaps?
Hello I was wondering about general considerations for performing overlap of genomic regions and doing Monte Carlo-type statistics. Below I have made a description of how I do it, unfortunately Im not...
View ArticlePicking Random Genomic Positions
I do have a set of TF binding coordinates and want to see if there is any significant overlap with an open chromatin annotation.Example of TF coord: chr1 19280 19298 chr1 245920 245938 chr2 97290 97308...
View ArticleReporting The Bam Reads Overlapping A Set Of Intervals With Bedtools
I am trying to use bedtools to pull out the reads falling directly within a set of BED coordinates. While this command does it successfully:intersectBed -abam mybam.bam -b intervals.gff -wa -wb -f 1 |...
View ArticleHelp With Exception When Using Bedtools Coveragebed With Paired Alignment....
I use bwa mem to align paired reads to few hundreds of microbial contigs; then I sort the alignment, and trying to get a coverage using bedtools genomecov -ibam alignments.paired.sorted.bam -bg...
View ArticleIntersectbed: Return Reads In Fraction In Input Files
I have a question with respect to intersectBED and multiple input files:Is it possible to return reads which are present in, say 8/10 input files, without fractioning the reads in smaller...
View ArticleClik here to view.
Extract Only Paired-End Reads That Map A Specific Interval
Hi,Is it possible to extract paired-end reads that map to a specific interval ( from a bam file ). I tried with intersectBed :intersectBed -abam align.bam -b interval.gff3 -wa > result.bamhere's the...
View ArticleConverting Gff To Bed With Bedtools?
I use bedtools's sortBed utility to sort BED files for various operations. It takes as input GFF files as well. However, when I feed it a GFF file as in:sortBed -i myfile.gffit outputs it as GFF, not...
View ArticleHeatmap Of Read Coverage Around Tsss
I am trying to plot a heatmap of read density around a feature of interest (TSSs) very common in genomics papers. something like this (B): However, I am struggling a bit in getting to look "right". A...
View ArticleHow Can I Merge Intervals ?
Hello everybody, I should be grateful if you would kindly help me de fix my problem. I have a table like that :Chromosome start end info1 info2 chr01 1 100 15 35 chr01 150 300 15 39 chr01 299 750 16...
View ArticleQuestion about number of reads within intervals
Hi there, This question is very basic but I need to ensure that I'm going on the right way. I need to calculate the number of reads falling inside my bed intervals and the number of reads falling...
View ArticleConvert Bamtobed Score
Hey,just a short question....is there a possibility to set the score in the bed file to "1" an not to the the alignment score?? arguments -tag and -ed only use BAM alignment tags... ?!? :/Cheers!
View ArticleWhich Of The Genes Are Enriched With Repeat Elements
I would like to know which of my genes are enriched with repeats of LINE/SINE/ERV etc. elements. I have a bam file and the repeats in bed format. As far as I know BAM files contains aligned data for...
View ArticleCounting Features In A Bed File
I have a file in the following BED formatChr1 1022071 1022105 + Chr1 1022071 1022105 + Chr1 1022072 1022106 - Chr1 1022072 1022106 - Chr1 1022072 1022106 - Chr1 1022072 1022106 -I am trying get the...
View ArticleSplice Junction file intersection with genome annotation
Hello, I have a tab delimited format Splice Junction file and the file looks something like this: chr1 11212 12009 1 1 0 0 2 48 chr1 11672 12009 1 1 0 0 1...
View ArticleIntersectbed Tool Generating Empty File
I have used the Bedtools command intersectBed to check the overlap between two bed files. A is my INDEL file and B is my Reference file. But it is producing an empty output file. I thought the problem...
View Articleproblem with bedtools complement: how to extract strand information?
Hi everyone,I have a minor trouble when using bedtools complement. My gff3 file contains Chr start end strand ect. information. I want to use bedtools complement to extract other regions, however,...
View ArticleHow To Check Whole Genome With Bigwigsummary ?
Hi,I have question about bigwigsummary tools ,I have my start and end positions and my bigwig file but I want to check whole genome instead of chromosome by chromosome Is there any option to use this...
View Article