cutseq 0.0.44__tar.gz → 0.0.45__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
|
@@ -344,12 +344,13 @@ def pipeline_single(input1, output1, short1, untrimmed1, barcode, settings):
|
|
|
344
344
|
# step 7: trim polyA
|
|
345
345
|
if settings.trim_polyA:
|
|
346
346
|
pA_max_errors = 0.15
|
|
347
|
+
pA_max_length = 100
|
|
347
348
|
if barcode.strand == "+":
|
|
348
349
|
modifiers.append(
|
|
349
350
|
AdapterCutter(
|
|
350
351
|
[
|
|
351
352
|
NonInternalBackAdapter(
|
|
352
|
-
sequence="A" *
|
|
353
|
+
sequence="A" * pA_max_length, max_errors=pA_max_errors
|
|
353
354
|
)
|
|
354
355
|
]
|
|
355
356
|
)
|
|
@@ -359,7 +360,7 @@ def pipeline_single(input1, output1, short1, untrimmed1, barcode, settings):
|
|
|
359
360
|
AdapterCutter(
|
|
360
361
|
[
|
|
361
362
|
NonInternalFrontAdapter(
|
|
362
|
-
sequence="T" *
|
|
363
|
+
sequence="T" * pA_max_length, max_errors=pA_max_errors
|
|
363
364
|
)
|
|
364
365
|
]
|
|
365
366
|
)
|
|
@@ -580,20 +581,21 @@ def pipeline_paired(
|
|
|
580
581
|
# step 7: trim polyA
|
|
581
582
|
if settings.trim_polyA:
|
|
582
583
|
pA_max_errors = 0.15
|
|
584
|
+
pA_max_length = 100
|
|
583
585
|
if barcode.strand == "+":
|
|
584
586
|
modifiers.append(
|
|
585
587
|
(
|
|
586
588
|
AdapterCutter(
|
|
587
589
|
[
|
|
588
590
|
NonInternalBackAdapter(
|
|
589
|
-
sequence="A" *
|
|
591
|
+
sequence="A" * pA_max_length, max_errors=pA_max_errors
|
|
590
592
|
)
|
|
591
593
|
]
|
|
592
594
|
),
|
|
593
595
|
AdapterCutter(
|
|
594
596
|
[
|
|
595
597
|
NonInternalFrontAdapter(
|
|
596
|
-
sequence="T" *
|
|
598
|
+
sequence="T" * pA_max_length, max_errors=pA_max_errors
|
|
597
599
|
)
|
|
598
600
|
]
|
|
599
601
|
),
|
|
@@ -605,14 +607,14 @@ def pipeline_paired(
|
|
|
605
607
|
AdapterCutter(
|
|
606
608
|
[
|
|
607
609
|
NonInternalFrontAdapter(
|
|
608
|
-
sequence="T" *
|
|
610
|
+
sequence="T" * pA_max_length, max_errors=pA_max_errors
|
|
609
611
|
)
|
|
610
612
|
]
|
|
611
613
|
),
|
|
612
614
|
AdapterCutter(
|
|
613
615
|
[
|
|
614
616
|
NonInternalBackAdapter(
|
|
615
|
-
sequence="A" *
|
|
617
|
+
sequence="A" * pA_max_length, max_errors=pA_max_errors
|
|
616
618
|
)
|
|
617
619
|
]
|
|
618
620
|
),
|
|
File without changes
|
|
File without changes
|