Skip to content

QIAseq Immune Repertoire RNA Library kit

Here we will discuss how to process TCR cDNA libraries obtained with QIAseq Immune Repertoire RNA Library kit.

Data libraries

This tutorial uses the data from the following publication: "Tumor-infiltrating mast cells are associated with resistance to anti-PD-1 therapy.",Somasundaram R et al., Nat Commun, 2021 Jan 12;12(1):346 doi:10.1038/s41467-020-20600-7

A humanized (Hu)-mouse melanoma model was created by injecting fetal liver-derived CD34+ cells and implanting autologous thymus in immune-deficient NOD-scid IL2Rγnull (NSG) mice. Reconstituted Hu-mice were challenged with HLA-matched melanomas and treated with anti-PD-1, which resulted in restricted tumor growth but not complete regression. RNA was isolated from spleen and tumor tissues obtained from pre- and post-therapy (anti-PD-1) mice. cDNA synthesis was performed using QIAseq Immune Repertoire RNA Library prep kit (QIAGEN). Enrichment of complete TCR variable region by gene-specific primers targeting constant regions and molecular indexing (UMIs) for accurate and sensitive TCR clonotype and repertoire diversity assessment. Final libraries were sequenced on Illumina's NextSeq500 using the 300 cycles mid output sequencing kit, 261 cycles for read 1 and 41 cycles for read 2.

On the scheme bellow you can see structure of cDNA library. UMI is located in the first 12 bp of R2.

All data may be downloaded directly from SRA (PRJEB44566) using e.g. SRA Explorer.

Use aria2c for efficient download of the full dataset with the proper filenames:

download.sh
mkdir -p raw
aria2c -c -s 16 -x 16 -k 1M -j 8 -i download-list.txt
download-list.txt
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/087/SRR12991387/SRR12991387_1.fastq.gz
  out=raw/mice_tumor_1_R1.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/087/SRR12991387/SRR12991387_2.fastq.gz
  out=raw/mice_tumor_1_R2.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/088/SRR12991388/SRR12991388_1.fastq.gz
  out=raw/mice_tumor_2_R1.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/088/SRR12991388/SRR12991388_2.fastq.gz
  out=raw/mice_tumor_2_R2.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/090/SRR12991390/SRR12991390_1.fastq.gz
  out=raw/mice_spleen_2_R1.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/090/SRR12991390/SRR12991390_2.fastq.gz
  out=raw/mice_spleen_2_R2.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/089/SRR12991389/SRR12991389_1.fastq.gz
  out=raw/mice_spleen_1_R1.fastq.gz
ftp://ftp.sra.ebi.ac.uk/vol1/fastq/SRR129/089/SRR12991389/SRR12991389_2.fastq.gz
  out=raw/mice_spleen_1_R2.fastq.gz

Upstream analysis

MiXCR has a dedicated preset for this protocol, thus analysing the data is as easy as:

mixcr analyze qiagen-human-rna-tcr-umi-qiaseq \
    raw/mice_tumor_1_R1.fastq.gz \
    raw/mice_tumor_1_R2.fastq.gz \
    results/mice_tumor_1

Running the command above will generate the following files:

> ls result/

# human-readable reports 
mice_tumor_1.report
# raw alignments (highly compressed binary file)
mice_tumor_1.vdjca
# alignments with refined UMI barcode sequences 
mice_tumor_1.refined.vdjca
# TRA, TRB CDR3 clonotypes (highly compressed binary file)
mice_tumor_1.clns
# TRA,TRB CDR3 clonotypes exported in tab-delimited txt
mice_tumor_1.clonotypes.TRA.tsv
mice_tumor_1.clonotypes.TRB.tsv

While .clns file holds all data and is used for downstream analysis using mixcr postanalisis, the output .txt clonotype table will contain exhaustive information about each clonotype as well:

See first 100 records from FebControl1.clones.IGH.tsv clonotype table
cloneId cloneCount uniqueTagCountUMI cloneFraction targetSequences targetQualities allVHitsWithScore allDHitsWithScore allJHitsWithScore allCHitsWithScore allVAlignments allDAlignments allJAlignments allCAlignments nSeqFR1 minQualFR1 nSeqCDR1 minQualCDR1 nSeqFR2 minQualFR2 nSeqCDR2 minQualCDR2 nSeqFR3 minQualFR3 nSeqCDR3 minQualCDR3 nSeqFR4 minQualFR4 aaSeqFR1 aaSeqCDR1 aaSeqFR2 aaSeqCDR2 aaSeqFR3 aaSeqCDR3 aaSeqFR4 refPoints
0 869027 6008 0.150075 TGTGCCGTGAACATCAACTTCAACAAATTTTACTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(662.8) nan TRAJ21*00(262.2) TRAC*00(260.2) 264 278 297 0 14 70.0 nan 22 44 75 14 36 110.0 nan nan nan nan nan nan nan nan nan nan nan
1 306627 1903 0.0529522 TGTGCTTATAGGAGGGGCCAGTTCTATTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV38-2DV8*00(676) nan TRAJ49*00(216.7) TRAC*00(259.4) 273 287 309 0 14 70.0 nan 32 45 76 17 30 65.0 nan nan nan nan nan nan nan nan nan nan nan
2 209580 1416 0.0361929 TGTGCAGCAAGAAACGGGCAGGAGAGCACTTACTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-1*00(653.9) nan TRAJ5*00(272.8) TRAC*00(258.4) 267 278 300 0 11 55.0 nan 22 49 80 10 37 SC24A 121.0 nan nan nan nan nan nan nan nan nan nan nan
3 142220 2030 0.0245604 TGTGCCAGCAGTCCCGGGACAGGGGGCTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV27*00(658.1) TRBD1*00(70) TRBJ2-7*00(218.3) TRBC2*00(287.1) 270 282 307 0 12 60.0 9 23 36 12 26 70.0 23 39 67 26 42 80.0 nan nan nan nan nan
4 123769 831 0.021374 TGTGCTGTTAAGACAACTGACAGCTGGGGGAAATTGCAGTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV21*00(619.4) nan TRAJ24*00(307.8) TRAC*00(259.7) 267 275 299 0 8 40.0 nan 15 52 83 5 42 SC18TST20A 157.0 nan nan nan nan nan nan nan nan nan nan nan
5 106352 651 0.0183662 TGTGCCAGTAGTATCGTGGCAGGGCCTTTAACCGGGGAGCTGTTTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV19*00(564.5) TRBD2*00(31) TRBJ2-2*00(225.9) TRBC2*00(280.1) 270 284 307 0 14 70.0 25 34 48 17 26 SG27C 31.0 25 43 71 30 48 90.0 nan nan nan nan nan
6 96240 727 0.0166199 TGTGCTCTGAGTGGTAACACCGACAAGCTCATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV19*00(683.8) nan TRAJ34*00(262.2) TRAC*00(261.3) 273 286 310 0 13 65.0 nan 25 47 78 14 36 110.0 nan nan nan nan nan nan nan nan nan nan nan
7 77999 497 0.0134699 TGTGCCAGGCGGGGCTCAGGAACCTACAAATACATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV22*00(614.9) nan TRAJ40*00(276.8) TRAC*00(260.6) 258 263 291 0 5 25.0 nan 25 50 81 14 39 125.0 nan nan nan nan nan nan nan nan nan nan nan
8 67766 723 0.0117027 TGCGCCAGCAGCCAAGCGGACAGCTCTGGGGCCAACGTCCTGACTTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV4-1*00(588.9) nan TRBJ2-6*00(272.1) TRBC2*00(287.9) 270 286 307 0 16 80.0 nan 18 45 73 21 48 135.0 nan nan nan nan nan nan nan nan nan nan nan
9 60224 512 0.0104002 TGTGCTACGGAACTGAACGACTACAAGCTCAGCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(678.8) nan TRAJ20*00(257.7) TRAC*00(262.9) 264 275 297 0 11 55.0 nan 25 46 77 15 36 105.0 nan nan nan nan nan nan nan nan nan nan nan
10 44410 276 0.00766928 TGTGCCAGCAGCCAAGATCCCGACAGGGGAGCTTTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV4-2*00(660.6) TRBD1*00(40) TRBJ1-1*00(188) TRBC1*00(283.6) 270 289 307 0 19 95.0 14 22 36 21 29 40.0 30 40 68 29 39 50.0 nan nan nan nan nan
11 44001 286 0.00759865 TGCGCCAGCAGCCTGGACAGGGGGAATTCACCCCTCCACTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV4-3*00(604.9) TRBD1*00(50) TRBJ1-6*00(227) TRBC1*00(283) 270 283 307 0 13 65.0 13 23 36 14 24 50.0 27 45 73 24 42 90.0 nan nan nan nan nan
12 43498 245 0.00751179 TGTGCAGCCGGGCTCGGAAAGCTGATTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV25*00(653.7) nan TRAJ4*00(213) TRAC*00(260) 264 271 293 0 7 35.0 nan 40 52 83 18 30 60.0 nan nan nan nan nan nan nan nan nan nan nan
13 39304 239 0.00678751 TGTGCTGTGCTGGAGGAAACCAGTGGCTCTAGGTTGACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV20*00(634.9) nan TRAJ58*00(286) TRAC*00(251.6) 261 274 294 0 13 SA271T 51.0 nan 25 52 83 15 42 135.0 nan nan nan nan nan nan nan nan nan nan nan
14 36795 318 0.00635423 TGCGGCACAGCCACGGGAGGAGGAAACAAACTCACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV30*00(647.4) nan TRAJ10*00(292.1) TRAC*00(264.4) 261 271 294 0 10 50.0 nan 25 53 84 11 39 140.0 nan nan nan nan nan nan nan nan nan nan nan
15 33932 226 0.00585981 TGTGCCACCTGGGACGGGCCGGACCACTGGTTGGTTCAAGATATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV2*00(501.7) nan TRGJP1*00(284.8) TRGC1*00(352.5) 282 302 320 0 20 100.0 nan 22 46 80 22 46 120.0 nan nan nan nan nan nan nan nan nan nan nan
16 29400 200 0.00507716 TGTGCCGTGTCATCTGGTTCTGCAAGGCAACTGACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(625.7) nan TRAJ22*00(287.3) TRAC*00(258.5) 264 273 297 0 9 45.0 nan 22 52 83 9 39 ST24A 136.0 nan nan nan nan nan nan nan nan nan nan nan
17 29281 165 0.00505661 TGTGCCTTTATCGATAACTATGGTCAGAATTTTGTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV24*00(632.4) nan TRAJ26*00(285.8) TRAC*00(255.6) 267 278 297 0 11 55.0 nan 22 49 80 12 39 135.0 nan nan nan nan nan nan nan nan nan nan nan
18 28034 155 0.00484127 TGTGCAGAGAGTATAGGAGGAAGCTACATACCTACATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV5*00(647.4) nan TRAJ6*00(276.4) TRAC*00(258.4) 264 279 297 0 15 75.0 nan 26 51 82 14 39 125.0 nan nan nan nan nan nan nan nan nan nan nan
19 27314 201 0.00471693 TGTGCTGTGCAGGCCCGAGGAACCTACAAATACATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV20*00(676.5) nan TRAJ40*00(262.2) TRAC*00(256.4) 261 276 294 0 15 75.0 nan 28 50 81 17 39 110.0 nan nan nan nan nan nan nan nan nan nan nan
20 26762 363 0.0046216 TGTGCCAGCAGTTTCGGCAGGGGGTCCTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV27*00(665.2) nan TRBJ2-7*00(230.1) TRBC2*00(286.9) 270 284 307 0 14 70.0 nan 15 39 67 18 42 SA18GSC20G 92.0 nan nan nan nan nan nan nan nan nan nan nan
21 26666 167 0.00460502 TGTGCCAGCAGTTTTCTGCTCGAATCCTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV28*00(643.5) nan TRBJ2-7*00(228.1) TRBC2*00(282) 270 284 307 0 14 70.0 nan 21 39 67 24 42 90.0 nan nan nan nan nan nan nan nan nan nan nan
22 26299 289 0.00454164 TGCATCCTGAGGCCCCTATTATAACCAGGGAGGAAAGCTTATCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV26-2*00(590.9) nan TRAJ23*00(291.4) TRAC*00(260.5) 261 272 296 0 11 55.0 nan 24 52 83 18 46 140.0 nan nan nan nan nan nan nan nan nan nan nan
23 26019 340 0.00449329 TGTGCCAGCAGCTATGAGGGCCACCTTTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV7-2*00(631.8) TRBD2*00(35) TRBJ2-7*00(212.7) TRBC2*00(286.5) 273 286 310 0 13 65.0 27 34 48 15 22 35.0 24 39 67 27 42 75.0 nan nan nan nan nan
24 25677 186 0.00443423 TGTGCCGTGAAAAAAAACACCGACAAGCTCATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(633) nan TRAJ34*00(256.8) TRAC*00(260.6) 264 275 297 0 11 55.0 nan 26 47 78 15 36 105.0 nan nan nan nan nan nan nan nan nan nan nan
25 25361 173 0.00437966 TGTGCCAGCAGTCCCTTTGCGGGGGATGGAACAGATACGCAGTATTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV28*00(589.6) TRBD100(25),TRBD200(25) TRBJ2-3*00(228.2) TRBC2*00(283.5) 270 282 307 0 12 60.0 18 23 36 20 25 25.0;23 28 48 18 23 25.0 23 41 69 30 48
26 24540 145 0.00423788 TGTGCAGCACCTACGGGAGGAGGAAACAAACTCACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-1*00(631.6) nan TRAJ10*00(286.8) TRAC*00(258.3) 267 276 300 0 9 45.0 nan 26 53 84 12 39 135.0 nan nan nan nan nan nan nan nan nan nan nan
27 23553 309 0.00406743 TGTGCCAGCAGTTTTGGAGCAGGCTATGGCTACACCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV27*00(679.1) nan TRBJ1-2*00(218) TRBC1*00(285.5) 270 284 307 0 14 70.0 nan 24 40 68 23 39 80.0 nan nan nan nan nan nan nan nan nan nan nan
28 22400 170 0.00386832 TGTGCAATGAGGTTCTACATGGATAGCAACTATCAGTTAATCTGG NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV14DV4*00(588.4) nan TRAJ33*00(291.9) TRAC*00(260.5) 273 284 310 0 11 55.0 nan 18 46 77 17 45 140.0 nan nan nan nan nan nan nan nan nan nan nan
29 22239 158 0.00384051 TGTGCTGTGAACCCTTATAACCAGGGAGGAAAGCTTATCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV21*00(627.8) nan TRAJ23*00(291.4) TRAC*00(265.2) 267 277 299 0 10 50.0 nan 24 52 83 14 42 140.0 nan nan nan nan nan nan nan nan nan nan nan
30 22097 149 0.00381599 TGTGCCTCTCAGGATTCTGGGGGTTACCAGAAAGTTACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV24*00(614) nan TRAJ13*00(282.4) TRAC*00(258.8) 267 274 297 0 7 35.0 nan 23 52 83 13 42 145.0 nan nan nan nan nan nan nan nan nan nan nan
31 21877 131 0.003778 TGTGCCAGCAGTGTACTAGCGGGGGGGGACACCGGGGAGCTGTTTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV12-300(455.1),TRBV12-400(455) TRBD2*00(46) TRBJ2-2*00(236.5) TRBC2*00(285.9) 273 288 310 0 15 ST285G 61.0;273 288 310 0 15 ST285G 61.0 20 32 48 15 27 SA28G 46.0 23 43 71 28 48
32 21760 183 0.00375779 TGCGGCATCTCACAGAACACCGGTAACCAGTTCTATTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV30*00(625.2) nan TRAJ49*00(276.3) TRAC*00(260.2) 261 268 294 0 7 35.0 nan 20 45 76 14 39 125.0 nan nan nan nan nan nan nan nan nan nan nan
33 21134 121 0.00364969 TGCAGTGCTAGATGGGGGAATGCGGAGACCCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV20-1*00(578.6) TRBD1*00(25) TRBJ2-5*00(210.4) TRBC2*00(277.5) 279 291 313 0 12 60.0 18 23 36 13 18 25.0 25 40 68 24 39 75.0 nan nan nan nan nan
34 20371 294 0.00351792 TGTGCCAGCAGTGCGGGACTGCAAGAGACCCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV25-1*00(686.2) TRBD200(35),TRBD100(30) TRBJ2-5*00(228.4) TRBC2*00(291) 270 283 307 0 13 65.0 15 22 48 13 20 35.0;11 17 36 13 19 30.0 22 40 68 21 39
35 19293 175 0.00333176 TGTGCAATGAGAGAGGGCCTGGGCACCGACAAGCTCATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV14DV4*00(675.7) nan TRAJ34*00(247.2) TRAC*00(262) 273 292 310 0 19 95.0 nan 28 47 78 23 42 95.0 nan nan nan nan nan nan nan nan nan nan nan
36 18790 154 0.0032449 TGTGCCAGCAGCTTAGCTCCGGGACAACAAGAGACCCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV7-9*00(629.2) TRBD100(40),TRBD200(35) TRBJ2-5*00(227.1) TRBC2*00(283.5) 273 290 310 0 17 85.0 10 18 36 18 26 40.0;14 21 48 18 25 35.0 22 40 68 27 45
37 18663 105 0.00322296 TGTGCCACCTGGGACAGCACTCTTCGAATTATTATAAGAAACTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV3*00(509.7) nan TRGJ2*00(260.5) TRGC2*00(341.8) 282 299 320 0 17 85.0 nan 17 42 70 22 47 125.0 nan nan nan nan nan nan nan nan nan nan nan
38 18317 111 0.00316321 TGTGCCGCAGGATCTGGGGGTTACCAGAAAGTTACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(616) nan TRAJ13*00(272.7) TRAC*00(264) 264 271 297 0 7 35.0 nan 25 52 83 12 39 135.0 nan nan nan nan nan nan nan nan nan nan nan
39 18128 118 0.00313057 TGTGCCACCTGGGACGGGCTGATACCACTGGTTGGTTCAAGATATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV2*00(500) nan TRGJP1*00(295.5) TRGC1*00(356.5) 282 301 320 0 19 95.0 nan 20 46 80 21 47 130.0 nan nan nan nan nan nan nan nan nan nan nan
40 17946 136 0.00309914 TGTGCAGGGTTCTCATACTCTGGGGCTGGGAGTTACCAACTCACTTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV25*00(558) nan TRAJ28*00(325.9) TRAC*00(257.8) 264 273 293 0 9 45.0 nan 20 55 86 13 48 175.0 nan nan nan nan nan nan nan nan nan nan nan
41 17558 118 0.00303214 TGTCGCCACTCTGGGGCTGGGAGTTACCAACTCACTTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(617.3) nan TRAJ28*00(312) TRAC*00(259.5) 264 267 297 0 3 15.0 nan 23 55 86 7 39 160.0 nan nan nan nan nan nan nan nan nan nan nan
42 17547 155 0.00303024 TGTGCTACAGCTGCAGGCAACAAGCTAACTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(683.5) nan TRAJ17*00(277.5) TRAC*00(261.8) 264 272 297 0 8 40.0 nan 27 52 83 8 33 125.0 nan nan nan nan nan nan nan nan nan nan nan
43 17189 135 0.00296841 TGCGCCAGCAGCCAAGATTTAGTGGCAGCCTACACCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV4-3*00(646.8) nan TRBJ1-2*00(187.9) TRBC1*00(285.2) 270 288 307 0 18 90.0 nan 30 40 68 29 39 50.0 nan nan nan nan nan nan nan nan nan nan nan
44 16467 102 0.00284373 TGTGCCAGTAGTATCGGACAGGGGGCTAGAACCACTGAAGCTTTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV19*00(569.4) TRBD1*00(55) TRBJ1-1*00(219.8) TRBC1*00(280) 270 284 307 0 14 70.0 13 24 36 15 26 55.0 21 40 68 28 48 I25C 83.0 nan nan nan nan nan
45 16050 148 0.00277172 TGCGGCACAGGCGATGCTGGTGGTACTAGCTATGGAAAGCTGACATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV30*00(596.9) nan TRAJ52*00(326) TRAC*00(259.4) 261 271 294 0 10 50.0 nan 23 58 89 13 48 175.0 nan nan nan nan nan nan nan nan nan nan nan
46 15475 89 0.00267242 TGTGCAGGGCCCAAGGATAGCAACTATCAGTTAATCTGG NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV25*00(650.6) nan TRAJ33*00(282.9) TRAC*00(253.5) 264 276 293 0 12 60.0 nan 17 46 77 10 39 ST20A 131.0 nan nan nan nan nan nan nan nan nan nan nan
47 15339 115 0.00264893 TGTGCCGTGAAAGGGGCTAACGACTACAAGCTCAGCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV8-1*00(605.6) nan TRAJ20*00(266.5) TRAC*00(256) 270 281 304 0 11 55.0 nan 23 46 77 16 39 115.0 nan nan nan nan nan nan nan nan nan nan nan
48 15178 141 0.00262113 TGTGCAATGAGAGAGGTATACAACTTCAACAAATTTTACTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV14DV4*00(609.9) nan TRAJ21*00(287.3) TRAC*00(263.5) 273 289 310 0 16 80.0 nan 17 44 75 15 42 135.0 nan nan nan nan nan nan nan nan nan nan nan
49 15139 89 0.00261439 TGTGCCACCTGGGACGGGCGTAGTAGTGATTGGATCAAGACGTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV2*00(529.3) nan TRGJP2*00(289.5) TRGC2*00(346.7) 282 301 320 0 19 95.0 nan 21 46 80 20 45 125.0 nan nan nan nan nan nan nan nan nan nan nan
50 15128 108 0.00261249 TGTGCCGCTACGGGAGGAGGAAACAAACTCACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(641.4) nan TRAJ10*00(287.1) TRAC*00(259.3) 264 271 297 0 7 35.0 nan 26 53 84 9 36 135.0 nan nan nan nan nan nan nan nan nan nan nan
51 14503 86 0.00250456 TGCAGTGCCAGACCTCAGGGGGCCTGGGGCCAACGTCCTGACTTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV20-1*00(580.1) TRBD1*00(40) TRBJ2-6*00(251.1) TRBC2*00(286) 279 291 313 0 12 ST287C 46.0 16 24 36 15 23 40.0 22 45 73 23 46 115.0 nan nan nan nan nan
52 14035 82 0.00242374 TGTGCAGCAAGTCCAACTGACAGCTGGGGGAAATTGCAGTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-1*00(643.9) nan TRAJ24*00(295.6) TRAC*00(251.7) 267 279 300 0 12 60.0 nan 23 52 83 13 42 145.0 nan nan nan nan nan nan nan nan nan nan nan
53 13869 110 0.00239507 TGTGCTACGGTCTTTAATGCTGGTGGTACTAGCTATGGAAAGCTGACATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(584.3) nan TRAJ52*00(336) TRAC*00(261.4) 264 274 297 0 10 50.0 nan 21 58 89 14 51 185.0 nan nan nan nan nan nan nan nan nan nan nan
54 13467 91 0.00232565 TGTGCAGAGAACCCCCCCGGCCAGAAGCTGCTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-2*00(631.8) nan TRAJ16*00(227.6) TRAC*00(267.8) 267 278 300 0 11 55.0 nan 31 49 80 18 36 90.0 nan nan nan nan nan nan nan nan nan nan nan
55 13065 78 0.00225623 TGTGCAATGAGCGCTTACAATGCCAGACTCATGTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-3*00(673.1) nan TRAJ31*00(250.9) TRAC*00(247.6) 264 278 297 0 14 70.0 nan 26 46 77 16 36 100.0 nan nan nan nan nan nan nan nan nan nan nan
56 13023 82 0.00224898 TGTGCAGCAAGTATTAATTCAGGATACAGCACCCTCACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-1*00(659.2) nan TRAJ11*00(287.7) TRAC*00(253.9) 267 281 300 0 14 70.0 nan 19 49 80 12 42 SG21T 136.0 nan nan nan nan nan nan nan nan nan nan nan
57 12655 78 0.00218543 TGGGGAGGAACCTACAAATACATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(651.8) nan TRAJ40*00(262.1) TRAC*00(264.7) 264 266 297 0 2 10.0 nan 28 50 81 5 27 110.0 nan nan nan nan nan nan nan nan nan nan nan
58 12461 68 0.00215192 TGTGCCACCTGGGATGGACATTATAAGAAACTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV4*00(555.4) nan TRGJ2*00(220.9) TRGC2*00(349.9) 282 299 320 0 17 85.0 nan 25 42 70 19 36 85.0 nan nan nan nan nan nan nan nan nan nan nan
59 12435 82 0.00214743 TGTGCTTCAGGGGGAGGAAGCCAAGGAAATCTCATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV21*00(633.2) nan TRAJ42*00(286.9) TRAC*00(256.7) 267 273 299 0 6 30.0 nan 28 55 86 12 39 135.0 nan nan nan nan nan nan nan nan nan nan nan
60 12406 191 0.00214243 TGTGCCAGCAGCTGGACAGGGGGCTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV7-2*00(661.2) TRBD1*00(50) TRBJ2-7*00(217.9) TRBC2*00(288.1) 273 286 310 0 13 65.0 13 23 36 13 23 50.0 23 39 67 23 39 80.0 nan nan nan nan nan
61 12287 79 0.00212187 TGTGCTGCGTGGGATTATACCACTGGTTGGTTCAAGATATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV10*00(522.5) nan TRGJP1*00(300.1) TRGC1*00(343.9) 291 308 329 0 17 85.0 nan 13 46 80 8 42 ST17GI20T 139.0 nan nan nan nan nan nan nan nan nan nan nan
62 12226 94 0.00211134 TGTGCCAGCAGTGAGCAGGGTTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV6-5*00(670.2) TRBD1*00(25) TRBJ2-7*00(212.9) TRBC2*00(287.8) 270 282 307 0 12 60.0 16 21 36 15 20 25.0 24 39 67 21 36 75.0 nan nan nan nan nan
63 12151 84 0.00209839 TGTGCCAGCAGTTACGGTTGGACAGGGACGAGTGGCTACACCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV6-5*00(636.7) TRBD1*00(40) TRBJ1-2*00(202.6) TRBC1*00(292) 270 285 307 0 15 75.0 13 21 36 19 27 40.0 27 40 68 32 45 65.0 nan nan nan nan nan
64 12020 65 0.00207577 TGTGCCATCAGTGAGATGACAGGGGACTCCTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV10-3*00(616.9) TRBD1*00(30) TRBJ2-7*00(235.3) TRBC2*00(280.4) 270 285 307 0 15 75.0 14 20 36 17 23 30.0 16 39 67 23 45 DG19 98.0 nan nan nan nan nan
65 11897 73 0.00205452 TGTGCCACCTGGGACAGGCCTCCCGTGATTGGATCAAGACGTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRGV7*00(532.9) nan TRGJP2*00(266) TRGC2*00(347.5) 281 302 319 0 21 105.0 nan 26 46 80 24 44 100.0 nan nan nan nan nan nan nan nan nan nan nan
66 11541 77 0.00199305 TGTGCTGTCAACTGGGAATATGGAAACAAACTGGTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV41*00(646.3) nan TRAJ47*00(276.1) TRAC*00(256.2) 258 268 290 0 10 50.0 nan 21 46 77 14 39 125.0 nan nan nan nan nan nan nan nan nan nan nan
67 11497 81 0.00198545 TGTGCTACGGACGACGGCGACTACAAGCTCAGCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(699.3) nan TRAJ20*00(248) TRAC*00(257.1) 264 277 297 0 13 65.0 nan 27 46 77 17 36 95.0 nan nan nan nan nan nan nan nan nan nan nan
68 11333 100 0.00195713 TGTGCTCTGAGTGGGTATAACACCGACAAGCTCATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV19*00(668.2) nan TRAJ34*00(272.2) TRAC*00(260.5) 273 286 310 0 13 65.0 nan 23 47 78 15 39 120.0 nan nan nan nan nan nan nan nan nan nan nan
69 10912 68 0.00188442 TGTGCTACGTCCCGTGGAGGTAGCAACTATAAACTGACATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(645.6) nan TRAJ53*00(297.3) TRAC*00(261.7) 264 273 297 0 9 45.0 nan 26 55 86 13 42 145.0 nan nan nan nan nan nan nan nan nan nan nan
70 10911 98 0.00188425 TGTGCTACGGACACCCACTACAAGCTCAGCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(687.4) nan TRAJ20*00(237.8) TRAC*00(267.5) 264 276 297 0 12 60.0 nan 29 46 77 16 33 85.0 nan nan nan nan nan nan nan nan nan nan nan
71 10672 92 0.00184298 TGTGCAGCCCAAACCAGTTCTATTTT NNNNNNNNNNNNNNNNNNNNNNNNNN TRAV23DV6*00(685.4) nan TRAJ49*00(226.8) TRAC*00(261.6) 267 275 300 0 8 40.0 nan 30 45 76 11 26 75.0 nan nan nan nan nan nan nan nan nan nan nan
72 10598 88 0.0018302 TGTGCTACGGGGGGGTACAATAACAATGACATGCGCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(638.9) nan TRAJ43*00(277.1) TRAC*00(260.6) 264 274 297 0 10 50.0 nan 18 43 74 14 39 125.0 nan nan nan nan nan nan nan nan nan nan nan
73 10531 69 0.00181863 TGCATCGTCAGAGTCGCGGAATATGGAAACAAACTGGTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV26-1*00(650.4) nan TRAJ47*00(275.5) TRAC*00(262.6) 261 279 297 0 18 90.0 nan 21 46 77 17 42 125.0 nan nan nan nan nan nan nan nan nan nan nan
74 10438 86 0.00180257 TGTGCTGTGGCCTTGAACTTCAACAAATTTTACTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV21*00(653.6) nan TRAJ21*00(257.2) TRAC*00(259.6) 267 276 299 0 9 45.0 nan 23 44 75 15 36 105.0 nan nan nan nan nan nan nan nan nan nan nan
75 10252 83 0.00177045 TGTGCTACGGACGGAGGAGGTGCTGACGGACTCACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(664.7) nan TRAJ45*00(286.9) TRAC*00(258.2) 264 277 297 0 13 65.0 nan 28 55 86 12 39 135.0 nan nan nan nan nan nan nan nan nan nan nan
76 9998 78 0.00172658 TGCAGTGCTAGAAAGGGCGATATTCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV20-1*00(659.6) TRBD100(25),TRBD200(25) TRBJ2-4*00(197) TRBC2*00(287.3) 279 291 313 0 12 60.0 20 25 36 14 19 25.0;28 33 48 13 18 25.0 30 42 70 21 33
77 9848 80 0.00170068 TGTGCCAGCAGCGTAGGAGCCAGCGAGCAGTTCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV9*00(669.1) nan TRBJ2-1*00(197.6) TRBC2*00(277.9) 270 286 306 0 16 80.0 nan 30 42 70 24 36 60.0 nan nan nan nan nan nan nan nan nan nan nan
78 9712 76 0.00167719 TGTGCTACGGACGCGAACGACTACAAGCTCAGCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(697.5) nan TRAJ20*00(257.5) TRAC*00(264.2) 264 279 297 0 15 75.0 nan 25 46 77 15 36 105.0 nan nan nan nan nan nan nan nan nan nan nan
79 9542 51 0.00164783 TGTGCAGCAAGTATGGGTAACCAGGGAGGAAAGCTTATCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-1*00(640.3) nan TRAJ23*00(276.5) TRAC*00(261.2) 267 281 300 0 14 70.0 nan 27 52 83 17 42 125.0 nan nan nan nan nan nan nan nan nan nan nan
80 9426 82 0.0016278 TGTGTGGTGAACATCCCCGCATCAGGAGGAAGCTACATACCTACATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-1*00(610.5) nan TRAJ6*00(301) TRAC*00(256.6) 261 275 294 0 14 70.0 nan 21 51 82 18 48 150.0 nan nan nan nan nan nan nan nan nan nan nan
81 9385 60 0.00162072 TGTGCCAGCAGTCCCGGAGGGGGTCAGGGTTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV12-300(618.1),TRBV12-400(618.1) TRBD100(31),TRBD200(30) TRBJ2-7*00(212.4) TRBC2*00(279) 273 285 310 0 12 60.0;273 285 310 0 12 60.0 12 21 36 20 29 SA15T 31.0;26 32 48 15 21 30.0
82 9224 71 0.00159292 TGTGCTACGGTTCCAATTCAGGGAGCCCAGAAGCTGGTATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV17*00(644.5) nan TRAJ54*00(292.2) TRAC*00(260.3) 264 274 297 0 10 50.0 nan 21 49 80 14 42 140.0 nan nan nan nan nan nan nan nan nan nan nan
83 9138 41 0.00157807 TGTGCCTTTAACTCAAATTCCGGGTATGCACTCAACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV24*00(634.8) nan TRAJ41*00(301.1) TRAC*00(264.4) 267 277 297 0 10 50.0 nan 21 51 82 9 39 150.0 nan nan nan nan nan nan nan nan nan nan nan
84 8969 59 0.00154888 TGTGCTTATTATAGTGGAGGTAGCAACTATAAACTGACATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV19*00(622) nan TRAJ53*00(311.4) TRAC*00(252.4) 273 279 310 0 6 30.0 nan 23 55 86 10 42 160.0 nan nan nan nan nan nan nan nan nan nan nan
85 8959 56 0.00154715 TGTGCTGGGCTCAAAGCTGCAGGCAACAAGCTAACTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV21*00(632.3) nan TRAJ17*00(296.9) TRAC*00(262.5) 267 274 299 0 7 35.0 nan 23 52 83 10 39 145.0 nan nan nan nan nan nan nan nan nan nan nan
86 8627 68 0.00148982 TGTGCCAGCAGTCGTGGGACAGGGGGCATCGGGGAGCTGTTTTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV7-9*00(558.8) TRBD1*00(55) TRBJ2-2*00(217.9) TRBC2*00(284.8) 273 284 310 0 11 55.0 12 23 36 15 26 55.0 24 43 71 26 45 SC26T 81.0 nan nan nan nan nan
87 8383 51 0.00144768 TGTGCCAGCAGTTCCGCCGGGACAGCCCCCTACAATGAGCAGTTCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV12-300(595.2),TRBV12-400(580.9) TRBD1*00(45) TRBJ2-1*00(236.7) TRBC2*00(283.4) 273 286 310 0 13 65.0;273 286 310 0 13 65.0 10 19 36 16 25 45.0 22 42 70 28 48
88 8238 65 0.00142264 TGTGCTCTGAGTGAGGCATACAACTTCAACAAATTTTACTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV19*00(673.2) nan TRAJ21*00(277.3) TRAC*00(264.3) 273 290 310 0 17 85.0 nan 19 44 75 17 42 125.0 nan nan nan nan nan nan nan nan nan nan nan
89 8199 36 0.00141591 TGTGCCAGCAGTTTAAGGGGGGCGCGGAATCAGCCCCAGCATTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV12-300(606.5),TRBV12-400(606.5) TRBD1*00(40) TRBJ1-5*00(227) TRBC1*00(278.6) 273 288 310 0 15 75.0;273 288 310 0 15 75.0 18 26 36 17 25 40.0 24 42 70 27 45
90 8149 50 0.00140727 TGCGCCAGCAGCCAATTGGGGACAGGGCCCTGGATGAACACTGAAGCTTTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV4-3*00(570) TRBD100(45),TRBD200(40) TRBJ1-1*00(241.4) TRBC1*00(284.4) 270 285 307 0 15 75.0 12 21 36 18 27 45.0;28 36 48 23 31 40.0 19 40 68 33 54
91 7945 112 0.00137204 TGTGCCAGCAGTTTATCAGGGGGGCTCTACGAGCAGTACTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV27*00(678.5) nan TRBJ2-7*00(222.5) TRBC2*00(285.6) 270 287 307 0 17 85.0 nan 16 39 67 20 42 SA18GDC23 84.0 nan nan nan nan nan nan nan nan nan nan nan
92 7902 48 0.00136462 TGTGCCGTGAGGCCTCTAAATACTGGAGGCTTCAAAACTATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(611.5) nan TRAJ9*00(292.4) TRAC*00(261.3) 264 274 297 0 10 50.0 nan 22 50 81 17 45 140.0 nan nan nan nan nan nan nan nan nan nan nan
93 7888 42 0.0013622 TGTGCAGAGACGCGTTACCAGAAAGTTACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-2*00(673.1) nan TRAJ13*00(238.4) TRAC*00(260) 267 277 300 0 10 50.0 nan 32 52 83 13 33 100.0 nan nan nan nan nan nan nan nan nan nan nan
94 7867 52 0.00135857 TGTGCAATGAGAACCGGAACCTACAAATACATCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-3*00(658.8) nan TRAJ40*00(256.6) TRAC*00(260.1) 264 275 297 0 11 55.0 nan 29 50 81 15 36 105.0 nan nan nan nan nan nan nan nan nan nan nan
95 7742 47 0.00133699 TGTGCCGTGTACGACAGCTGGGGGAAATTGCAGTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV12-2*00(635.1) nan TRAJ24*00(271.5) TRAC*00(259.1) 264 276 297 0 12 SA273T 46.0 nan 28 52 83 12 36 120.0 nan nan nan nan nan nan nan nan nan nan nan
96 7667 45 0.00132403 TGTGCTTATAGGAGCGCGTATCGAACCAGTGGCTCTAGGTTGACCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV38-2DV8*00(636.8) nan TRAJ58*00(276.5) TRAC*00(250.6) 273 291 309 0 18 90.0 nan 27 52 83 23 48 125.0 nan nan nan nan nan nan nan nan nan nan nan
97 7603 52 0.00131298 TGCGCCAGCAGCCAAGGGTACAGGGGGAACACTGAAGCTTTCTTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV4-3*00(596.7) TRBD1*00(36) TRBJ1-1*00(231.7) TRBC1*00(288.7) 270 286 307 0 16 80.0 12 22 36 16 26 SG14T 36.0 21 40 68 26 45 95.0 nan nan nan nan nan
98 7469 48 0.00128984 TGTGCCAGCAGTTCCGGGACAACTAACTATGGCTACACCTTC NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRBV28*00(629.7) TRBD100(40),TRBD200(35) TRBJ1-2*00(237.7) TRBC1*00(282.9) 270 283 307 0 13 65.0 10 18 36 13 21 40.0;14 21 48 13 20 35.0 20 40 68 22 42
99 7407 50 0.00127913 TGTGCAGCAATCTGGTGGGGCTTTCAGAAACTTGTATTT NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN TRAV13-1*00(646.9) nan TRAJ8*00(257.2) TRAC*00(267.9) 267 277 300 0 10 50.0 nan 28 49 80 18 39 105.0 nan nan nan nan nan nan nan nan nan nan nan

In order to run the analysis for all samples in the project on Linux we can use GNU Parallel in the following way:

#!/usr/bin/env bash

mkdir -p results

ls /raw/*R1* |
    parallel -j 2 --line-buffer \
    "mixcr analyze qiagen-human-rna-tcr-umi-qiaseq \
    {} \
    {=s:R1:R2:=} \
    {=s:.*/:results/:;s:_R.*::=}"

Under the hood pipeline

Under the hood the command above actually executes the following pipeline:

align

Alignment of raw sequencing reads against reference database of V-, D-, J- and C- gene segments.

mixcr align \
    --species hsa \
    -p generic-amplicon-with-umi \
    -OvParameters.geneFeatureToAlign="VTranscriptWithout5UTRWithP" \
    -OvParameters.parameters.floatingLeftBound=false \
    -OjParameters.parameters.floatingRightBound=false \
    -OсParameters.parameters.floatingRightBound=true \
    --tag-pattern "^(R1:*) \ ^(UMI:N{12})N{12}(R2:*)" \
    --report results/mice_tumor_1.report.txt \
    --json-report results/mice_tumor_1.report.json \
     raw/mice_tumor_1_R1.fastq.gz \
     raw/mice_tumor_1_R2.fastq.gz \
     results/mice_tumor_1.vdjca

