Mercurial > repos > triasteran > ribogalaxy_umi_processing
comparison UMI_riboseq_processing/UMI.py @ 2:6958515efa76 draft
Uploaded
author | triasteran |
---|---|
date | Mon, 20 Jun 2022 07:27:23 +0000 |
parents | 5d0d5933d370 |
children | d27375bc4a1c |
comparison
equal
deleted
inserted
replaced
1:5d0d5933d370 | 2:6958515efa76 |
---|---|
18 header = lines[0] | 18 header = lines[0] |
19 seq = lines[1] | 19 seq = lines[1] |
20 sep = lines[2] | 20 sep = lines[2] |
21 qual = lines[3] | 21 qual = lines[3] |
22 trimmed_seq = seq[2:-5]+"\n" # fooprint + barcode | 22 trimmed_seq = seq[2:-5]+"\n" # fooprint + barcode |
23 UMI = seq[0:2]+seq[-5:] #7nt in total; 5'NN and last 3'NNNNN | 23 UMI = seq[0:2]+seq[-5:].strip('\n') #7nt in total; 5'NN and last 3'NNNNN |
24 split_header = header.split(" ") | 24 split_header = header.split(" ") |
25 new_header = split_header[0]+"_"+UMI+" "+split_header[1] | 25 new_header = split_header[0]+"_"+UMI+" "+split_header[1] |
26 if qual[-1:] == "\n": | 26 if qual[-1:] == "\n": |
27 new_qual = qual[2:-5]+"\n" | 27 new_qual = qual[2:-5]+"\n" |
28 else: | 28 else: |