diff --git a/metaSNV.py b/metaSNV.py index 2b590c98be30432f66f1be948c234eb7a1a50749..64be75aa8a822c5aa7956bf54a3d1ed0a3312c8a 100755 --- a/metaSNV.py +++ b/metaSNV.py @@ -82,9 +82,10 @@ def get_header(args): line = line.rstrip().split('\t') if len(line) != 3: continue - line[1] = line[1].replace('SN:','') - line[2] = line[2].replace('LN:','') - f.write(line[1]+'\t1\t'+line[2]+'\n') + if line[0] == "@SQ": + line[1] = line[1].replace('SN:','') + line[2] = line[2].replace('LN:','') + f.write(line[1]+'\t1\t'+line[2]+'\n') f.close() args.ctg_len = args.project_dir + '/bed_header' @@ -270,4 +271,3 @@ SOLUTION: make\n\n'''.format(basedir)) if __name__ == '__main__': main() - diff --git a/src/collapse_coverages.py b/src/collapse_coverages.py index f1d8f2fc3e4557531d9bd2f68370f7f0050b08d5..2ba7b71882a541ac629d6fa0b7f0336a31af9a8c 100644 --- a/src/collapse_coverages.py +++ b/src/collapse_coverages.py @@ -28,7 +28,7 @@ def write_matrix(cov, header, ofile): out.write('TaxId\t') out.write('\t'.join([header for _ in bamfiles])) out.write('\n') - for taxid in sorted(avg_cov.keys(), key=int): + for taxid in sorted(avg_cov.keys()): c = cov[taxid] out.write('{}\t'.format(taxid)) out.write('\t'.join([c[bf] for bf in bamfiles]))