File tree Expand file tree Collapse file tree 2 files changed +4
-2
lines changed
Expand file tree Collapse file tree 2 files changed +4
-2
lines changed Original file line number Diff line number Diff line change 11from Bio import SeqIO
2+ from importlib .metadata import version
23
34
45def fasta_header (args ):
56 with open (args .filename ) as fh :
67 rec = list (SeqIO .parse (fh , "fasta" ))
78
89 fasta_format = "fasta-2line" if args .linearise_fasta else "fasta"
10+ artic_version = version ("artic" )
911
1012 with open (args .filename , "w" ) as fh :
1113 for record in rec :
1214 chrom = record .id
13- record .id = f"{ args .samplename } / { chrom } /ARTIC/ { args . caller } "
15+ record .id = f"{ args .samplename } { chrom } _artic-network/fieldbioinformatics_ { artic_version } "
1416
1517 SeqIO .write (record , fh , fasta_format )
1618
Original file line number Diff line number Diff line change @@ -121,7 +121,7 @@ def run(parser, args):
121121 raise SystemExit (1 )
122122
123123 ## collect the primer pools
124- pools = set ([row ["PoolName" ] for row in read_bed_file (bed )])
124+ pools = set ([row ["PoolName" ] for row in read_bed_file (bed )] + [ "unmatched" ] )
125125
126126 ## create a holder to keep the pipeline commands in
127127 cmds = []
You can’t perform that action at this time.
0 commit comments