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:
mkdir -p raw
aria2c -c -s 16 -x 16 -k 1M -j 8 -i 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 byCDR3
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
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
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
}
}
}
}