NGS Library Molarity & Equimolar Pooling Calculator (Qubit ng/µL to nM)
Core Formula
Library molar concentration is determined jointly by the Qubit mass concentration and the average fragment length:
c (nM) = c_mass (ng/µL) × 10^6 / (660 × L_bp)
c_mass: mass concentration from Qubit fluorometric quantification (dsDNA HS kit), in ng/µLL_bp: average library fragment length from Bioanalyzer or TapeStation, in bp, including adapters660: average molecular weight of dsDNA per base pair, in g/mol
Derivation: dsDNA molecular weight = 660 × L g/mol; converting ng/µL to g/L requires multiplying by 10⁻³; dividing by molecular weight gives mol/L, then multiplying by 10⁹ converts to nM, which simplifies to the formula above.
Worked example (Lib-A, 5.20 ng/µL, 385 bp): 5.20 × 10⁶ / (660 × 385) = 5,200,000 / 254,100 ≈ 20.5 nM.
Why Do You Need Both Mass Concentration and Fragment Length?
Molarity measures the number of molecules, not mass. Two libraries both at 5 ng/µL — the 200 bp library contains roughly twice as many molecules as the 400 bp library. Pooling by equal mass gives the shorter library twice the template molecules, resulting in a two-fold difference in sequencing depth — the most common cause of failed equimolar pooling.
On “Average Fragment Length”
You must use the complete library length including adapters, i.e., the average peak from the Bioanalyzer electropherogram, not the insert size. Illumina standard TruSeq adapters add approximately 120 bp to the insert (Y-adapter total on both ends): a 250 bp insert library has a complete fragment of approximately 370 bp. Enter the number reported by the instrument directly — no manual addition or subtraction needed.
Equimolar Pooling Protocol
The standard workflow has two steps:
- Dilute all libraries to the same working concentration W (typically 10 nM) using TE buffer or 10 mM Tris-HCl (pH 8.0):
V_TE = V_stock × (c / W − 1)
- Take equal volumes from each diluted library and combine — this gives the equimolar pool.
If a library’s stock concentration is already below the target working concentration and cannot be diluted to it, lower the target working concentration until all libraries can reach the same level.
Illumina Loading Concentration Reference
After pooling, further dilute to the loading concentration (after alkaline denaturation to single-stranded DNA):
| Instrument | Typical loading concentration |
|---|---|
| NovaSeq 6000 | 300–400 pM |
| NextSeq 500/550 | 1.5–2.0 pM |
| MiSeq | 6–20 pM |
Follow the Denature and Dilute Guide for each instrument; the table above is for reference only.
Common Mistakes
Using NanoDrop instead of Qubit: NanoDrop relies on A260 absorbance and cannot distinguish target dsDNA from free nucleotides or RNA carry-over, typically overestimating by 10–50%, leading to lower-than-intended loading concentration and insufficient cluster density.
Entering insert size instead of library fragment length: Systematically underestimates library concentration by ~30%, causing shorter-fragment libraries to be over-represented when pooling by molar ratio.
Applicability
- For dsDNA libraries only (Illumina TruSeq, TruSeq Nano, Nextera, and other dsDNA library prep products)
- 660 g/mol/bp is an approximation averaged across all four bases; for libraries with extreme GC bias (< 30% or > 70% GC) the error is approximately 1–2%, negligible for pooling decisions
- When the fragment length distribution is broad (e.g., ATAC-seq with multiple peaks), use the mode of the main peak rather than the arithmetic mean for more accurate results
- Not applicable to ssDNA libraries (small RNA sequencing) or direct RNA quantification
FAQ
Two libraries have the same ng/µL — why are their nM values different?
Because molarity measures the number of molecules, and each molecule's mass depends on its length. At the same 5 ng/µL, a 200 bp fragment library contains roughly twice as many molecules as a 400 bp library. Pooling by equal mass (ng) will give the shorter library twice the sequencing depth.
Should I subtract the adapter length from the Bioanalyzer fragment size?
No. The average peak reported by Bioanalyzer or TapeStation is already the complete library length including adapters — enter it directly. Illumina TruSeq adapters add approximately 120 bp to the insert, and that is already included in the reported value.
Why use Qubit for quantification instead of NanoDrop?
NanoDrop measures A260 absorbance and cannot distinguish target double-stranded DNA from free nucleotides or RNA carry-over, often overestimating by 10–50%. Use the Qubit dsDNA HS kit for library quantification.
Where does the factor 660 come from?
The average residue mass of the four deoxyribonucleotides is approximately 309 g/mol/nt; two complementary strands total approximately 618 g/mol/bp, rounded to 660 g/mol/bp after hydration correction — consistent with Illumina technical documentation and standard practice in the field.
What if a library's concentration is below the target working concentration?
A library below the target cannot be diluted to that concentration; the tool will flag it with a warning. Lower the target working concentration (e.g., from 10 nM to 4 nM) until all libraries can reach the same concentration before pooling.
Related tools
SDS-PAGE Gel Preparation Calculator (Resolving Gel & Stacking Gel)
Enter acrylamide percentage and per-gel volume; calculate acrylamide, Tris, SDS, APS, TEMED, and water for both gel layers, with multi-gel batch scaling.
Transfection Volume Calculator — Lipofectamine 3000 / PEI Scaled by Growth Area
Scale Lipofectamine 3000 or PEI transfection volumes by growth area; gives a master mix with 10% overage that reproduces manufacturer protocol table values.
Centrifuge RCF (g) ↔ RPM Converter — Relative Centrifugal Force Calculator
Convert between relative centrifugal force (× g) and rotor speed (rpm) in both directions using rotor radius, with common reference values.
Buy me a coffee