What are the steps to extract a particular sequence from a WGS raw read

I am new to bioinformatics and recently I started a research project to compare promoters of a particular gene (not decided yet). I have several raw reads of WGS deposited in ENA site. How should I use these raw reads to extract a particular region? For example, sequence of the rice LEAFY gene.

genome next-gen sequencing

• 3,638 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by rus2dil •

1 answer

Download the data and align it against the rice reference genome. Then use the coordinates of the gene you are looking for and extract those sequences from the BAM file using samtools.

• 0 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by Ashutosh Pandey

Thank you very much Ashutosh Pandey. I think using galaxy platform I can do this? Right?

• 0 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by rus2dil •

Yes, upload FASTQ, align to ref genome. Then use SAM manipulation tools and a BED file to extract the region.

• 0 views • link

written 11.5 years ago by Ram

Thank you RamRS. Should this BED file only contain the sequence of interest? If so, how to create it from FASTQ or FATSA file. Because I have the FASTA file of the gene of interest downloaded from the reference genome.

• 0 views • link

updated 4.3 years ago by Ram

• written 11.5 years ago by rus2dil •

BED is for co-ordinates of your gene of interest. Use UCSC genome browser and beware the index differences in tools (Check out this post for more details)

• 0 views • link

written 11.5 years ago by Ram

More posts like this

Prokka annotation
written 2.1 years ago by Cancerplse •

Hey everyone! I was using Prokka recently to find unique genes for several different species from the same genus. I extracted sequences from the NCBI …
Database of Embryo-Lethal or Known Knockout Phenotype Genes in Rice
written 3.5 years ago by sviatoslav.kendall

I am looking for a database of rice genes with known knock-out phenotypes. The only example I know of something like this existing is the …
GSEA with Clusterprofiler
written 6.2 years ago by Van9 •

Recently I started a project in which I have to do bioinformatics and I am new at this. I have to do a GSEA analysis …
HMP Data Portal Samples vs Files
written 6.6 years ago by dinz_412 •

I'm quite new to Bioinformatics. I recently started familiarizing myself with [the Human Microbiome Project Data Portal](https://portal.hmpdacc.org/). In the portal,there are samples and files listed …
Forum: Up-to-date RNA-Seq Analysis Training/Courses/Papers
written 7.9 years ago by kullboys •

Hello everyone I am a PhD student with biology background. I recently inherit a RNA Sequencing project from another PhD student in my lab. We …
Finding zebrafish promoters ENSEMBL
written 8.4 years ago by Seigfried

Hello I have a list of ENSEMBL zebrafish genes whose promoters I want to extract and they are way too many to manually do it. …
From aligened sliced and filtered bam to single sequence in fasts format
written 11.4 years ago by rus2dil •

Hello, I have mapped 6 paired illumina raw reads against a reference genome and filtered by mapping quality and chromosome number 6. Then sliced using …
Extracting the sequence of the promoter of a gene using WGS raw data - Need to learn whole process
written 11.5 years ago by rus2dil •

I have WGS unaligned raw reads. I need to extract the sequence of the promoter region of a particular gene for further analysis. I would …
Finding Repeat Genes In Unused Reads
written 14.6 years ago by Lee Katz

<p>Hi, I would like to find the sequence of a repeat gene in my WGS reads. I have raw reads from both 454 and Illumina, …
Where Can I Get Whole Human Genomes?
written 15.7 years ago by Jai Puri •

<p>Hi , </p> <p>i am an Bioinformatics Master Student and i want to do a little bit of Research on my own. I decided to …