More posts like this
-
Extracting sequences from tBlastn results(BLAST standalone)
written by bpaddy •I am using tblastn to search protein queries against custom eukaryotic genome assemblies that lack protein annotations. One suggestion was to parse the hit coordinates …
-
What is a literal sequence in context of freebayes?
written by jernej.kralj96 •Hello, I need to describe the program called freebayes in my materials and methods section. I would like to know if anyone can explain what …
-
Obtaining nucleotide sequences from a set of protein identifiers
written by hopf •I've done a bit of searching but haven't seen this specific problem raised (or answered) before. Apologies if it has. I'm interested in a particular …
-
How to compare different groups p value and fold change
written by F.eghbalpour.amt •Hi dear all, I am trying to compare my RNA-seq data sets.how can I compare expression profile between groups as following: Case treated, Case untreated, …
-
discovering novel regulators
written by jiwpark00Hi, So maybe this is too far-fetched question but I have a question. I was thinking about whether or not this is possible...Let's say you …
-
How to make my nucleotide sequence as a BLAST database against my targeted sequence
written by suhalsuhaimi •hi all, I am currently working with consensus sequence of scaffold/matrix attachment regions (S/MARs) from Chinese hamster ovary (CHO) cells. Is it possible to make …
-
Open Source Or Free Pharmacophore Screening Software
written by Flow<p>I have used MOE, Discovery Studio and Ligand Scout for Pharmacophore screening of large databases using clusters and was happy with that. But now I …
-
Getting Nucleotide Based From Ucsc Or Ensembl?
written by Ashwin •<blockquote> <p><strong>Possible Duplicate:</strong> <br/> <a href="http://biostar.stackexchange.com/questions/56/how-to-get-the-sequence-of-a-genomic-region-from-ucsc">How to get the sequence of a genomic region from UCSC?</a> </p> </blockquote> <p>Hi all, I wanted to know if …
-
Variants Without Rs#
written by Ashwin •<p>Hi all, I am performing variant detection using SAMTools and 1K Genomes data and generating a VCF file. The issue is that i have several …
-
Blast Database Sequences Length
written by Fungsten •<p>Hi all</p> <p>I am doing blast queries using genomes as databases. I would like to know if is possible to filter out database sequences by …
Do you mean "decoy"?
Even if that is meant to read 'decoy' I'm not sure the question is answerable
Found little abstract about it on this page http://www.matrixscience.com/help/decoy_help.html. It loosely suggests DECOY search would lead you to screen your alignment results for false positives. It also gives perl implementation of the method(http://www.matrixscience.com/downloads/decoy.pl.gz). But this is been done with peptide sequences. All it does is it checks for random sequences & reverse of the targets & screens them accordingly. Was wondering if it makes sense to impose same rules for DNA sequence search??
This IS a real question. Decoy sequences were introduced by Heng Li and Richard Durbin and the guys at the Broad as part of the 1000 Genomes project. The idea is to add extra sequence (satellite, known BACs, HuRef) to the reference to absorb reads with do not truly belong/map to the reference, but are often mapped there by mappers because that's the best place they can put them. Heng has shown this has the effect of reducing false SNP calls by some (small) percentage. - it's not generally possible to translate back to a reference coordinate, and I dont know what you mean by screening