bedtools no strand information - fasta file · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

bedtools no strand information - fasta file

Hello,

I want to extract some regions based on my bed file where I have chr, start, end (no strand information).

I am using bedtools getfasta, but if I don't specify the -s option, does it give me in output the reverse complement if the sequence is on the - strand?

thank you in advance

Best

bedtools fasta sequence

• 1,312 views • link

updated 3.1 years ago by Ram

• written 3.6 years ago by pingu77 •

1 answer

The manual says that the strand information is ignored by default (unless you specify -s). So if you do not have strand information, the sequences will be returned in forward orientation. To be 100% sure it is always convenient to run a small test with a dummy fasta sequence and dummy bed file (as the examples you can see in the manual).

• 0 views • link

written 3.6 years ago by iraun

Log in to answer this question.

More posts like this

How to replace the sequence in the genome.fasta with reverse complement in specific chromosome and …
written 2.6 years ago by Ap1438 •

#!/bin/bash # Input BED file with regions from multiple chromosomes input_bed="wall_region.bed" # genome file genome="Wall_10_chr.fasta" # Iterate through unique chromosomes in the BED file for …
How to correctly use bedtools merge for annotated .bed files?
written 3.1 years ago by beginnertobioinformatics •

Hello, I have two annotated .bed files that each contain 26 columns-- the first 3 columns are the standard chr number, start position, and end …
Bedtools Merge -s strandedness option
written 3.8 years ago by Jendo •

Hello, I am trying to combine features using bedtools merge. I have been using the -s option to enforce strandedness when combining the features. I'd …
Subset genomic intervals by other intervals
written 5.2 years ago by Filago •

Hello, I have a bed-like dataframe in R and want to subset it by another bed-like file e.g.: Bed file 1: Chr Start End 1 …
Bedtools maskFastaFromBed does not mask on the minus strand
written 6.5 years ago by Biostar •

Hello, I have an unmasked genome and its repeat file in bed/gff format. I tried to softmask the genome using maskFastaFromBed from bedtools (v2.26), but …
Fasta extraction from bed file
written 7.0 years ago by baurumon •

hello, how can i extract fasta in reverse order. i have bed file where start position is greater than the stop. this could come from …
How to get regions with mapped reads
written 7.8 years ago by nanoide

Hello all, So I would like to get the sequence of the genomic regions with mapped reads after an RNA-seq analysis. First, I have mapped …
Bedtools with name and full header
written 10.4 years ago by bharata1803

Hello, So, I have a bed file and a fasta reference file. I want to generate a fasta file corresponding to the bed file. I …
how to retreive sequences given a start and end site.
written 11.7 years ago by Affan

So I have a gff that has information about MEF2 transcription factor binding sites. So given a start and end site, 19815641 - 19815654 on …
problem with bedtools complement: how to extract strand information?
written 11.8 years ago by biolab

Hi everyone, I have a minor trouble when using bedtools complement. My gff3 file contains Chr start end strand etc. information. I want to use …