We will be offering mothur and R workshops throughout 2019. Learn more.

Get.seqs

From mothur
Revision as of 12:06, 1 July 2009 by Pschloss (Talk | contribs)

Jump to: navigation, search

Given a list of accession numbers (i.e. sequence names) and one or more file formats, generate a new file that contains only those sequences.


Options

accnos, fasta, name, group, alignreport; each takes a file name


Required

accnos and one of fasta/name/group/alignreport


Output

*.pick.*


Algorithm

  1. read accnos file into a set<string> container, close the file
  2. read through the file to be parsed and for each entry, if the sequence name is in the set<string> container:
    • spit the data out to the new file and delete the entry from the set<string> container
    • otherwise do nothing