Extracting 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 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 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 ArticleGTF2/GFF3 "feature" types and expression analysis
Hi, I aligned a few samples using STAR to the genome provided in the Illumina iGenomes UCSC hg19 bundle (here) -- I used the provided gene feature (gtf2) file as is.  Now, my motive is to calculate the...
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 ArticleGetting The Average Coverage From The Coverage Counts At Each Depth.
Hi, I have read quite a few posts here about coverage already. But I still had a few questions. I have a BAM file I'm trying to find the coverage of it (typically like say 30X). So, I decided to use...
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 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 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 ArticleExtract coverage per feature from a bam and bed to a file
Hi,  a simple task.. or should be. I need to extract the average coverage per feature in a bam file. I have a genbank and bed file for the reference the bam was mapped to. if I map with e.g. Geneous...
View ArticleClik here to view.
Profile Coverage Of Rnaseq Samples?
Hi all,I have a quick question:How can I visualize aligned paired-end reads from RNAseq datasets in UCSC browser?I already mapped the reads and assembled the transcripts with Tophat/Cufflinks but I'm...
View ArticleRenaming SNPs or SNP matching
This should be easy to do by now, but... we have SNP data from an Illumina exome array given to us in PLINK format. The BIM file looks like this:1 exm2253575 0 881627 G A 1 exm269 0 881918 A G 1 exm340...
View ArticleGetting All Reads That Align To A Region In Compact Bed Format Using Bedtools?
I'm trying to find all the reads (by name) from a BAM file that align to various regions in a bed file. Right now I can do this with bedtools using intersectBed:intersectBed -abam reads.bam -wo -f 1 -b...
View Articlebedtools 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 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 ArticleSnps Comparison
Hello,I would like to compare SNPs from different methods:number of SNPsSNPs postion (position where method A has SNPs but not B and vice versa. Where both have SNPs)I would be interested to get a...
View ArticleHow To Install Bedtools In A User Directory
I am trying to install Bedtools in a user directory, however I looked at the manual for its makefile, and there is no such argument like "--prefix" for me to change. Is there a way to install all...
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 ArticleCreating Bed File For Lncrna Using Gencode Gtf File
Hi all,I want to get the bed file of lncRNA based on GENCODE GTF fileI download the file "gencode.v16.long_noncoding_RNAs.gtf.gz", and extract the chr, start, end info from the file, then I use...
View ArticleClik here to view.
Converting Bam To Bedgraph For Viewing On Ucsc?
I'm trying to go from a BAM file to a representation viewable in UCSC, ideally bedGraph. I am trying to use Bedtools's genomeCoverage like this: genomeCoverageBed -ibam accepted_hits.sorted.bam -bg...
View Article