Skip to content

Commit a53f0fa

Browse files
committed
added list of checks to the seqspec tool
1 parent a6f633d commit a53f0fa

File tree

1 file changed

+22
-1
lines changed

1 file changed

+22
-1
lines changed

docs/SEQSPEC_TOOL.md

+22-1
Original file line numberDiff line numberDiff line change
@@ -64,7 +64,28 @@ seqspec check [-h] [-o OUT] yaml
6464
- optionally, `-o OUT` can be used to write the output to a file.
6565
- `yaml` corresponds to the `seqspec` file.
6666

67-
A list of possible errors are shown below:
67+
A list of checks performed:
68+
69+
1. Check that modalities are unique
70+
2. Check that `region_id`s of he first level of the `library_spec` correspond to modalities (and that there is one `region` per `modality`).
71+
3. Check that the onlist files exist (either locally or remotely)
72+
4. Check that the `read_id`s in the `sequence_spec` are unique.
73+
5. Check that Read files exist (either locally or remotely).
74+
6. Check that Read `primer_id`s, strand pairs are unique across all reads.
75+
7. Check that the `region_id`s are unique across all regions in the `library_spec`.
76+
8. Check that the read modalities are in assay list of modalities.
77+
9. Check that the `primer_id`s for the reads in the `sequence_spec` exist as a `region_id`s in the `library_spec`.
78+
10. Check that the `primer_id`s for the reads exist as `region_ids` for the "leaves" of the `library_spec`.
79+
11. Check `sequence_type` and region annotation consistencies:
80+
81+
- if a region has a sequence type "fixed" then it should not contain subregions
82+
- if a region has a sequence type "joiend" then it should contain subregions
83+
- if a region has a sequence type "random" then it should not contain subregions and `sequence` should be all X's
84+
- if a region has a sequence type "onlist" then it should have an onlist object
85+
86+
12. Check that the `min_len` is less than the `max_len`.
87+
13. Check that the lengths of the `sequence` in every region object is between the `min_len` and `max_len` length.
88+
14. Check that the number of files in each `File` object for all `Read` objects are the same length.
6889

6990
```bash
7091
# The "assay" value was not specified in the spec

0 commit comments

Comments
 (0)