Can BLAST perfrom a search including IUPAC ambiguity codes?

More Bradley Austin Weller's questions See All

Is there a solution to this conundrum on freefall into a black hole?

18 January 2021 1,582 3 View

Is there a way to batch-convert mass spectrometry cdf files to mzML files?

I have a set of roughly one hundred cdf files from a mass spectrometry experiment. I would like to analyse them using MetaboAnalystR, which requires mzML files as input. I know I can convert each...

28 March 2020 1,794 6 View

What is the best methodology for working with temporal data holes and deep neural networking?

Any data science and/or remote sensing experts have any insight pertaining to multivariate regression and deep neural networking with large feature set inputs and variable temporal data holes...

03 February 2020 5,246 2 View

Why are surrogates only added at the time of extraction, not at the time of sample collection and processing?

I am working on a study on pharmaceutical fate and uptake in plants. We are wondering about analyte degradation during the initial transport, processing, and storage of the soil and plant samples,...

30 October 2019 3,540 1 View

Which Protein Viewer do you prefer?

There are many different software tools available to visualize and sometimes to edit protein structure files, eg. in PDB format. I am using a Windows System (Windows 10). Thank you for any suggestion.

30 September 2019 5,046 7 View

Is there a method or reference that shows how (or what is) the oxygen transmission rate (OTR) of water?

I have found multiple references on solubility however I can find nothing on transmissibility. The later is usually reserved for polymer films but in theory it should be determinable for water as...

05 August 2019 4,776 3 View

How is depth information encoded in the frequency of spectral interferograms obtained by Spectral domain optical coherence tomography?

I am trying to understand how the depth information in SD-OCT A-scans is obtained over distances far in excess of the source coherence length. From what I have gathered so far, it seems that the...

10 June 2019 6,783 2 View

How can I improve laser penetration on a confocal microscope in aortic tissue?

I am currently doing an experiment to image thick (~100 micron) cut pieces of aorta and am having a few problems. To imagine the shape of my tissue, I am cutting rings or bands of aorta and...

01 August 2018 920 3 View

What is the expected number of connected components in an Erdos-Renyi graph?

Let G(N,p) be an Erdos-Renyi graph, where N is the number of vertices, and p is the probability that two distinct vertices form an edge. If I plot 1-b0/N over log(p), then I obtain a curve which...

23 November 2016 6,443 3 View

Database of reagent off-target effects?

I remember seeing a recent website/database that listed commonly used cell culture and research chemical agents and their known off-target effects. Can anyone help me re-discover the site?

17 April 2016 5,311 2 View

error in running a model in Abaqus

Excessive rotation of nodes in node set ErrNodeExcessRotation-Step1

17 August 2021 0 0 View

error in modelling model in Abaqus

Excessive rotation of nodes in node set ErrNodeExcessRotation-Step1

17 August 2021 0 0 View

MD MOE RMSD

How to calculate the RMSD values for a MD simulation using MOE?

07 August 2021 0 0 View

MD MOE RMSD

How to calculate the RMSD values for a MD simulation using MOE?

07 August 2021 0 0 View

MGL tools

20 June 2021 0 0 View

fatal error

When I tried to energy minimization my system, I got fatal error as below. Fatal error: Atomtype opls_116 not found Although I've already added this line: ; include water #include "oplsaa.ff/spc.itp" to [molecultype] directive in my topology.

16 June 2021 0 0 View

fatal error

16 June 2021 0 0 View

fatal error

16 June 2021 0 0 View

fatal error

16 June 2021 0 0 View

i have some problem with chtMultiRegionSimpleFoam

--> FOAM FATAL ERROR: Maximum number of iterations exceeded From function Foam::scalar Foam::species::thermo::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo = Foam::species::thermo >, Foam::sensibleEnthalpy>] in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::heRhoThermo >, Foam::sensibleEnthalpy> > > >::calculate() at ??:? #3 Foam::heRhoThermo >, Foam::sensibleEnthalpy> > > >::correct() at ??:? #4 ? at ??:? #5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #6 ? at ??:?

24 March 2021 0 0 View

Salvador Ramirez-Flandes

Hello Bradley,

Although you can include those ambiguous codes in the input for BLAST, they will be treated as mismatches in the alignment. So, they are not totally ignored, instead you can think of it as if they were gaps.

You can read more about this here:

https://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastDocs&DOC_TYPE=BlastHelp

I hope that helps.

Best regards

Abhishek Kumar

problem can be resolve by biopython alphabets

http://biopython.org/DIST/docs/tutorial/Tutorial.html

and if u look for the Substitution Matrix for IUPAC Nucleotide Ambiguity Codes - NUC.4.4

ftp://ftp.ncbi.nih.gov/blast/matrices/NUC.4.4

Bradley Austin Weller

Thank you for the responses. It turns out that BLAST will include (and 'interpret') IUPAC ambiguity codes for a nucleotide BLAST, but will discount them for a protein BLAST (where codons are needed rather than nuceotide bases) - makes sense!