action.skip

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 TGTGCCGTGAACATCAACTTCAACAAATTTTACTTT 45 nan nan nan nan nan nan nan CAVNINFNKFYF nan :::::::::0:1:14:::::14:-2:36:::
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 TGTGCTTATAGGAGGGGCCAGTTCTATTTT 45 nan nan nan nan nan nan nan CAYRRGQFYF nan :::::::::0:-2:14:::::17:-12:30:::
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 TGTGCAGCAAGAAACGGGCAGGAGAGCACTTACTTTT 45 nan nan nan nan nan nan nan CAARNG_RRALTF nan :::::::::0:-2:11:::::10:-2:37:::
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 nan nan nan nan nan nan TGTGCCAGCAGTCCCGGGACAGGGGGCTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSPGTGGYEQYF nan :::::::::0:-5:12:12:3:-1:26:26:-3:42:::
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 TGTGCTGTTAAGACAACTGACAGCTGGGGGAAATTGCAGTTT 45 nan nan nan nan nan nan nan CAVKTTDSWGKLQF nan :::::::::0:-4:8:::::5:5:42:::
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 nan nan nan nan nan nan TGTGCCAGTAGTATCGTGGCAGGGCCTTTAACCGGGGAGCTGTTTTTT 45 nan nan nan nan nan nan nan CASSIVAGPLTGELFF nan :::::::::0:-3:14:17:-9:2:26:30:-5:48:::
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 TGTGCTCTGAGTGGTAACACCGACAAGCTCATCTTT 45 nan nan nan nan nan nan nan CALSGNTDKLIF nan :::::::::0:-4:13:::::14:-5:36:::
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 TGTGCCAGGCGGGGCTCAGGAACCTACAAATACATCTTT 45 nan nan nan nan nan nan nan CARRGSGTYKYIF nan :::::::::0:-8:5:::::14:-5:39:::
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 TGCGCCAGCAGCCAAGCGGACAGCTCTGGGGCCAACGTCCTGACTTTC 45 nan nan nan nan nan nan nan CASSQADSSGANVLTF nan :::::::::0:-1:16:::::21:2:48:::
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 TGTGCTACGGAACTGAACGACTACAAGCTCAGCTTT 45 nan nan nan nan nan nan nan CATELNDYKLSF nan :::::::::0:-2:11:::::15:-5:36:::
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 nan nan nan nan nan nan TGTGCCAGCAGCCAAGATCCCGACAGGGGAGCTTTCTTT 45 nan nan nan nan nan nan nan CASSQDPDRGAFF nan :::::::::0:2:19:21:-2:-2:29:29:-10:39:::
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 nan nan nan nan nan nan TGCGCCAGCAGCCTGGACAGGGGGAATTCACCCCTCCACTTT 45 nan nan nan nan nan nan nan CASSLDRGNSPLHF nan :::::::::0:-4:13:14👎-1:24:24:-7:42:::
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 TGTGCAGCCGGGCTCGGAAAGCTGATTTTT 45 nan nan nan nan nan nan nan CAAGLGKLIF nan :::::::::0:-2:7:::::18:-20:30:::
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 TGTGCTGTGCTGGAGGAAACCAGTGGCTCTAGGTTGACCTTT 45 nan nan nan nan nan nan nan CAVLEETSGSRLTF nan :::::::::0:0:13:::::15:-5:42:::
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 TGCGGCACAGCCACGGGAGGAGGAAACAAACTCACCTTT 45 nan nan nan nan nan nan nan CGTATGGGNKLTF nan :::::::::0:-3:10:::::11:-5:39:::
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 TGTGCCACCTGGGACGGGCCGGACCACTGGTTGGTTCAAGATATTT 45 nan nan nan nan nan nan nan CATWDGPD_TGWFKIF nan :::::::::0:2:20:::::22:-2:46:::
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 TGTGCCGTGTCATCTGGTTCTGCAAGGCAACTGACCTTT 45 nan nan nan nan nan nan nan CAVSSGSARQLTF nan :::::::::0:-4:9:::::9:-2:39:::
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 TGTGCCTTTATCGATAACTATGGTCAGAATTTTGTCTTT 45 nan nan nan nan nan nan nan CAFIDNYGQNFVF nan :::::::::0:1:11:::::12:-2:39:::
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 TGTGCAGAGAGTATAGGAGGAAGCTACATACCTACATTT 45 nan nan nan nan nan nan nan CAESIGGSYIPTF nan :::::::::0:2:15:::::14:-6:39:::
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 TGTGCTGTGCAGGCCCGAGGAACCTACAAATACATCTTT 45 nan nan nan nan nan nan nan CAVQARGTYKYIF nan :::::::::0:2:15:::::17:-8:39:::
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 TGTGCCAGCAGTTTCGGCAGGGGGTCCTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSFGRGSYEQYF nan :::::::::0:-3:14:::::18:5:42:::
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 TGTGCCAGCAGTTTTCTGCTCGAATCCTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSFLLESYEQYF nan :::::::::0:-3:14:::::24:-1:42:::
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 TGCATCCTGAGGCCCCTATTATAACCAGGGAGGAAAGCTTATCTTC 45 nan nan nan nan nan nan nan CILRPLL*_QGGKLIF nan :::::::::0:-4:11:::::18:-4:46:::
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 nan nan nan nan nan nan TGTGCCAGCAGCTATGAGGGCCACCTTTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSYEGHLYEQYF nan :::::::::0:-4:13:15:-11:2:22:27:-4:42:::
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 TGTGCCGTGAAAAAAAACACCGACAAGCTCATCTTT 45 nan nan nan nan nan nan nan CAVKKNTDKLIF nan :::::::::0:-2:11:::::15:-6:36:::
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||90.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTCCCTTTGCGGGGGATGGAACAGATACGCAGTATTTT 45 nan nan nan nan nan nan nan CASSPFAGDGTDTQYF nan :::::::::0:-5:12:20:-6:-1:25:30:-3: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 TGTGCAGCACCTACGGGAGGAGGAAACAAACTCACCTTT 45 nan nan nan nan nan nan nan CAAPTGGGNKLTF nan :::::::::0:-4:9:::::12:-6:39:::
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 TGTGCCAGCAGTTTTGGAGCAGGCTATGGCTACACCTTC 45 nan nan nan nan nan nan nan CASSFGAGYGYTF nan :::::::::0:-3:14:::::23:-4:39:::
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 TGTGCAATGAGGTTCTACATGGATAGCAACTATCAGTTAATCTGG 45 nan nan nan nan nan nan nan CAMRFYMDSNYQLIW nan :::::::::0:-6:11:::::17:2:45:::
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 TGTGCTGTGAACCCTTATAACCAGGGAGGAAAGCTTATCTTC 45 nan nan nan nan nan nan nan CAVNPYNQGGKLIF nan :::::::::0:-2:10:::::14:-4:42:::
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 TGTGCCTCTCAGGATTCTGGGGGTTACCAGAAAGTTACCTTT 45 nan nan nan nan nan nan nan CASQDSGGYQKVTF nan :::::::::0:-3:7:::::13:-3:42:::
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||100.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTGTACTAGCGGGGGGGGACACCGGGGAGCTGTTTTTT 45 nan nan nan nan nan nan nan CASSVLAGGDTGELFF nan :::::::::0:-2:15:15:-4:0:27:28:-3: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 TGCGGCATCTCACAGAACACCGGTAACCAGTTCTATTTT 45 nan nan nan nan nan nan nan CGISQNTGNQFYF nan :::::::::0:-6:7:::::14:0:39:::
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 nan nan nan nan nan nan TGCAGTGCTAGATGGGGGAATGCGGAGACCCAGTACTTC 45 nan nan nan nan nan nan nan CSARWGNAETQYF nan :::::::::0:-2:12:13:-6:-1:18:24:-5:39:::
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||90.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTGCGGGACTGCAAGAGACCCAGTACTTC 45 nan nan nan nan nan nan nan CASSAGLQETQYF nan :::::::::0:-4:13:13:1:-10:20:21:-2: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 TGTGCAATGAGAGAGGGCCTGGGCACCGACAAGCTCATCTTT 45 nan nan nan nan nan nan nan CAMREGLGTDKLIF nan :::::::::0:2:19:::::23:-8:42:::
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||90.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGCTTAGCTCCGGGACAACAAGAGACCCAGTACTTC 45 nan nan nan nan nan nan nan CASSLAPGQQETQYF nan :::::::::0:0:17:18:2:-6:26:27:-2: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 TGTGCCACCTGGGACAGCACTCTTCGAATTATTATAAGAAACTCTTT 45 nan nan nan nan nan nan nan CATWDSTL_NYYKKLF nan :::::::::0:-1:17:::::22:3:47:::
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 TGTGCCGCAGGATCTGGGGGTTACCAGAAAGTTACCTTT 45 nan nan nan nan nan nan nan CAAGSGGYQKVTF nan :::::::::0:-6:7:::::12:-5:39:::
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 TGTGCCACCTGGGACGGGCTGATACCACTGGTTGGTTCAAGATATTT 45 nan nan nan nan nan nan nan CATWDGLI_TGWFKIF nan :::::::::0:1:19:::::21:0:47:::
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 TGTGCAGGGTTCTCATACTCTGGGGCTGGGAGTTACCAACTCACTTTC 45 nan nan nan nan nan nan nan CAGFSYSGAGSYQLTF nan :::::::::0:0:9:::::13:0:48:::
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 TGTCGCCACTCTGGGGCTGGGAGTTACCAACTCACTTTC 45 nan nan nan nan nan nan nan CRHSGAGSYQLTF nan :::::::::0:-10:3:::::7:-3:39:::
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 TGTGCTACAGCTGCAGGCAACAAGCTAACTTTT 45 nan nan nan nan nan nan nan CATAAGNKLTF nan :::::::::0:-5:8:::::8:-7:33:::
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 TGCGCCAGCAGCCAAGATTTAGTGGCAGCCTACACCTTC 45 nan nan nan nan nan nan nan CASSQDLVAAYTF nan :::::::::0:1:18:::::29:-10:39:::
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 nan nan nan nan nan nan TGTGCCAGTAGTATCGGACAGGGGGCTAGAACCACTGAAGCTTTCTTT 45 nan nan nan nan nan nan nan CASSIGQGARTTEAFF nan :::::::::0:-3:14:15👎0:26:28👎48:::
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 TGCGGCACAGGCGATGCTGGTGGTACTAGCTATGGAAAGCTGACATTT 45 nan nan nan nan nan nan nan CGTGDAGGTSYGKLTF nan :::::::::0:-3:10:::::13:-3:48:::
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 TGTGCAGGGCCCAAGGATAGCAACTATCAGTTAATCTGG 45 nan nan nan nan nan nan nan CAGPKDSNYQLIW nan :::::::::0:3:12:::::10:3:39:::
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 TGTGCCGTGAAAGGGGCTAACGACTACAAGCTCAGCTTT 45 nan nan nan nan nan nan nan CAVKGANDYKLSF nan :::::::::0:-3:11:::::16:-3:39:::
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 TGTGCAATGAGAGAGGTATACAACTTCAACAAATTTTACTTT 45 nan nan nan nan nan nan nan CAMREVYNFNKFYF nan :::::::::0:-1:16:::::15:3:42:::
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 TGTGCCACCTGGGACGGGCGTAGTAGTGATTGGATCAAGACGTTT 45 nan nan nan nan nan nan nan CATWDGRSSDWIKTF nan :::::::::0:1:19:::::20:-1:45:::
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 TGTGCCGCTACGGGAGGAGGAAACAAACTCACCTTT 45 nan nan nan nan nan nan nan CAATGGGNKLTF nan :::::::::0:-6:7:::::9:-6:36:::
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 nan nan nan nan nan nan TGCAGTGCCAGACCTCAGGGGGCCTGGGGCCAACGTCCTGACTTTC 45 nan nan nan nan nan nan nan CSARPQGA_GANVLTF nan :::::::::0:-2:12:15:-4:0:23:23:-2:46:::
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 TGTGCAGCAAGTCCAACTGACAGCTGGGGGAAATTGCAGTTT 45 nan nan nan nan nan nan nan CAASPTDSWGKLQF nan :::::::::0:-1:12:::::13:-3:42:::
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 TGTGCTACGGTCTTTAATGCTGGTGGTACTAGCTATGGAAAGCTGACATTT 45 nan nan nan nan nan nan nan CATVFNAGGTSYGKLTF nan :::::::::0:-3:10:::::14:-1:51:::
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 TGTGCAGAGAACCCCCCCGGCCAGAAGCTGCTCTTT 45 nan nan nan nan nan nan nan CAENPPGQKLLF nan :::::::::0:-2:11:::::18:-11:36:::
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 TGTGCAATGAGCGCTTACAATGCCAGACTCATGTTT 45 nan nan nan nan nan nan nan CAMSAYNARLMF nan :::::::::0:1:14:::::16:-6:36:::
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 TGTGCAGCAAGTATTAATTCAGGATACAGCACCCTCACCTTT 45 nan nan nan nan nan nan nan CAASINSGYSTLTF nan :::::::::0:1:14:::::12:1:42:::
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 TGGGGAGGAACCTACAAATACATCTTT 45 nan nan nan nan nan nan nan WGGTYKYIF nan :::::::::0:-11:2:::::5:-8:27:::
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 TGTGCCACCTGGGATGGACATTATAAGAAACTCTTT 45 nan nan nan nan nan nan nan CATWDGHYKKLF nan :::::::::0:-1:17:::::19:-5:36:::
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 TGTGCTTCAGGGGGAGGAAGCCAAGGAAATCTCATCTTT 45 nan nan nan nan nan nan nan CASGGGSQGNLIF nan :::::::::0:-6:6:::::12:-8:39:::
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 nan nan nan nan nan nan TGTGCCAGCAGCTGGACAGGGGGCTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSWTGGYEQYF nan :::::::::0:-4:13:13👎-1:23:23:-3:39:::
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 TGTGCTGCGTGGGATTATACCACTGGTTGGTTCAAGATATTT 45 nan nan nan nan nan nan nan CAAWDYTTGWFKIF nan :::::::::0:-1:17:::::8:7:42:::
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 nan nan nan nan nan nan TGTGCCAGCAGTGAGCAGGGTTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSEQGYEQYF nan :::::::::0:-5:12:15:-4:-3:20:21:-4:36:::
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 nan nan nan nan nan nan TGTGCCAGCAGTTACGGTTGGACAGGGACGAGTGGCTACACCTTC 45 nan nan nan nan nan nan nan CASSYGWTGTSGYTF nan :::::::::0:-2:15:19👎-3:27:32:-7:45:::
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 nan nan nan nan nan nan TGTGCCATCAGTGAGATGACAGGGGACTCCTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CAISEMTGDSYEQYF nan :::::::::0:-2:15:17:-2:-4:23:23:4:45:::
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 TGTGCCACCTGGGACAGGCCTCCCGTGATTGGATCAAGACGTTT 45 nan nan nan nan nan nan nan CATWDRP_RDWIKTF nan :::::::::0:3:21:::::24:-6:44:::
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 TGTGCTGTCAACTGGGAATATGGAAACAAACTGGTCTTT 45 nan nan nan nan nan nan nan CAVNWEYGNKLVF nan :::::::::0:-2:10:::::14:-1:39:::
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 TGTGCTACGGACGACGGCGACTACAAGCTCAGCTTT 45 nan nan nan nan nan nan nan CATDDGDYKLSF nan :::::::::0:0:13:::::17:-7:36:::
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 TGTGCTCTGAGTGGGTATAACACCGACAAGCTCATCTTT 45 nan nan nan nan nan nan nan CALSGYNTDKLIF nan :::::::::0:-4:13:::::15:-3:39:::
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 TGTGCTACGTCCCGTGGAGGTAGCAACTATAAACTGACATTT 45 nan nan nan nan nan nan nan CATSRGGSNYKLTF nan :::::::::0:-4:9:::::13:-6:42:::
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 TGTGCTACGGACACCCACTACAAGCTCAGCTTT 45 nan nan nan nan nan nan nan CATDTHYKLSF nan :::::::::0:-1:12:::::16:-9:33:::
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 TGTGCAGCCCAAACCAGTTCTATTTT 45 nan nan nan nan nan nan nan CAAQ_QFYF nan :::::::::0:-5:8:::::11:-10:26:::
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 TGTGCTACGGGGGGGTACAATAACAATGACATGCGCTTT 45 nan nan nan nan nan nan nan CATGGYNNNDMRF nan :::::::::0:-3:10:::::14:2:39:::
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 TGCATCGTCAGAGTCGCGGAATATGGAAACAAACTGGTCTTT 45 nan nan nan nan nan nan nan CIVRVAEYGNKLVF nan :::::::::0:2:18:::::17:-1:42:::
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 TGTGCTGTGGCCTTGAACTTCAACAAATTTTACTTT 45 nan nan nan nan nan nan nan CAVALNFNKFYF nan :::::::::0:-3:9:::::15:-3:36:::
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 TGTGCTACGGACGGAGGAGGTGCTGACGGACTCACCTTT 45 nan nan nan nan nan nan nan CATDGGGADGLTF nan :::::::::0:0:13:::::12:-8:39:::
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||60.0 nan nan nan nan nan nan nan nan nan nan nan TGCAGTGCTAGAAAGGGCGATATTCAGTACTTC 45 nan nan nan nan nan nan nan CSARKGDIQYF nan :::::::::0:-2:12:14:-8:1:19:21:-10: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 TGTGCCAGCAGCGTAGGAGCCAGCGAGCAGTTCTTC 45 nan nan nan nan nan nan nan CASSVGASEQFF nan :::::::::0:0:16:::::24:-10:36:::
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 TGTGCTACGGACGCGAACGACTACAAGCTCAGCTTT 45 nan nan nan nan nan nan nan CATDANDYKLSF nan :::::::::0:2:15:::::15:-5:36:::
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 TGTGCAGCAAGTATGGGTAACCAGGGAGGAAAGCTTATCTTC 45 nan nan nan nan nan nan nan CAASMGNQGGKLIF nan :::::::::0:1:14:::::17:-7:42:::
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 TGTGTGGTGAACATCCCCGCATCAGGAGGAAGCTACATACCTACATTT 45 nan nan nan nan nan nan nan CVVNIPASGGSYIPTF nan :::::::::0:1:14:::::18:-1:48:::
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 24|39|67|30|45||75.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTCCCGGAGGGGGTCAGGGTTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSPGGGQGYEQYF nan :::::::::0:-5:12:20:0:-3:29:30:-4:45:::
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 TGTGCTACGGTTCCAATTCAGGGAGCCCAGAAGCTGGTATTT 45 nan nan nan nan nan nan nan CATVPIQGAQKLVF nan :::::::::0:-3:10:::::14:-1:42:::
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 TGTGCCTTTAACTCAAATTCCGGGTATGCACTCAACTTC 45 nan nan nan nan nan nan nan CAFNSNSGYALNF nan :::::::::0:0:10:::::9:-1:39:::
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 TGTGCTTATTATAGTGGAGGTAGCAACTATAAACTGACATTT 45 nan nan nan nan nan nan nan CAYYSGGSNYKLTF nan :::::::::0:-11:6:::::10:-3:42:::
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 TGTGCTGGGCTCAAAGCTGCAGGCAACAAGCTAACTTTT 45 nan nan nan nan nan nan nan CAGLKAAGNKLTF nan :::::::::0:-5:7:::::10:-3:39:::
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 nan nan nan nan nan nan TGTGCCAGCAGTCGTGGGACAGGGGGCATCGGGGAGCTGTTTTTT 45 nan nan nan nan nan nan nan CASSRGTGGIGELFF nan :::::::::0:-6:11:15:0:-1:26:26:-4:45:::
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||100.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTTCCGCCGGGACAGCCCCCTACAATGAGCAGTTCTTC 45 nan nan nan nan nan nan nan CASSSAGTAPYNEQFF nan :::::::::0:-4:13:16:2:-5:25:28:-2: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 TGTGCTCTGAGTGAGGCATACAACTTCAACAAATTTTACTTT 45 nan nan nan nan nan nan nan CALSEAYNFNKFYF nan :::::::::0:0:17:::::17:1:42:::
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.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTTTAAGGGGGGCGCGGAATCAGCCCCAGCATTTT 45 nan nan nan nan nan nan nan CASSLRGARNQPQHF nan :::::::::0:-2:15:17:-6:2:25:27:-4: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||105.0 nan nan nan nan nan nan nan nan nan nan nan TGCGCCAGCAGCCAATTGGGGACAGGGCCCTGGATGAACACTGAAGCTTTCTTT 45 nan nan nan nan nan nan nan CASSQLGTGPWMNTEAFF nan :::::::::0:-2:15:18:0:-3:27:33:1: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 TGTGCCAGCAGTTTATCAGGGGGGCTCTACGAGCAGTACTTC 45 nan nan nan nan nan nan nan CASSLSGGLYEQYF nan :::::::::0:0:17:::::20:4:42:::
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 TGTGCCGTGAGGCCTCTAAATACTGGAGGCTTCAAAACTATCTTT 45 nan nan nan nan nan nan nan CAVRPLNTGGFKTIF nan :::::::::0:-3:10:::::17:-2:45:::
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 TGTGCAGAGACGCGTTACCAGAAAGTTACCTTT 45 nan nan nan nan nan nan nan CAETRYQKVTF nan :::::::::0:-3:10:::::13:-12:33:::
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 TGTGCAATGAGAACCGGAACCTACAAATACATCTTT 45 nan nan nan nan nan nan nan CAMRTGTYKYIF nan :::::::::0:-2:11:::::15:-9:36:::
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 TGTGCCGTGTACGACAGCTGGGGGAAATTGCAGTTT 45 nan nan nan nan nan nan nan CAVYDSWGKLQF nan :::::::::0:-1:12:::::12:-8:36:::
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 TGTGCTTATAGGAGCGCGTATCGAACCAGTGGCTCTAGGTTGACCTTT 45 nan nan nan nan nan nan nan CAYRSAYRTSGSRLTF nan :::::::::0:2:18:::::23:-7:48:::
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 nan nan nan nan nan nan TGCGCCAGCAGCCAAGGGTACAGGGGGAACACTGAAGCTTTCTTT 45 nan nan nan nan nan nan nan CASSQGYRGNTEAFF nan :::::::::0:-1:16:16:0:-2:26:26👎45:::
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||100.0 nan nan nan nan nan nan nan nan nan nan nan TGTGCCAGCAGTTCCGGGACAACTAACTATGGCTACACCTTC 45 nan nan nan nan nan nan nan CASSSGTTNYGYTF nan :::::::::0:-4:13:13:2:-6:21:22:0: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 TGTGCAGCAATCTGGTGGGGCTTTCAGAAACTTGTATTT 45 nan nan nan nan nan nan nan CAAIWWGFQKLVF nan :::::::::0:-3:10:::::18:-8:39:::

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
      }
    }
  }
}