#743 closed task (performed)
Implement tool for SAI filtered sequence export
Reported by: | westram | Owned by: | westram |
---|---|---|---|
Priority: | major | Milestone: | arb7.0 |
Component: | External_tools | Version: | SVN |
Keywords: | silva | Cc: |
Description (last modified by westram)
Specification:
- load ARB database
- filter sequences by existing SAI(s)
- export remaining columns of all sequences to FASTA file
- skip sequences with < --min-bases base positions
- logged skipped species to stderr (using field 'name')
CLI :
Switches……………………………. | Description |
---|---|
--db "database.arb" | name of ARB input database |
--ali "aliname" | name of alignment in ARB |
--fasta "flat.file" | name of FASTA output file |
--id "ACI" | ACI defining FASTA header line (default: "readdb(name)") (see ACI manual) |
--min-bases NUM | do not export sequences with < NUM bases left after applying filter |
--count-bases "ACGTU" | list of base-characters to count for --min-bases (default: "A..Z") |
(the following arguments may be specified multiple times) | |
--filterby "SAIname" | name of SAI used to filter sequence data |
--pass [allbut] "chars" | characters that will forward column to output if found in SAI |
--block [allbut] "chars" | characters that will block column from output if found in SAI |
- --pass and --block may only occur after and always apply to the last preceding --filterby.
- the optional parameter [allbut] inverts the character set, i.e. if only one filter is specified the following are equivalent:
- --pass "ABC" ⇔ --block allbut "ABC"
- --block "XYZ" ⇔ --pass allbut "XYZ"
Filter combination:
OR | --filterby "PV1" --pass allbut ".-=0123" --filterby "PV2" --pass allbut ".-=012345" |
AND | --filterby "PV1" --block ".-=0123" --filterby "PV2" --block ".-=012345" |
Assuming the SAIs contain positional variability:
- OR: will export all columns where PV1 is ≥ 4 or PV2 is ≥ 6
- AND: will export all columns where PV1 is ≥ 4 and PV2 is ≥ 6
Change History (10)
comment:1 Changed 8 years ago by westram
- Description modified (diff)
comment:2 Changed 8 years ago by westram
- Description modified (diff)
comment:3 Changed 8 years ago by westram
- Description modified (diff)
comment:4 Changed 8 years ago by westram
- Milestone changed from wishlist to r17q2
- Owner changed from devel to westram
- Status changed from new to assigned
comment:5 Changed 8 years ago by westram
- Status changed from assigned to _started
comment:6 Changed 7 years ago by westram
- Description modified (diff)
- proposes add. switch --count-bases
comment:7 Changed 7 years ago by westram
- Description modified (diff)
- Resolution set to performed
- Status changed from _started to closed
by [16119]
comment:8 Changed 7 years ago by westram
- Resolution performed deleted
- Status changed from closed to _started
- add sequence post-processing
comment:9 Changed 7 years ago by westram
- Resolution set to performed
- Status changed from _started to closed
by [16234]
comment:10 Changed 3 years ago by westram
- Milestone changed from r17q2 to arb7.0
Note: See
TracTickets for help on using
tickets.