Opened 11 years ago
Last modified 11 years ago
#596 closed enhancement
Import adapter peak information when registering libraries — at Version 1
Reported by: | Nicklas Nordborg | Owned by: | Nicklas Nordborg |
---|---|---|---|
Priority: | major | Milestone: | Reggie v2.15.2 |
Component: | net.sf.basedb.reggie | Keywords: | |
Cc: |
Description (last modified by )
The library registration wizard currently import values for the library size (in 180-600BP) region and concentration measurements (from Qubit). This information is used to calculate the molarity (stored as CA_Molarity annotation).
While this is fine most of the time it can yield incorrect results (too low molarity) if the adapter peak (which is in the 120-140BP region) is high compared to the library peak.
When calculating the final molarity the size used should be a weighted average between the average adapter size and the average library size. The lab procedure has now been updated to include more information in the Caliper well table file:
Region[120-140] Size [BP]
: The average size of the adaptorsRegion[120-140] Molarity (nmol/l)
: The molarity of the adaptors (must be scaled by a factor of 10!)Region[180-600] Size [BP]
: The average size of the library fragments (this is the same column as before)Region[180-600] Molarity (nmol/l)
: The molarity of the library (must be scaled by a factor of 10!)
In principle the Qubit concentration values are not needed, but since we believe that they have higher accuracy, we use the molarity values from Caliper only to calculate a weighted size and fraction of adapters:
S(w) = S(a)*M(a) + S(l)*M(l) / (M(a)+M(l)) F(a) = M(a) / (M(a)+M(l))
Then, we use the weighted size and the Qubit concentration to calculate the final molarity as before.
We need some new annotation types:
CA_Size
: Measured size from caliper (Region[180-600] Size [BP]
)CA_Size_Adpt
: Measured size from caliper (Region[120-140] Size [BP]
)CA_Molarity
: Measured molarity from caliper (Region[180-600] Molarity (nmol/l)
). NOTE! This conflicts with the existing annotation that is the calculated value but this should be stored in a new annotation instead.CA_Molarity_Adpt
: Measured molarity from caliper (Region[120-140] Molarity (nmol/l)
).Library_Size_Est
: Weighted sizeLibrary_Molarity_Est
: Calculated molarity (the old CA_Molarity annotation)Library_Frac_Adpt
: Fraction of adapter molecules in percent.
This information will then be used in the pooling wizard to detect problems. See other ticket #597.