code: def blastnprep(self):
          
          # Populate threads for each gene, genome combination
          for sample in self.metadata:
              if sample.general.bestassemblyfile != 'NA':
                  #
                  # sample[self.analysistype].alleleresults = GenObject()
                  sample[self.analysistype].closealleles = dict()
                  sample[self.analysistype].mismatches = dict()
                  sample[self.analysistype].alignmentlength = dict()
                  sample[self.analysistype].subjectlength = dict()
                  sample[self.analysistype].queryid = dict()
                  sample[self.analysistype].start = dict()
                  sample[self.analysistype].end = dict()
                  sample[self.analysistype].queryseq = dict()
                  if type(sample[self.analysistype].allelenames) == list:
                      for allele in sample[self.analysistype].combinedalleles:
                          # Add each fasta/allele file combination to the threads
                          self.runblast(sample.general.bestassemblyfile, allele, sample)
docs: Setup blastn analyses