Skip to content
Snippets Groups Projects
Commit fed7c783 authored by Holly Ruess's avatar Holly Ruess
Browse files

fix pool and pseudo

parent 7022e417
Branches
Tags
1 merge request!12Resolve "Fix Pool and Pseudoreps"
......@@ -133,7 +133,7 @@ def self_psuedoreplication(tag_file, prefix, paired):
# Convert read pairs to reads into standard tagAlign file
for i, index in enumerate([0, 1]):
string_index = '0' + str(index)
string_index = '.0' + str(index)
steps = ['cat %s' % (splits_prefix + string_index)]
if paired:
steps.extend([r"""awk 'BEGIN{OFS="\t"}{printf "%s\t%s\t%s\tN\t1000\t%s\n%s\t%s\t%s\tN\t1000\t%s\n",$1,$2,$3,$9,$4,$5,$6,$10}'"""])
......@@ -272,7 +272,9 @@ def main():
# Update tagAlign with single end data
if paired:
design_new_df['tag_align'] = design_new_df['se_tag_align']
design_new_df.drop(labels='se_tag_align', axis=1, inplace=True)
# design_new_df.drop(labels='se_tag_align', axis=1, inplace=True)
design_new_df.drop(labels='tag_align', axis=1, inplace=True)
# Check controls against cutoff_ratio
# if so replace with pool_control
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment