@@ -61,12 +61,12 @@ intervals (R style) using the default setting of ``base="0"``. If
61
61
62
62
library(airr)
63
63
64
- example_data <- file.path(".. ", "tests", "data-tests ", "toy_data.tsv ")
64
+ example_data <- system.file("extdata ", "rearrangement-example.tsv.gz ", package="airr ")
65
65
basename(example_data)
66
66
67
67
::
68
68
69
- ## [1] "toy_data .tsv"
69
+ ## [1] "rearrangement-example .tsv.gz "
70
70
71
71
.. code :: r
72
72
@@ -83,27 +83,24 @@ intervals (R style) using the default setting of ``base="0"``. If
83
83
84
84
::
85
85
86
- ## # A tibble: 6 x 44
87
- ## rearrangement_id rearrangement_se… sequence_id sequence rev_comp
88
- ## <chr> <chr> <chr> <chr> <lgl>
89
- ## 1 IVKNQEJ01BVGQ6 1 IVKNQEJ01B… GGCCCAGGACTGGTG… TRUE
90
- ## 2 IVKNQEJ01AQVWS 1 IVKNQEJ01A… GGCCCAGGACTGGTG… TRUE
91
- ## 3 IVKNQEJ01AOYFZ 1 IVKNQEJ01A… GGCCCAGGACTGGTG… TRUE
92
- ## 4 IVKNQEJ01EI5S4 1 IVKNQEJ01E… GGCCCAGGACTGGTG… TRUE
93
- ## 5 IVKNQEJ01DGRRI 1 IVKNQEJ01D… GGCCCAGGACTGGTG… TRUE
94
- ## 6 IVKNQEJ01APN5N 1 IVKNQEJ01A… GGCCCAGGACTGGTG… TRUE
95
- ## # ... with 39 more variables: productive <lgl>, sequence_alignment <chr>,
96
- ## # germline_alignment <chr>, v_call <chr>, d_call <chr>, j_call <chr>,
97
- ## # c_call <chr>, junction <chr>, junction_length <int>,
98
- ## # junction_aa <chr>, v_score <dbl>, d_score <dbl>, j_score <dbl>,
99
- ## # c_score <dbl>, v_cigar <chr>, d_cigar <chr>, j_cigar <chr>,
100
- ## # c_cigar <chr>, v_identity <dbl>, v_evalue <dbl>, d_identity <dbl>,
101
- ## # d_evalue <dbl>, j_identity <dbl>, j_evalue <dbl>,
86
+ ## # A tibble: 6 x 33
87
+ ## sequence_id sequence rev_comp productive stop_codon vj_in_frame v_call
88
+ ## <chr> <chr> <lgl> <lgl> <lgl> <lgl> <chr>
89
+ ## 1 SRR765688.7… NNNNNNNN… FALSE TRUE FALSE TRUE IGHV2…
90
+ ## 2 SRR765688.3… NNNNNNNN… FALSE TRUE FALSE TRUE IGHV5…
91
+ ## 3 SRR765688.3… NNNNNNNN… FALSE TRUE FALSE TRUE IGHV7…
92
+ ## 4 SRR765688.3… NNNNNNNN… FALSE TRUE FALSE TRUE IGHV7…
93
+ ## 5 SRR765688.4… NNNNNNNN… FALSE TRUE FALSE TRUE IGHV7…
94
+ ## 6 SRR765688.1… NNNNNNNN… FALSE FALSE TRUE TRUE IGHV2…
95
+ ## # ... with 26 more variables: d_call <chr>, j_call <chr>, c_call <chr>,
96
+ ## # sequence_alignment <chr>, germline_alignment <chr>, junction <chr>,
97
+ ## # junction_aa <chr>, v_cigar <chr>, d_cigar <chr>, j_cigar <chr>,
98
+ ## # junction_length <int>, np1_length <int>, np2_length <int>,
102
99
## # v_sequence_start <dbl>, v_sequence_end <int>, v_germline_start <dbl>,
103
100
## # v_germline_end <int>, d_sequence_start <dbl>, d_sequence_end <int>,
104
101
## # d_germline_start <dbl>, d_germline_end <int>, j_sequence_start <dbl>,
105
102
## # j_sequence_end <int>, j_germline_start <dbl>, j_germline_end <int>,
106
- ## # np1_length <int>, np2_length <int>, duplicate_count <int>
103
+ ## # duplicate_count <int>
107
104
108
105
Writing AIRR formatted files
109
106
----------------------------
0 commit comments