Problems Extracting Non-Snps From A Vcf File
Hello,In an SNP analysis, I am trying to extract those editing sites no found in the dbSNPs vcf file I have downloaded a couple of files (All SNPs and Common/Medical SNPs) from...
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 ArticleRandom shuffling of features leaving gene models intact
I am looking for a tool that can randomly shuffle gff features into intergenic regions, but leaving the gene-models 'intact', so that at least all features of a gene are placed on the same contig and...
View ArticleHow Can I Compare And Merge Bed Files
I have three bed files with chrNo, start, end position and type. I need to compare each chrNo, start and end position of one file with 2 other files and write the common one in a new file. Can any one...
View ArticleSamtools or Bedtools: How to filter a bam file with a bed file using strand...
HiI would like to filter a bam file, keeping only reads overlapping with genomic intervals from a bed file. I used samtools for this:samtools view -b -h -L bedfile.bed bamfile.bam However the -L option...
View ArticleHow To Get Fasta Format Using Fastafrombed Or How To Turn Linearized Fasta To...
I extracted sequences with fastaFromBed and have no complains about the BEDTools which is really awesome thing.Otherwise extracted sequences look like this:>chr19:13985513-13985622...
View ArticleHow To Count Genes In Genomic Regions Using A Gtf/Gff3 And A Bed File Of Regions
I'd like to count the number of unique genes in a gff file falling within a list of genomic regions. With bedtools I can count the number of regions within the gff which is almost what I want, but not...
View ArticleDoes Bedops Have A Command Similar To The Bedtools Makewindows?
With bedtools you can make genomic windows from a genome file or a bed fileinput.bedchr1 1000000 1500000 chr3 500000 900000[prompt]$ bedtools makewindows -b input.bed -w 250000chr1 1000000 1250000 chr1...
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 ArticleIntersectbed Overlap
Hi,I've a question about intersectBed. Is it possible to extract only alignment like this :chromosome =============================================================== BED/BAM A ==============...
View ArticleBedtools Genomecoveragebed Usage : How To Create A Genome File?
I am using BEDTOOLS and the following command to get the coverage file:$ ./genomeCoverageBed -ibam ~/GG_project/trim/ecoli.bam -g > ~/GG_project/trim/coveragewhere ecoli.bam is my sorted bam file,...
View ArticleBed File Of Mapq Sliding Window On A Bam File?
There may already be a recipe for this, so asking first before reinventing the wheel: I would like to create a bed file where the score is the average mapQ from the reads of the input.bam file. I think...
View ArticleCalculate reciprocal overlap for thousands of samples
I have around 20k samples with BED files. How can I calculate reciprocal overlap for each segment? I want to find all segments with 50% reciprocal overlap or better.
View ArticleBedtools subtract not dealing well with large datasets
I am using bedtools subtract with large datasets and it keeps crashing, giving the following error terminate called after throwing an instance of 'std::bad_alloc'Is there a way to get over this...
View ArticleSimple Redirection, I/O Problem With Bedtools
Hi Guys, Just a quick question. Its more of a Bash question rather than Bioinformatics, with Bedtools in question.I mostly pipe the bedtools I/O. Here's a general scenario :sed 1d fileA.bed |...
View ArticleHow To Find The Nearest Gene To A Retrotransposon Insert?
Hi,I have a BED file with the position of retrotransposons in the mouse genome and I would like to find the nearest gene, the distance to that gene and whether it is on the + or - strand. There are so...
View ArticleGenomecoveragebed - Bedtool For Reporting Per Base Genome Coverage
Hi Everyone I would nedd some help on genomeCoverageBed tool. This tools when used for finding per base genome coverage uses an option -d. I am actually interested in finding read counts for each base...
View Articlemacs and bedtools
HelloI have MACS2 output and now looking for peaks which are situated in introns. I have bed file with introns from USCS for my species. What file with peaks should I use for bedtools intersection?...
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 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 Article