Option --report is specified here explicitly.

--species hsa
determines the organism species (hsa for Homo Sapiens). Note we choose Human because mice were humanized.
-p generic-amplicon-with-umi
a preset of MiXCR parameters for amplicon data with UMIs.
-OvParameters.geneFeatureToAlign="VTranscriptWithout5UTRWithP"
Sets a V gene feature to align. Check gene features for more info.
-OvParameters.parameters.floatingLeftBound=false
Results in a global alignment algorithm for V gene left bound. We use it because we don't have any primers covering V gene coding sequence.
-OjParameters.parameters.floatingRightBound=false
Results in a global alignment algorithm for J gene right bound, because reverse primers are located in C-gene region.
-OcParameters.parameters.floatingRightBound=true
Results in a local alignment algorithm for C gene right bound, because reverse primers are located in C-gene region.
--tagPattern "^(R1:*) \ ^(UMI:N{12})N{12}(R2:*)"
tag pattern specifies the location of the UMI barcode according to the library structure.

refineTagsAndSort

Corrects sequencing and PCR errors inside barcode sequences. This step does extremely important job by correcting artificial diversity caused by errors in barcodes. In the considered example project it corrects only sequences of UMIs.

mixcr refineTagsAndSort \
    --report results/mice_tumor_1.report.txt \
    --json-report results/mice_tumor_1.report.json \
    results/mice_tumor_1.vdjca \
    results/mice_tumor_1.refined.vdjca

