bedtools intersect - something wrong with chromosome numbers >= 10?
Hi!I have an alignment (.bam) of reads to mm9 genome. I sorted it with samtools sort, so that later I can use -sorted key with bedtools. I also created a .bed-file with regions of interest, in which I...
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 ArticleExtracting Genomic Coverage Information Across Different Samples
Hello, I have 3 bam files that i wanted to compare against each other. For example i have reference file with 10,000 sequences. I have paired end reads sequenced for 3 different samples. 1) Sample 1 is...
View ArticleBedtools on Cygwin problem.
Hi  I'm trying to install the latest release of Bedtools via Cygwin but there's a weird error during process. I know this isn't the best solution, but I do not have an other choice. Perhaps anyone...
View ArticleClik here to view.
How To Explain Uneven Coverage Of A Dna Seqment Obtained Via Pcr Amplification.
Experiment: deep sequencing for mutants in 700nt fragment.the fragment of dna was preamplified by primers flanking the fragment followed by hiseq.per base coverage was calculated by coverageBed -d...
View ArticleRaw Counts From Cufflinks Output
Hi, I want to ask how to get the raw counts from the output of cufflinks. One way to do this is to use the fpkm.raw counts = FPKM * (length of that transcript/1000) * (# of mapped reads / 1e6)The FPKM...
View ArticleBedtools Compare Multiple Bed Files?
I've been dealing with comparison between two bed files using intersectBed -a -b command. I'm just wondering, is there any commands in Bedtools which can help us compare multiple bed files? Say, I have...
View ArticleFastafrombed Problem
hi, I try this tools from BedTools but it doesnt work!$ cat testgenome404.fa>chr1 AAAAAAAACCCCCCCCCCCCCGCTACTGGGGGGGGGGGGGGGGGG $ cat test.bed chr1 5 10 $ ./fastaFromBed -fi testgenome404.fa -bed...
View ArticleConverting Sam Files To Bam Files - Reproduce Results Nature Paper:...
I want to reproduce the results that people achieved in the following Nature paper: Transcriptome genetics using second generation sequencing in a Caucasian...
View ArticleTool: Bedtools: Analyzing Genomic Features
All practicing bioinformaticians will face problems that require them to compare, query and select genomic features across an entire genome. As it happens efficient interval representation and query is...
View ArticleWhat Is The Best Way To Run Bedtools In Parallel With Blocking
Say I am working on a server with a shared file system and 4 quad core nodes (I/O is not an issue, 16 cores total). I want to run coverageBed across 20 files. Currently I have a shell script that would...
View ArticleDoes Bedtools Intersect -V Consider Unmapped Reads "As Not In B"
bedtools intersect -v -abam my.bam -b myregions.gff > notinmyregions.bamwould we see reads with 4 in the FLAG field - i.e. unmapped reads in notinmyregions.bam
View Articlebedtools: extracting no coverage regions
Hello,I am not sure if this has been answered before as I looked and couldn't find a simple answer.I have a bam file, and all I want is to annotated all regions with 0 coverage in bed format. Is that...
View ArticleGetting Rna Sequences From Gff And Fa Files
Hi. I have a folder full of .fa files, and a .gff. The gff file contains information about which loci look like they code for RNA sequences. The .fa contain the DNA sequences for a set of human...
View ArticleBedtools Multicov Need A Bam Index File Specification Option
bedtools version 2.16.2 multicov used to compute the multiple sample coverage given a feature file(gtf bed).format: bedtools multicov -bams alin1.bam aln2.bam .. -bed capturRegion.bed...
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 ArticleMerging/Intersecting Different Gene Annotations - Should I Extend Coordinates?
I want to create gene data-set (as big as possible), hence I am using several gene annotations. However, genes in different annotations overlap (it's the same gene). For reducing biases I overlap...
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 ArticleRemove Intronic Regions in .BAM
HiI have a .BAM file which contains discordantly and concordantly mapped mate-pairs. I used bedtools Pairtobed to extract the mate-pairs which both show overlap with targeted regions (Illumina target...
View ArticleTo Group Items In Bed Files
For example, we now have a bed file:chr1 23455 45678 chr1 23446 45663 chr1 23449 45669 chr1 30000 31000Is there anyway to group the first three lines, while leaving the last line alone? I know Bedtools...
View Article