Mercurial > repos > estrain > microrunqc

comparison median_size.py @ 2:25a92dfb780a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author estrain
date Fri, 10 Mar 2023 13:22:58 +0000
parents a53acd38d77e
children
comparison
equal deleted inserted replaced
1:2dc1074c9a6b 2:25a92dfb780a
1 #!/usr/bin/env
2
3 ## Errol Strain (estrain@gmail.com)
4 ## calculate median insert size from sam file
5
6 import numpy as np
7
8 def get_data(infile):
9 lengths = []
10 for line in infile:
11 if line.startswith('@'):
12 pass
13 else:
14 line = line.rsplit()
15 length = int(line[8])
16 if length > 0:
17 lengths.append(length)
18 else:
19 pass
20 return lengths
21
22 if __name__ == "__main__":
23 import sys
24 lengths = get_data(sys.stdin)
25 md = int(np.median(lengths))
26 print(md)