assemble

Assembles alignments into clonotypes and applies several layers of errors correction:

  • quality-dependent correction for sequencing errors
  • PCR-error correction by clustering
  • UMI-based error correction)

Check mixcr assemble for more information.

mixcr assemble \
    -OassemblingFeatures="CDR3"
    -OseparateByJ=true \
    -OseparateByV=true \
    --report results/mice_tumor_1.report.txt \
    --json-report results/mice_tumor_1.report.json \
    results/mice_tumor_1.refined.vdjca \
    results/mice_tumor_1.clns

Options --report and --json-report are specified here explicitly so that the report files will be appended with assembly report.

-OassemblingFeatures=CDR3
By default qiaseq-human-tcr-cdr3 preset assembles clones by CDR3 sequence.
-separateByJ: true
Separate clones with the same assembling feature, but different V-genes.
-separateByV: true
Separate clones with the same assembling feature, but different J-genes.

export

Exports clonotypes from .clns file into human-readable tables.

mixcr exportClones \
    -c TRA \
    -uniqueTagCount UMI \
    results/mice_tumor_1.clns \
    results/mice_tumor_1.clones.TRA.txt
-с TRB
defines a specific chain to be exported.
-uniqueTagCount UMI
adds a column with the number of UMIs for each clone.

Quality control

