__version__ = "v1.0"
__copyright__ = "Copyright 2022"
__license__ = "MIT"
__lab__ = "Adam Cribbs lab"
[docs]class umi(object):
def __init__(self, ):
self.umi_unit_len_fixed = 12
self.umi_unit_pattern = 1
self.pcr_errs, self.seq_errs = self.errors()
print(self.pcr_errs)
print(self.seq_errs)
[docs] def errors(self, ):
pcr_errs = []
seq_errs = []
e = 1e-5
while e < 3e-1:
pcr_errs.append(e)
seq_errs.append(e)
if 5 * e < 3e-1:
pcr_errs.append(2.5 * e)
pcr_errs.append(5 * e)
pcr_errs.append(7.5 * e)
seq_errs.append(2.5 * e)
seq_errs.append(5 * e)
seq_errs.append(7.5 * e)
e = 10 * e
pcr_errs.append(0.2)
seq_errs.append(0.2)
pcr_errs.append(0.3)
seq_errs.append(0.3)
# print(pcr_errs)
# print(seq_errs)
return pcr_errs, seq_errs
[docs] def pcrNums(self, ):
pass
[docs] def pcrErrs(self, ):
pass
[docs] def seqErrs(self, ):
pass
[docs] def umiLens(self, ):
pass
[docs] def amplRates(self, ):
pass
if __name__ == "__main__":
p = umi()
# print(p.pcrNums())
# print(p.pcrErrs())
# print(p.seqErrs())
# print(p.umiLens())
# print(p.amplRates())