Now when we have all files processed lets perform Quality Control. That can be easily done using mixcr exportQc function.

mixcr exportQc align results/*.clns figs/alignQc.pdf

align QC

Judging by the alignment report we can tell some major issues occurred in sample preparation. The main problems are the absence of J hits and reads that don't have TCR sequence at all. We can take the worst samples and realign it in a way that will allow us manually investigate these issues.

By default, MiXCR removed non target reads during alignment. We will realign reads for one of the samples (e.g. mice_tumor_2) and extract not aligned reads (--not-aligned-R1, --not-aligned-R2) to separate FASTQ files for manual inspection. See mixcr align for more details. Parameters -OallowNoCDR3PartAlignments=true and -OallowPartialAlignments=true preserve reads that lack CDR3 or V\J genes.

Bellow is the complete command:

mkdir -p debug

mixcr align \
    --species hsa \
    -p default_4.0 \
    -OvParameters.geneFeatureToAlign="VTranscriptWithout5UTRWithP" \
    -OvParameters.parameters.floatingLeftBound=false \
    -OjParameters.parameters.floatingRightBound=false \
    -OcParameters.parameters.floatingRightBound=true \
    --tag-pattern "^(R1:*) \ ^(UMI:N{12})N{12}(R2:*)" \
    --not-aligned-R1 debug/mice_tumor_2_notAligned_R1.fastq \
    --not-aligned-R2 debug/mice_tumor_2_notAligned_R2.fastq \
    -OallowNoCDR3PartAlignments=true \
    -OallowPartialAlignments=true \
    --report debug/mice_tumor_2.report.txt \
    --json-report debug/mice_tumor_2.report.json \
     raw/mice_tumor_2_R1.fastq.gz \
     raw/mice_tumor_2_R2.fastq.gz \
     debug/mice_tumor_2.vdjca

Now we can export raw alignments using the following command and manually investigate alignments that lack J gene and/or CDR3 region:

mixcr exportAlignmentsPretty
    --skip 1000 \
    --limit 100 \
    debug/mice_tumor_2.vdjca \
    debug/mice_tumor_2_debug.alignments.txt
--skip 1000
Skips first 1000 alignments
--limit 100
exports only 100 alignments

Bellow you can see one of the records from the output file with alignments. We can clearly see that J jene and CDR3 sequence are absent in the original reads.

>>> Read ids: 124



Quality   72222727772277777677777252776
Target0 0 GGCAGGGCTGGGGGTCAGGGAAGAATGAT 28  Score


                                      FR2><CDR2        CDR2><FR3
                _  Q  S  P  E  L  I  M  S  I  Y  S  N  C  D  K  D  A  R  R  F  T  P  Q  L  K  K
    Quality     45265222522225247662622522422222652472225462667462625252524525522262722676662622
    Target1   0 GGCAAAGCCCTGAGCTGATAATGTCTATATACTCCAATTGTGACAAAGATGCTCGAAGGTTTACACCACAGCTCAAAAAA 79   Score
 TRAV7-2*00 181 ggAaaagcccCAagAtgCtGatgtctatCtTctccaatGgtgaAaaGgaAgAAGgCagAttCacaGTTcaCctcaaTaaa 260  222
TRAV7D-2*00 181 ggAaaagcccCAagAtgCtGatgtctatCtTctccaatGgtgaAaaGgaAgAAGgCagAttCacaGTTcaCctcaaTaaa 260  222

                                                                                  V>
                                                                           FR3><CDR3
                 A  S  K  Y  G  S  L  L  I  R  D  S  Q  P  R  D  A  A  T  Y  L  C  A  V  N  I  N
    Quality     27577525762222752262267672725254625577527722252227777245772646266762464522272767
    Target1  80 GCCAGCAAGTATGGTTCGCTGCTCATCAGAGACTCCCAGCCCAGGGATGCAGCCACCTACCTCTGTGCCGTGAACATCAA 159  Score
 TRAV7-2*00 261 gccagcCTgCatACttcCctgcAcatcagagactcccagcccagTgaCTcTgcTCTctacctctgtgc             328  222
TRAV7D-2*00 261 gccagcCTgCatACttcCctgcAcatcagagactcccagcccagTgaCTcTgcTCTctacctctgtgc             328  222

                            CDR3><FR4
                              <J                            FR4>
                F  N  K  F  Y  F  G  S  G  T  K  L  N  V  K  P _
  Quality     62272777722222277777762777777777277722767777727777777767777777777777777777777777
  Target1 160 CTTCAACAAATTTTACTTTGGATCTGGGACCAAACTCAATGTAAAACCAAATATCCAGAACCCTGACCCTGCCGTGTACC 239  Score
TRAJ13*00  43                 tttggaActgggacAaaactcCaAgtCGTTccaa                               76   58

Resulting mice_tumor_2_notAligned_R1.fastq and mice_tumor_2_notAligned_R2.fastq files can be manually inspected with BLAST.

Finally, lets generate and look at chain usage report.

mixcr exportQc chainUsage results/*.clns figs/chainUsage.pdf

chain usage QC

We see that, as expected, samples consist of TRA and TRB chains.

Reports

Finally, MiXCR provides a very convenient way to look at the reports generated at ech step. Every .vdjca, .clns and .clna file holds all the reports for every MiXCR function that has been applied to this sample. E.g. in our case .clns file contains reports for mixcr align and mixcr assemble. To output this report use mixcr exportReports as shown bellow. Note --json parameter will output a JSON-formatted report.

mixcr exportReports \
    results/mice_tumor_1.clns \
    results/mice_tumor_1.report.txt
mixcr exportReports \
    --json \
    results/mice_tumor_1.clns \
    results/mice_tumor_1.report.json
Show report file
============== Align Report ==============
Input file(s): /raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R1.fastq.gz,/raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R2.fastq.gz
Output file(s): results/mice_tumor_1.vdjca
Version: ; built=Mon Sep 26 10:55:18 CEST 2022; rev=8c998df1ab; lib=repseqio.v2.0
Command line arguments: --report results/mice_tumor_1.align.report.txt --json-report results/mice_tumor_1.align.report.json --preset local:qiaseq-human-tcr-cdr3 /raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R1.fastq.gz /raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R2.fastq.gz results/mice_tumor_1.vdjca
Analysis time: 0ns
Total sequencing reads: 8434205
Successfully aligned reads: 6386814 (75.73%)
Chimeras: 21 (0%)
Paired-end alignment conflicts eliminated: 6804 (0.08%)
Alignment failed, no hits (not TCR/IG?): 907328 (10.76%)
Alignment failed because of absence of V hits: 682699 (8.09%)
Alignment failed because of absence of J hits: 397253 (4.71%)
No target with both V and J alignments: 42539 (0.5%)
Alignment failed because of low total score: 17572 (0.21%)
Overlapped: 105392 (1.25%)
Overlapped and aligned: 32073 (0.38%)
Alignment-aided overlaps: 6062 (18.9%)
Overlapped and not aligned: 73319 (0.87%)
No CDR3 parts alignments, percent of successfully aligned: 896 (0.01%)
Partial aligned reads, percent of successfully aligned: 25755 (0.4%)
 chains: 1 (0%)
 non-functional: 0 (0%)
TRA chains: 4345983 (68.05%)
TRA non-functional: 508098 (11.69%)
TRB chains: 1615970 (25.3%)
TRB non-functional: 65867 (4.08%)
TRD chains: 34445 (0.54%)
TRD non-functional: 4804 (13.95%)
TRG chains: 390383 (6.11%)
TRG non-functional: 240407 (61.58%)
IGH chains: 6 (0%)
IGH non-functional: 0 (0%)
IGL chains: 5 (0%)
IGL non-functional: 1 (20%)
Realigned with forced non-floating bound: 16669750 (197.64%)
Realigned with forced non-floating right bound in left read: 5429 (0.06%)
Realigned with forced non-floating left bound in right read: 5429 (0.06%)
============== RefineTagsAndSort Report ==============
Input file(s): results/mice_tumor_1.vdjca
Output file(s): results/mice_tumor_1.refined.vdjca
Version: ; built=Mon Sep 26 10:55:18 CEST 2022; rev=8c998df1ab; lib=repseqio.v2.0
Command line arguments: --report results/mice_tumor_1.refine.report.txt --json-report results/mice_tumor_1.refine.report.json results/mice_tumor_1.vdjca results/mice_tumor_1.refined.vdjca
Analysis time: 0ns
Time spent in correction: 0ns
UMI input diversity: 394970
UMI output diversity: 146618 (37.12%)
UMI input reads: 6386814
UMI output count: 6386809 (100%)
UMI mean reads per tag: 16.17
UMI input core diversity: 394493 (99.88%)
UMI input core reads: 6386335 (99.99%)
UMI directly corrected diversity: 248347 (62.88%)
UMI directly corrected reads: 304865 (4.77%)
UMI diversity filtered by tag quality: 5 (0%)
UMI reads filtered by tag quality: 5 (0%)
UMI diversity filtered by whitelist: 0 (0%)
UMI recursively corrected: 29628
Number of output records: 6028293 (94.39%)
Filter report:
  Number of groups: 146618
  Number of groups accepted: 44866 (30.6%)
  Total records weight: 6386809.0
  Records weight accepted: 6028293.0 (94.39%)
  Operator #0:
    Effective threshold: 14.0
============== Assemble Report ==============
Input file(s): results/mice_tumor_1.refined.vdjca
Output file(s): results/mice_tumor_1.clns
Version: ; built=Mon Sep 26 10:55:18 CEST 2022; rev=8c998df1ab; lib=repseqio.v2.0
Command line arguments: --report results/mice_tumor_1.assemble.report.txt --json-report results/mice_tumor_1.assemble.report.json results/mice_tumor_1.refined.vdjca results/mice_tumor_1.clns
Analysis time: 0ns
Number of input groups: 44866
Number of input alignments: 6028293
Number of output pre-clonotypes: 43546
Number of clonotypes per group:
  0: + 1430 (3.19%) = 1430 (3.19%)
  1: + 43192 (96.41%) = 44622 (99.6%)
  2: + 177 (0.4%) = 44799 (100%)
Number of core alignments: 5825926 (96.64%)
Discarded core alignments: 178462 (3.06%)
Empirically assigned alignments: 852 (0.01%)
Empirical assignment conflicts: 1 (0%)
UMI+VJ-gene empirically assigned alignments: 853 (0.01%)
VJ-gene empirically assigned alignments: 0 (0%)
UMI empirically assigned alignments: 0 (0%)
Number of ambiguous UMIs: 177
Number of ambiguous V-genes: 18
Number of ambiguous UMI+V/J-gene combinations: 18
Unassigned alignments: 196243 (3.26%)
Final clonotype count: 3733
Average number of reads per clonotype: 1551.2
Reads used in clonotypes, percent of total: 5790634 (68.66%)
Reads used in clonotypes before clustering, percent of total: 5826317 (69.08%)
Number of reads used as a core, percent of used: 5824845 (99.97%)
Mapped low quality reads, percent of used: 1472 (0.03%)
Reads clustered in PCR error correction, percent of used: 35683 (0.61%)
Reads pre-clustered due to the similar VJC-lists, percent of used: 755 (0.01%)
Reads dropped due to the lack of a clone sequence, percent of total: 23905 (0.28%)
Reads dropped due to a too short clonal sequence, percent of total: 0 (0%)
Reads dropped due to low quality, percent of total: 0 (0%)
Reads dropped due to failed mapping, percent of total: 461 (0.01%)
Reads dropped with low quality clones, percent of total: 0 (0%)
Clonotypes eliminated by PCR error correction: 260
Clonotypes dropped as low quality: 0
Clonotypes pre-clustered due to the similar VJC-lists: 4
TRA chains: 2169 (58.1%)
TRA non-functional: 354 (16.32%)
TRB chains: 1091 (29.23%)
TRB non-functional: 62 (5.68%)
TRD chains: 47 (1.26%)
TRD non-functional: 8 (17.02%)
TRG chains: 426 (11.41%)
TRG non-functional: 245 (57.51%)
{
  "type": "alignerReport",
  "commandLine": "--report results/mice_tumor_1.align.report.txt --json-report results/mice_tumor_1.align.report.json --preset local:qiaseq-human-tcr-cdr3 /raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R1.fastq.gz /raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R2.fastq.gz results/mice_tumor_1.vdjca",
  "inputFiles": [
    "/raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R1.fastq.gz",
    "/raw/PRJNA674781_QiaSeqTCR/mice_tumor_1_R2.fastq.gz"
  ],
  "outputFiles": [
    "results/mice_tumor_1.vdjca"
  ],
  "version": "; built=Mon Sep 26 10:55:18 CEST 2022; rev=8c998df1ab; lib=repseqio.v2.0",
  "trimmingReport": null,
  "totalReadsProcessed": 8434205,
  "aligned": 6386814,
  "notAligned": 2047391,
  "notAlignedReasons": {
    "NoJHits": 397253,
    "LowTotalScore": 17572,
    "NoBarcode": 0,
    "VAndJOnDifferentTargets": 42539,
    "NoVHits": 682699,
    "NoHits": 907328,
    "NoCDR3Parts": 0
  },
  "chimeras": 21,
  "overlapped": 105392,
  "alignmentAidedOverlaps": 6062,
  "overlappedAligned": 32073,
  "overlappedNotAligned": 73319,
  "pairedEndAlignmentConflicts": 6804,
  "vChimeras": 0,
  "jChimeras": 0,
  "chainUsage": {
    "type": "chainUsage",
    "chimeras": 21,
    "total": 6386814,
    "chains": {
      "": {
        "total": 1,
        "nonFunctional": 0,
        "isOOF": 0,
        "hasStops": 0
      },
      "TRA": {
        "total": 4345983,
        "nonFunctional": 508098,
        "isOOF": 456455,
        "hasStops": 51643
      },
      "TRB": {
        "total": 1615970,
        "nonFunctional": 65867,
        "isOOF": 46247,
        "hasStops": 19620
      },
      "TRD": {
        "total": 34445,
        "nonFunctional": 4804,
        "isOOF": 4036,
        "hasStops": 768
      },
      "TRG": {
        "total": 390383,
        "nonFunctional": 240407,
        "isOOF": 225458,
        "hasStops": 14949
      },
      "IGH": {
        "total": 6,
        "nonFunctional": 0,
        "isOOF": 0,
        "hasStops": 0
      },
      "IGL": {
        "total": 5,
        "nonFunctional": 1,
        "isOOF": 1,
        "hasStops": 0
      }
    }
  },
  "realignedWithForcedNonFloatingBound": 16669750,
  "realignedWithForcedNonFloatingRightBoundInLeftRead": 5429,
  "realignedWithForcedNonFloatingLeftBoundInRightRead": 5429,
  "noCDR3PartsAlignments": 896,
  "partialAlignments": 25755,
  "tagReport": {
    "type": "tagReport"
  }
}
{
  "type": "refineTagsAndSort",
  "commandLine": "--report results/mice_tumor_1.refine.report.txt --json-report results/mice_tumor_1.refine.report.json results/mice_tumor_1.vdjca results/mice_tumor_1.refined.vdjca",
  "inputFiles": [
    "results/mice_tumor_1.vdjca"
  ],
  "outputFiles": [
    "results/mice_tumor_1.refined.vdjca"
  ],
  "version": "; built=Mon Sep 26 10:55:18 CEST 2022; rev=8c998df1ab; lib=repseqio.v2.0",
  "correctionReport": {
    "inputRecords": 6386814,
    "outputRecords": 6028293,
    "steps": [
      {
        "tagName": "UMI",
        "inputGroups": 1,
        "inputDiversity": 394970,
        "inputCount": 6386814,
        "coreDiversity": 394493,
        "coreCount": 6386335,
        "directlyCorrectedDiversity": 248347,
        "directlyCorrectedCount": 304865,
        "filteredDiversity": 5,
        "filteredCount": 5,
        "recursivelyCorrected": 29628,
        "diversityFilteredByWhitelist": 0,
        "outputDiversity": 146618,
        "outputCount": 6386809
      }
    ],
    "filterReport": {
      "type": "filter_groups_report",
      "groupingKeys": [
        "UMI"
      ],
      "numberOfGroups": 146618,
      "numberOfGroupsAccepted": 44866,
      "totalWeight": 6386809,
      "totalWeightAccepted": 6028293,
      "operatorReports": [
        {
          "type": "generic_hist_report",
          "threshold": 14
        }
      ],
      "metricHists": [
        {
          "metric": {
            "type": "group_metric_sum_weight",
            "reportHist": {
              "log": true,
              "binNumber": 0,
              "minBinWidth": 0.2
            }
          },
          "hist": {
            "bins": [
              {
                "from": 1,
                "to": 1.5848931924611136,
                "weight": 48029
              },
              {
                "from": 1.5848931924611136,
                "to": 2.51188643150958,
                "weight": 6986
              },
              {
                "from": 2.51188643150958,
                "to": 3.981071705534973,
                "weight": 8072
              },
              {
                "from": 3.981071705534973,
                "to": 6.309573444801933,
                "weight": 19829
              },
              {
                "from": 6.309573444801933,
                "to": 10,
                "weight": 10844
              },
              {
                "from": 10,
                "to": 15.848931924611142,
                "weight": 10546
              },
              {
                "from": 15.848931924611142,
                "to": 25.11886431509581,
                "weight": 7380
              },
              {
                "from": 25.11886431509581,
                "to": 39.810717055349734,
                "weight": 4391
              },
              {
                "from": 39.810717055349734,
                "to": 63.09573444801933,
                "weight": 5943
              },
              {
                "from": 63.09573444801933,
                "to": 100,
                "weight": 5518
              },
              {
                "from": 100,
                "to": 158.48931924611142,
                "weight": 4183
              },
              {
                "from": 158.48931924611142,
                "to": 251.18864315095823,
                "weight": 5180
              },
              {
                "from": 251.18864315095823,
                "to": 398.1071705534973,
                "weight": 7798
              },
              {
                "from": 398.1071705534973,
                "to": 630.9573444801937,
                "weight": 1894
              },
              {
                "from": 630.9573444801937,
                "to": 1000,
                "weight": 24
              },
              {
                "from": 1000,
                "to": 1584.893192461114,
                "weight": 0
              },
              {
                "from": 1584.893192461114,
                "to": 2511.886431509582,
                "weight": 0
              },
              {
                "from": 2511.886431509582,
                "to": 3981.0717055349733,
                "weight": 0
              },
              {
                "from": 3981.0717055349733,
                "to": 6309.573444801937,
                "weight": 0
              },
              {
                "from": 6309.573444801937,
                "to": 10000,
                "weight": 0
              },
              {
                "from": 10000,
                "to": 15848.93192461114,
                "weight": 0
              },
              {
                "from": 15848.93192461114,
                "to": 25118.864315095823,
                "weight": 0
              },
              {
                "from": 25118.864315095823,
                "to": 39810.71705534977,
                "weight": 1
              }
            ],
            "collectionSpec": {
              "log": true,
              "binNumber": 0,
              "minBinWidth": 0.2
            }
          }
        }
      ]
    }
  }
}
{
  "type": "assemblerReport",
  "commandLine": "--report results/mice_tumor_1.assemble.report.txt --json-report results/mice_tumor_1.assemble.report.json results/mice_tumor_1.refined.vdjca results/mice_tumor_1.clns",
  "inputFiles": [
    "results/mice_tumor_1.refined.vdjca"
  ],
  "outputFiles": [
    "results/mice_tumor_1.clns"
  ],
  "version": "; built=Mon Sep 26 10:55:18 CEST 2022; rev=8c998df1ab; lib=repseqio.v2.0",
  "preCloneAssemblerReport": {
    "type": "preCloneAssemblerReport",
    "inputGroups": 44866,
    "inputAlignments": 6028293,
    "clonotypes": 43546,
    "clonotypesPerGroup": {
      "0": 1430,
      "1": 43192,
      "2": 177
    },
    "coreAlignments": 5825926,
    "discardedCoreAlignments": 178462,
    "empiricallyAssignedAlignments": 852,
    "vjEmpiricallyAssignedAlignments": 0,
    "umiEmpiricallyAssignedAlignments": 0,
    "gatEmpiricallyAssignedAlignments": 853,
    "empiricalAssignmentConflicts": 1,
    "unassignedAlignments": 196243,
    "umiConflicts": 177,
    "gatConflicts": 18,
    "geneConflicts": {
      "Variable": 18
    },
    "coreClonotypesDroppedByTagSuffix": 0,
    "coreAlignmentsDroppedByTagSuffix": 0
  },
  "totalReadsProcessed": 8434205,
  "initialClonesCreated": 3997,
  "readsDroppedNoTargetSequence": 23905,
  "readsDroppedTooShortClonalSequence": 0,
  "readsDroppedLowQuality": 0,
  "coreReads": 5824845,
  "readsDroppedFailedMapping": 461,
  "lowQualityRescued": 1472,
  "clonesClustered": 260,
  "readsClustered": 35683,
  "clones": 3733,
  "clonesDroppedAsLowQuality": 0,
  "clonesPreClustered": 4,
  "readsPreClustered": 755,
  "readsInClones": 5790634,
  "readsInClonesBeforeClustering": 5826317,
  "readsDroppedWithLowQualityClones": 0,
  "clonalChainUsage": {
    "type": "chainUsage",
    "chimeras": 0,
    "total": 3733,
    "chains": {
      "TRA": {
        "total": 2169,
        "nonFunctional": 354,
        "isOOF": 328,
        "hasStops": 26
      },
      "TRB": {
        "total": 1091,
        "nonFunctional": 62,
        "isOOF": 56,
        "hasStops": 6
      },
      "TRD": {
        "total": 47,
        "nonFunctional": 8,
        "isOOF": 8,
        "hasStops": 0
      },
      "TRG": {
        "total": 426,
        "nonFunctional": 245,
        "isOOF": 234,
        "hasStops": 11
      }
    }
  }
}