MiLaboratories Human TCR RNA Multiplex kit
This kit allows to obtain TCR alpha and beta repertoires for different types of available RNA material, with high sensitivity and UMI-based accuracy.
Bellow you can see the structure of cDNA library.
The data for this tutorial consists of two samples that represent two replicas of the same biological sample. Total RNA was isolated from PBMC and 25ng were used for each cDNA synthesis. Two cDNA libraries were prepared for each sample (one for TCR alpha and one for TCR betta) according to MiLaboratories Human TCR Repertoire RNA Multiplex Kit protocol. TCR sequencing was performed on an Illumina Miseq sequencer using the 300-cycle Miseq reagent kit(Illumina) with pair-end, 2x150 base pair reads.
All data may be downloaded using the script bellow.
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
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRA_FS115_1_S149_R1_001.fastq.gz
out=raw/Multi_TRA_FS115_1_R1.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRA_FS115_1_S149_R2_001.fastq.gz
out=raw/Multi_TRA_FS115_1_R2.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRA_FS115_2_S150_R1_001.fastq.gz
out=raw/Multi_TRA_FS115_2_R1.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRA_FS115_2_S150_R2_001.fastq.gz
out=raw/Multi_TRA_FS115_2_R2.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRB_FS115_1_S153_R1_001.fastq.gz
out=raw/Multi_TRB_FS115_1_R1.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRB_FS115_1_S153_R2_001.fastq.gz
out=raw/Multi_TRB_FS115_1_R2.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRB_FS115_2_S154_R1_001.fastq.gz
out=raw/Multi_TRB_FS115_2_R1.fastq.gz
https://cdn.milaboratory.com/data/milab-tcr-mulriplex/Multi_TRB_FS115_2_S154_R2_001.fastq.gz
out=raw/Multi_TRB_FS115_2_R2.fastq.gz
Upstream analysis
One-line solution
MiXCR has a dedicated preset for this protocol, thus analysing the data is as easy as:
mixcr analyze milab-human-rna-tcr-umi-multiplex \
raw/Multi_TRA_FS115_1_R1.fastq.gz \
raw/Multi_TRA_FS115_1_R2.fastq.gz \
results/Multi_TRA_FS115_1
Running the command above will generate the following files:
> ls result/
# human-readable reports
Multi_TRA_FS115_1.report
# raw alignments (highly compressed binary file)
Multi_TRA_FS115_1.vdjca
# alignments with refined UMI barcode sequences
Multi_TRA_FS115_1.refined.vdjca
# TRA, TRB CDR3 clonotypes (highly compressed binary file)
Multi_TRA_FS115_1.clns
# TRA,TRB CDR3 clonotypes exported in tab-delimited txt
Multi_TRA_FS115_1.clonotypes.TRA.tsv
Multi_TRA_FS115_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 | 125484 | 1724 | 0.0456243 | TGTGCGGAGCTAGCGGGGATATCTAACTTTGGAAATGAGAAATTAACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(277.6) | nan | TRAJ48*00(315.7) | TRAC*00(64.5) | 264 | 269 | 297 | 0 | 5 | 25.0 | nan | 19 | 52 | 83 | 18 | 51 | 165.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
1 | 83058 | 1140 | 0.0301988 | TGTGCCGTGACCGTCTGGTGGGGAGCCCAGAAGCTGGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(292.3) | nan | TRAJ54*00(264.6) | TRAC*00(64.5) | 264 | 274 | 297 | 0 | 10 | 50.0 | nan | 27 | 49 | 80 | 20 | 42 | 110.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
2 | 60550 | 596 | 0.0220152 | TGTGCTGTGCTGTTTATGAACAGAGATGACAAGATCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(253.2) | nan | TRAJ30*00(303.9) | TRAC*00(70.7) | 270 | 279 | 306 | 0 | 9 | 45.0 | nan | 13 | 46 | 77 | 9 | 42 | SC18T | 151.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
3 | 47245 | 1101 | 0.0171777 | TGTGCAGGGAAAAGTGGCAACACAGGCAAACTAATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV25*00(273.2) | nan | TRAJ37*00(278.7) | TRAC*00(65) | 264 | 273 | 293 | 0 | 9 | 45.0 | nan | 26 | 51 | 82 | 14 | 39 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
4 | 43079 | 474 | 0.015663 | TGTGCTGGGAAAACTGGGGCAAACAACCTCTTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(244.6) | nan | TRAJ36*00(283.4) | TRAC*00(89.5) | 270 | 280 | 306 | 0 | 10 | ST277G | 36.0 | nan | 22 | 48 | 79 | 10 | 36 | 130.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
5 | 28006 | 1059 | 0.0101826 | TGTGCTGTGAGTGATCGTTCAGGAGGAGGTGCTGACGGACTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-6*00(273.1) | nan | TRAJ45*00(308.5) | TRAC*00(64.6) | 270 | 286 | 304 | 0 | 16 | 80.0 | nan | 24 | 55 | 86 | 17 | 48 | 155.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
6 | 24675 | 493 | 0.0089715 | TGTGCAGCAATCAACACCGGTAACCAGTTCTATTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV23DV6*00(262.1) | nan | TRAJ49*00(273.3) | TRAC*00(83.1) | 267 | 280 | 300 | 0 | 13 | SG277T | 51.0 | nan | 21 | 45 | 76 | 12 | 36 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
7 | 24629 | 877 | 0.00895478 | TGTGCTGTGAGTGTCCCTTCAGGAGGAAGCTACATACCTACATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-6*00(258.9) | nan | TRAJ6*00(288.1) | TRAC*00(64.9) | 270 | 283 | 304 | 0 | 13 | 65.0 | nan | 24 | 51 | 82 | 18 | 45 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
8 | 24116 | 955 | 0.00876826 | TGTGCTGTCGACACCGACAAGCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV2*00(265.1) | nan | TRAJ34*00(253.2) | TRAC*00(101.7) | 249 | 260 | 283 | 0 | 11 | SG257C | 41.0 | nan | 27 | 47 | 78 | 10 | 30 | 100.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
9 | 21145 | 366 | 0.00768804 | TGTGTGGTGAAATTTGGGTCTGGTTCTGCAAGGCAACTGACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-1*00(296.4) | nan | TRAJ22*00(292.8) | TRAC*00(64.6) | 261 | 272 | 294 | 0 | 11 | 55.0 | nan | 16 | 52 | 83 | 8 | 45 | SC23GST24GI25G | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
10 | 18878 | 253 | 0.00686379 | TGTGCCGTGAAGGAGGATAGCAACTATCAGTTAATCTGG | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(298.1) | nan | TRAJ33*00(277.9) | TRAC*00(64.5) | 264 | 275 | 297 | 0 | 11 | 55.0 | nan | 21 | 46 | 77 | 14 | 39 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
11 | 16805 | 275 | 0.00611008 | TGTGCCGTAGGGGGGGGGCCAAGCAGTGCTTCCAAGATAATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(283.7) | nan | TRAJ3*00(273.5) | TRAC*00(64.6) | 264 | 272 | 297 | 0 | 8 | 40.0 | nan | 27 | 51 | 82 | 21 | 45 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
12 | 16598 | 457 | 0.00603481 | TGTGCTGTGGGTGCACTCTAACTTTGGAAATGAGAAATTAACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-3*00(271.1) | nan | TRAJ48*00(303.3) | TRAC*00(64.6) | 270 | 284 | 304 | 0 | 14 | 70.0 | nan | 22 | 52 | 83 | 16 | 46 | 150.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
13 | 16490 | 233 | 0.00599555 | TGTGCCCCTCGAGAGAGCTATGGTCAGAATTTTGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV24*00(241.9) | nan | TRAJ26*00(263.8) | TRAC*00(69.5) | 267 | 273 | 297 | 0 | 6 | 30.0 | nan | 27 | 49 | 80 | 17 | 39 | 110.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
14 | 14361 | 77 | 0.00522147 | TGCCTCGTGGGGGGGCGCGGTGGTGCTACAAACAAGCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(262.7) | nan | TRAJ32*00(288.2) | TRAC*00(66.2) | 261 | 272 | 297 | 0 | 11 | 55.0 | nan | 28 | 55 | 86 | 18 | 45 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
15 | 13243 | 361 | 0.00481498 | TGTGCTAAACTTACCAGTGGCTCTAGGTTGACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV17*00(260.3) | nan | TRAJ58*00(272.8) | TRAC*00(74.2) | 264 | 271 | 297 | 0 | 7 | 35.0 | nan | 28 | 52 | 83 | 12 | 36 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
16 | 12894 | 151 | 0.00468809 | TGTGCTGTGAGTGGATGGCCAGAAGCTGCTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-4*00(270.8) | nan | TRAJ16*00(257.9) | TRAC*00(101.9) | 270 | 283 | 304 | 0 | 13 | 65.0 | nan | 28 | 49 | 80 | 13 | 34 | 105.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
17 | 10482 | 115 | 0.00381112 | TGTGCTGTGAGAGATACTTGGGGAAGCCAAGGAAATCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(277) | nan | TRAJ42*00(276.8) | TRAC*00(64.7) | 270 | 284 | 306 | 0 | 14 | 70.0 | nan | 27 | 55 | 86 | 18 | 45 | DA30 | 123.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
18 | 9959 | 459 | 0.00362096 | TGTGCAGCCATTTCTGGTGGCTACAATAAGCTGATTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV13-1*00(261.9) | nan | TRAJ4*00(298.3) | TRAC*00(65.5) | 267 | 275 | 300 | 0 | 8 | 40.0 | nan | 23 | 52 | 83 | 10 | 39 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
19 | 9450 | 159 | 0.0034359 | TGTGCTCTGACCGTTGGATTCAATGACATGCGCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV9-2*00(231.9) | nan | TRAJ43*00(234.5) | TRAC*00(109.7) | 267 | 277 | 301 | 0 | 10 | 50.0 | nan | 27 | 43 | 74 | 20 | 36 | 80.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
20 | 8077 | 115 | 0.00293669 | TGTGCAGGAAAGGAAACGGGGGAAACCAGTGGCTCTAGGTTGACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV27*00(267.4) | nan | TRAJ58*00(288) | TRAC*00(64.7) | 264 | 273 | 294 | 0 | 9 | 45.0 | nan | 25 | 52 | 83 | 21 | 48 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
21 | 8055 | 112 | 0.00292869 | TGTGCCGGGCAGGGCGGATCTGAAAAGCTGGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(279.6) | nan | TRAJ57*00(288.2) | TRAC*00(64.5) | 264 | 271 | 297 | 0 | 7 | 35.0 | nan | 25 | 52 | 83 | 9 | 36 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
22 | 6582 | 162 | 0.00239313 | TGTGCTGGCCCCCCCCAGAGAGGAGGAAACAAACTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV35*00(265.4) | nan | TRAJ10*00(268.1) | TRAC*00(64.7) | 264 | 272 | 296 | 0 | 8 | 40.0 | nan | 30 | 53 | 84 | 19 | 42 | 115.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
23 | 6360 | 106 | 0.00231241 | TGTGTGGTGAACTAACCGGGGAGGTGGAGGCTTCAAAACTATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-1*00(302.1) | nan | TRAJ9*00(263.3) | TRAC*00(64.4) | 261 | 273 | 294 | 0 | 12 | 60.0 | nan | 28 | 50 | 81 | 24 | 46 | 110.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
24 | 6312 | 176 | 0.00229496 | TGTGCTGTGTCGTATTCAGGAGGAGGTGCTGACGGACTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-3*00(246.9) | nan | TRAJ45*00(322.8) | TRAC*00(64.7) | 270 | 279 | 304 | 0 | 9 | 45.0 | nan | 21 | 55 | 86 | 11 | 45 | 170.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
25 | 6036 | 274 | 0.00219461 | TGTGCAGCAAGCCGGACCCAGGGAGCCCAGAAGCTGGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV13-1*00(275.3) | nan | TRAJ54*00(273.4) | TRAC*00(64.7) | 267 | 278 | 300 | 0 | 11 | 55.0 | nan | 25 | 49 | 80 | 18 | 42 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
26 | 5863 | 22 | 0.00213171 | TGCCTCGTCTACAACTTCAACAAATTTTACTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(247.4) | nan | TRAJ21*00(273.8) | TRAC*00(103.9) | 261 | 269 | 297 | 0 | 8 | 40.0 | nan | 17 | 44 | 75 | 6 | 33 | SA19C | 121.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
27 | 5472 | 72 | 0.00198955 | TGTGCTGTCGGTTCTGGGGGTTACCAGAAAGTTACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(247.5) | nan | TRAJ13*00(292.3) | TRAC*00(76.9) | 270 | 278 | 306 | 0 | 8 | 40.0 | nan | 24 | 52 | 83 | 11 | 39 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
28 | 4780 | 56 | 0.00173794 | TGTGCTGTGAGAGACGCGCAAGGAGGAGGTGCTGACGGACTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(280.8) | nan | TRAJ45*00(293.6) | TRAC*00(64.5) | 270 | 285 | 306 | 0 | 15 | 75.0 | nan | 27 | 55 | 86 | 20 | 48 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
29 | 4778 | 116 | 0.00173722 | TGTGCAATGATCGGGAACACAGGCTTTCAGAAACTTGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-3*00(277.1) | nan | TRAJ8*00(293.5) | TRAC*00(64.7) | 264 | 277 | 297 | 0 | 13 | SG274T | 51.0 | nan | 21 | 49 | 80 | 14 | 42 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
30 | 4648 | 209 | 0.00168995 | TGTGCAGCAAGTATGTTCCCAGGAAGCCAAGGAAATCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV13-1*00(291.4) | nan | TRAJ42*00(278.9) | TRAC*00(64.7) | 267 | 281 | 300 | 0 | 14 | 70.0 | nan | 30 | 55 | 86 | 20 | 45 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
31 | 4550 | 81 | 0.00165432 | TGTGTGGTGAACCCTAGGGACCAGTTCTATTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-1*00(303) | nan | TRAJ49*00(223.4) | TRAC*00(65.5) | 261 | 273 | 294 | 0 | 12 | 60.0 | nan | 31 | 45 | 76 | 19 | 33 | 70.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
32 | 4520 | 64 | 0.00164341 | TGCATCCTGGACGGCAAGCAGGGCGGATCTGAAAAGCTGGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-2*00(267) | nan | TRAJ57*00(287.9) | TRAC*00(64.9) | 261 | 270 | 296 | 0 | 9 | 45.0 | nan | 25 | 52 | 83 | 18 | 45 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
33 | 4399 | 52 | 0.00159942 | TGTGCTGTGAGAGCGTATTTAGCAGGAAACACACCTCTTGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(271.9) | nan | TRAJ29*00(270.3) | TRAC*00(64.5) | 270 | 283 | 306 | 0 | 13 | 65.0 | nan | 26 | 49 | 80 | 22 | 45 | 115.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
34 | 4388 | 196 | 0.00159542 | TGTGCAGCAAGTTACCTCTATGGGAACAACAGACTCGCTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV13-1*00(284.2) | nan | TRAJ7*00(278.3) | TRAC*00(64.7) | 267 | 283 | 300 | 0 | 15 | DA279 | 63.0 | nan | 23 | 48 | 79 | 17 | 42 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
35 | 4359 | 37 | 0.00158487 | TGCATCGTCAGAGTCGCGGATCGGCAGTGCTTCCAAGATAATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-1*00(307.1) | nan | TRAJ3*00(271.1) | TRAC*00(64.8) | 261 | 279 | 297 | 0 | 18 | 90.0 | nan | 28 | 51 | 82 | 23 | 46 | 115.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
36 | 4355 | 71 | 0.00158342 | TGTGCCGGGGGGGCGAATTATGGTGGTGCTACAAACAAGCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(282.5) | nan | TRAJ32*00(323) | TRAC*00(64.7) | 264 | 271 | 297 | 0 | 7 | 35.0 | nan | 21 | 55 | 86 | 14 | 48 | 170.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
37 | 4320 | 34 | 0.00157069 | TGTGCAATGGAGGGGGATAGCAGCTATAAATTGATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV14DV4*00(257.5) | nan | TRAJ12*00(278.6) | TRAC*00(70.5) | 273 | 290 | 310 | 0 | 16 | DA282SA286G | 54.0 | nan | 24 | 49 | 80 | 14 | 39 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
38 | 3833 | 39 | 0.00139363 | TGTGCTGTGAGACCCATTGATGGAGGAAGCCAAGGAAATCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(268.2) | nan | TRAJ42*00(302.3) | TRAC*00(64.5) | 270 | 282 | 306 | 0 | 12 | 60.0 | nan | 23 | 55 | 86 | 15 | 48 | I26G | 148.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
39 | 3491 | 28 | 0.00126928 | TGCATCGTCAGAGTCTTATACTCTGGGGCTGGGAGTTACCAACTCACTTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-1*00(297.7) | nan | TRAJ28*00(323) | TRAC*00(64.8) | 261 | 276 | 297 | 0 | 15 | 75.0 | nan | 11 | 55 | 86 | 8 | 51 | SA17CSG19TDC20 | 175.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
40 | 3321 | 41 | 0.00120747 | TGTGCCGTCTTTTCGGGCACTGCCAGTAAACTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(282.3) | nan | TRAJ44*00(274.1) | TRAC*00(64.2) | 264 | 272 | 297 | 0 | 8 | 40.0 | nan | 28 | 52 | 83 | 15 | 39 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
41 | 3298 | 12 | 0.00119911 | TGCCTCGTGGGTGTTCTGGGGGTTACCAGAAAGTTACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(272.6) | nan | TRAJ13*00(292.8) | TRAC*00(73.1) | 261 | 274 | 297 | 0 | 13 | 65.0 | nan | 24 | 52 | 83 | 13 | 41 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
42 | 3254 | 29 | 0.00118311 | TGTGCAATGAGAGAGAACACGGGCAGGAGAGCACTTACTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV14DV4*00(276.6) | nan | TRAJ5*00(283.4) | TRAC*00(65.6) | 273 | 288 | 310 | 0 | 15 | 75.0 | nan | 23 | 49 | 80 | 16 | 42 | 130.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
43 | 3099 | 83 | 0.00112676 | TGTGTGGTGAGCGACAGAGGCTCAACCCTGGGGAGGCTATACTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV10*00(274) | nan | TRAJ18*00(322.8) | TRAC*00(64.8) | 267 | 280 | 300 | 0 | 13 | 65.0 | nan | 21 | 55 | 86 | 11 | 45 | 170.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
44 | 2946 | 76 | 0.00107113 | TGTGCAGAGAATATTGCCCAGGCAGGAACTGCTCTGATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV13-2*00(301.4) | nan | TRAJ15*00(278.4) | TRAC*00(64.6) | 267 | 281 | 300 | 0 | 14 | 70.0 | nan | 24 | 49 | 80 | 17 | 42 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
45 | 2868 | 57 | 0.00104277 | TGTGCAGCAAGCCTTAGGACAGGAGGAAGCTACATACCTACATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV23DV6*00(269.4) | nan | TRAJ6*00(283.4) | TRAC*00(64.1) | 267 | 279 | 300 | 0 | 12 | 60.0 | nan | 25 | 51 | 82 | 19 | 45 | 130.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
46 | 2803 | 11 | 0.00101913 | TGCCTCGTGGGTGACATAGCAGCTATAAATTGATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(290.4) | nan | TRAJ12*00(268.7) | TRAC*00(81.9) | 261 | 278 | 297 | 0 | 17 | 85.0 | nan | 26 | 49 | 80 | 15 | 38 | 115.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
47 | 2774 | 36 | 0.00100859 | TGTGCCGTGAATCGCTCTGACATGCGCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(299.9) | nan | TRAJ43*00(219) | TRAC*00(74) | 264 | 275 | 297 | 0 | 11 | 55.0 | nan | 30 | 43 | 74 | 17 | 30 | 65.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
48 | 2540 | 10 | 0.00092351 | TGCCTCGTGGGTCCTCTTTATAACCAGGGAGGAAAGCTTATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(266.4) | nan | TRAJ23*00(298.1) | TRAC*00(66.3) | 261 | 273 | 297 | 0 | 12 | 60.0 | nan | 23 | 52 | 83 | 16 | 45 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
49 | 2526 | 24 | 0.00091842 | TGTGCTGTGAGAGACACCGCCCGGGTCGTGGAGGAGAGCACTTACTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(288) | nan | TRAJ5*00(243.7) | TRAC*00(64.5) | 270 | 286 | 306 | 0 | 16 | 80.0 | nan | 31 | 49 | 80 | 31 | 49 | 90.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
50 | 2499 | 12 | 0.000908603 | TGTGCCGTCCGAGGCTCCCGGAATAACAATGCCAGACTCATGTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV39*00(277.6) | nan | TRAJ31*00(303.8) | TRAC*00(62.4) | 264 | 272 | 297 | 0 | 8 | 40.0 | nan | 18 | 46 | 77 | 17 | 45 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
51 | 2475 | 9 | 0.000899877 | TGCCTCGTGGGTCCACTTTATAACCAGGGAGGAAAGCTTATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(266.2) | nan | TRAJ23*00(299) | TRAC*00(66.4) | 261 | 273 | 297 | 0 | 12 | 60.0 | nan | 23 | 52 | 83 | 16 | 45 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
52 | 2380 | 31 | 0.000865337 | TGTGCCGTGAATTCAGGATACAGCACCCTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(301.6) | nan | TRAJ11*00(292.4) | TRAC*00(64.7) | 264 | 275 | 297 | 0 | 11 | 55.0 | nan | 20 | 49 | 80 | 7 | 36 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
53 | 2335 | 46 | 0.000848975 | TGTGCTCTAGACGACGACTACAAGCTCAGCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV6*00(320.3) | nan | TRAJ20*00(252.4) | TRAC*00(65.1) | 267 | 279 | 300 | 0 | 12 | 60.0 | nan | 26 | 46 | 77 | 13 | 33 | 100.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
54 | 2309 | 81 | 0.000839522 | TGTGCTTCTCCTGGGGCTGGGAGTTACCAACTCACTTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-2*00(338.2) | nan | TRAJ28*00(298.1) | TRAC*00(64.7) | 261 | 267 | 295 | 0 | 6 | 30.0 | nan | 26 | 55 | 86 | 10 | 39 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
55 | 2237 | 69 | 0.000813344 | TGTGCTGTGAGATGGGCGGGAGGAAGCTACATACCTACATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-2*00(361.4) | nan | TRAJ6*00(272.7) | TRAC*00(64.6) | 261 | 273 | 295 | 0 | 12 | 60.0 | nan | 27 | 51 | 82 | 18 | 42 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
56 | 2188 | 220 | 0.000795528 | TGTGCCGTGAGACGAGGTGGCAACACAGGCAAACTAATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-1*00(291.8) | nan | TRAJ37*00(278.2) | TRAC*00(64.6) | 270 | 280 | 304 | 0 | 10 | 50.0 | nan | 26 | 51 | 82 | 17 | 42 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
57 | 1983 | 40 | 0.000720993 | TGTGCTCTACGGGCGTACAGCAGTGCTTCCAAGATAATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV6*00(304) | nan | TRAJ3*00(301) | TRAC*00(64.8) | 267 | 276 | 300 | 0 | 9 | 45.0 | nan | 19 | 51 | 82 | 9 | 42 | I23C | 148.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
58 | 1961 | 20 | 0.000712994 | TGTGCTGTGTCCCAGGGGGATGGCCAGAAGCTGCTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-4*00(247.8) | nan | TRAJ16*00(258.8) | TRAC*00(78.5) | 270 | 279 | 304 | 0 | 9 | 45.0 | nan | 28 | 49 | 80 | 18 | 39 | 105.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
59 | 1961 | 64 | 0.000712994 | TGTGCTGTGAGAGATAGCAACTATCAGTTAATCTGG | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-2*00(377.3) | nan | TRAJ33*00(272.5) | TRAC*00(64.7) | 261 | 276 | 295 | 0 | 15 | 75.0 | nan | 22 | 46 | 77 | 12 | 36 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
60 | 1920 | 21 | 0.000698087 | TGTGCCCAGGGGACTGGAGGCTTCAAAACTATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(275.1) | nan | TRAJ9*00(272.8) | TRAC*00(64.7) | 264 | 270 | 297 | 0 | 6 | 30.0 | nan | 26 | 50 | 81 | 12 | 36 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
61 | 1918 | 25 | 0.000697359 | TGTGCCGTGAAAAGGCCCACAGGCTTTCAGAAACTTGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(297.5) | nan | TRAJ8*00(278.9) | TRAC*00(64.5) | 264 | 275 | 297 | 0 | 11 | 55.0 | nan | 24 | 49 | 80 | 17 | 42 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
62 | 1900 | 35 | 0.000690815 | GTATTTCCAAGGGGTTCAGGGAGCCCAGAAGCTGGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV27*00(225.3) | nan | TRAJ54*00(283.3) | TRAC*00(65.3) | 264 | 264 | 294 | 0 | 0 | 0.0 | nan | 23 | 49 | 80 | 14 | 40 | 130.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
63 | 1853 | 19 | 0.000673726 | TGTGCTGTGAGTGAAGGTGCAGGCAACATGCTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-4*00(270.8) | nan | TRAJ39*00(264.3) | TRAC*00(79) | 270 | 284 | 304 | 0 | 14 | 70.0 | nan | 30 | 52 | 83 | 17 | 39 | 110.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
64 | 1822 | 32 | 0.000662455 | TGTGCTCTGAGTGAAGAAGGTGGTACTAGCTATGGAAAGCTGACATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV9-2*00(248.4) | nan | TRAJ52*00(303.7) | TRAC*00(65.3) | 267 | 281 | 301 | 0 | 14 | 70.0 | nan | 28 | 58 | 89 | 18 | 48 | 150.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
65 | 1783 | 16 | 0.000648275 | TGTGCTGTGAGAGACGGCGGGGGAAACAAACTGGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(282) | nan | TRAJ47*00(243.6) | TRAC*00(77.1) | 270 | 285 | 306 | 0 | 15 | 75.0 | nan | 28 | 46 | 77 | 21 | 39 | 90.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
66 | 1748 | 51 | 0.00063555 | TGTGCAGAGACCGCCGTCAATACTGGAGGCTTCAAAACTATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV13-2*00(279.6) | nan | TRAJ9*00(288.4) | TRAC*00(64.7) | 267 | 277 | 300 | 0 | 10 | 50.0 | nan | 23 | 50 | 81 | 18 | 45 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
67 | 1742 | 23 | 0.000633368 | TGTGTTGTGAGTGAATGGCAGGCAGGAACTGCTCTGATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-2*00(274.9) | nan | TRAJ15*00(272.6) | TRAC*00(69.8) | 270 | 284 | 304 | 0 | 14 | 70.0 | nan | 25 | 49 | 80 | 18 | 42 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
68 | 1673 | 15 | 0.000608281 | TGTGCAATGAGAGAGGGCGGATTCACGGGAGGAGGAAACAAACTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV14DV4*00(296.2) | nan | TRAJ10*00(297.1) | TRAC*00(64.5) | 273 | 291 | 310 | 0 | 18 | 90.0 | nan | 24 | 53 | 84 | 22 | 51 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
69 | 1610 | 22 | 0.000585375 | TGCATCCTGAGAGACCGGGGCGGGGCAAACAACCTCTTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-2*00(296.7) | nan | TRAJ36*00(257.9) | TRAC*00(64.8) | 261 | 276 | 296 | 0 | 15 | 75.0 | nan | 27 | 48 | 79 | 21 | 42 | 105.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
70 | 1557 | 101 | 0.000566105 | TGTGCTCTGAGTGAGGAGAAGAGGGATAGCAGCTATAAATTGATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV19*00(297.2) | nan | TRAJ12*00(279.6) | TRAC*00(64.8) | 273 | 289 | 310 | 0 | 16 | 80.0 | nan | 21 | 49 | 80 | 20 | 48 | ST23G | 126.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
71 | 1535 | 110 | 0.000558106 | TGTGCTCTGAGTGAGGCGTGGGGAGCCCAGAAGCTGGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV19*00(309.2) | nan | TRAJ54*00(263.4) | TRAC*00(64.9) | 273 | 291 | 310 | 0 | 18 | 90.0 | nan | 27 | 49 | 80 | 20 | 42 | 110.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
72 | 1535 | 9 | 0.000558106 | TGCATCGTCAGACCTCAGGAATATGGAAACAAACTGGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-1*00(278.1) | nan | TRAJ47*00(278) | TRAC*00(64.8) | 261 | 273 | 297 | 0 | 12 | 60.0 | nan | 21 | 46 | 77 | 17 | 42 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
73 | 1500 | 23 | 0.00054538 | TGTGCCGTGAACAATCCAAACTATGGTCAGAATTTTGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(307.5) | nan | TRAJ26*00(272.5) | TRAC*00(64.6) | 264 | 277 | 297 | 0 | 13 | 65.0 | nan | 25 | 49 | 80 | 18 | 42 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
74 | 1494 | 8 | 0.000543199 | TGCCTCTGCGAGGGAGCAGGCAAATCAACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(243.4) | nan | TRAJ27*00(246.2) | TRAC*00(104.2) | 261 | 267 | 297 | 0 | 6 | 30.0 | nan | 30 | 48 | 79 | 15 | 33 | 90.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
75 | 1458 | 46 | 0.00053011 | TGCGCTGTGAGGCGGAGCTGGGGGAAATTGCAGTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-1*00(358.3) | nan | TRAJ24*00(258) | TRAC*00(64.9) | 261 | 272 | 295 | 0 | 11 | 55.0 | nan | 31 | 52 | 83 | 15 | 36 | 105.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
76 | 1448 | 43 | 0.000526474 | TGTGCTGCCATGGATAGCAACTATCAGTTAATCTGG | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-2*00(343.2) | nan | TRAJ33*00(297.7) | TRAC*00(64.8) | 261 | 268 | 295 | 0 | 7 | 35.0 | nan | 17 | 46 | 77 | 7 | 36 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
77 | 1446 | 7 | 0.000525746 | TGCCTCGTGGGTGGAGGATTCATGGAATATGGAAACAAACTGGTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(270.4) | nan | TRAJ47*00(296.6) | TRAC*00(64.5) | 261 | 274 | 297 | 0 | 13 | 65.0 | nan | 14 | 46 | 77 | 17 | 48 | DC18 | 143.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
78 | 1436 | 10 | 0.000522111 | TGCATCGTCAGACGACTGTCTAACGACTACAAGCTCAGCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-1*00(278) | nan | TRAJ20*00(273.7) | TRAC*00(65.1) | 261 | 273 | 297 | 0 | 12 | 60.0 | nan | 22 | 46 | 77 | 18 | 42 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
79 | 1435 | 18 | 0.000521747 | TGTGCTGTGAGTTCCTCAGGGACACAGGCTTTCAGAAACTTGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-4*00(263.6) | nan | TRAJ8*00(285.4) | TRAC*00(64.5) | 270 | 282 | 304 | 0 | 12 | 60.0 | nan | 17 | 49 | 80 | 15 | 47 | ST20GSA22G | 132.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
80 | 1434 | 19 | 0.000521383 | TGTGCCTTTTCAGATGGCCAGAAGCTGCTCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV24*00(258.9) | nan | TRAJ16*00(287.7) | TRAC*00(95.4) | 267 | 276 | 297 | 0 | 9 | 45.0 | nan | 22 | 49 | 80 | 6 | 33 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
81 | 1382 | 15 | 0.000502477 | TGTGCTGTGAGAGACATGGGGCCAACGGGCAGGAGAGCACTTACTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV3*00(297.2) | nan | TRAJ5*00(273.7) | TRAC*00(64.8) | 270 | 288 | 306 | 0 | 18 | 90.0 | nan | 25 | 49 | 80 | 24 | 48 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
82 | 1361 | 15 | 0.000494842 | TGTGTTGTGACTCAGGGAACTGGAGGCTTCAAAACTATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-2*00(251.6) | nan | TRAJ9*00(272.9) | TRAC*00(69) | 270 | 280 | 304 | 0 | 10 | 50.0 | nan | 26 | 50 | 81 | 18 | 42 | 120.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
83 | 1340 | 101 | 0.000487206 | TGTGCAGAGGGTTTTCTAACGACTACAAGCTCAGCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV5*00(277.3) | nan | TRAJ20*00(279.2) | TRAC*00(64.8) | 264 | 276 | 297 | 0 | 12 | SA273G | 46.0 | nan | 21 | 46 | 77 | 13 | 38 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
84 | 1329 | 4 | 0.000483207 | TGCCTAATAAAGGAAATGGGAGCCCAGAAGCTGGTATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV40*00(242.2) | nan | TRAJ54*00(262.7) | TRAC*00(65) | 249 | 254 | 283 | 0 | 5 | ST251C | 11.0 | nan | 27 | 49 | 80 | 17 | 39 | 110.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
85 | 1316 | 37 | 0.00047848 | TGTGCTGTGGGTGCTCATGGAGGAAGCCAAGGAAATCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV8-3*00(273.4) | nan | TRAJ42*00(298.7) | TRAC*00(64.8) | 270 | 287 | 304 | 0 | 17 | SG284T | 71.0 | nan | 26 | 55 | 86 | 16 | 45 | 145.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
86 | 1304 | 41 | 0.000474117 | TGTGCTGTGATGGATAGCAACTATCAGTTAATCTGG | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-2*00(355.9) | nan | TRAJ33*00(287.3) | TRAC*00(64.6) | 261 | 271 | 295 | 0 | 10 | 50.0 | nan | 19 | 46 | 77 | 9 | 36 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
87 | 1300 | 7 | 0.000472663 | TGCCTCGTGGTCCCCTCCGAGCGGAGAAAGCTTATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(259.5) | nan | TRAJ23*00(226.4) | TRAC*00(78.3) | 261 | 271 | 297 | 0 | 10 | 50.0 | nan | 34 | 52 | 83 | 22 | 39 | DG38 | 73.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
88 | 1294 | 37 | 0.000470481 | TGCGCTGTGAGTCGGGATAGCAGCTATAAATTGATCTTC | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV1-1*00(356.3) | nan | TRAJ12*00(278) | TRAC*00(64.6) | 261 | 272 | 295 | 0 | 11 | 55.0 | nan | 24 | 49 | 80 | 14 | 39 | 125.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
89 | 1285 | 54 | 0.000467209 | TGTGCTGTGGTGGCTCTAGCAACACAGGCAAACTAATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV2*00(272.8) | nan | TRAJ37*00(279.2) | TRAC*00(65) | 249 | 265 | 283 | 0 | 16 | SA259TSA262C | 52.0 | nan | 20 | 51 | 82 | 10 | 41 | SG27A | 141.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan |
90 | 1284 | 19 | 0.000466845 | TGTGTGGTGAACCCCACCGACAAGCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-1*00(303.8) | nan | TRAJ34*00(248.4) | TRAC*00(65.5) | 261 | 273 | 294 | 0 | 12 | 60.0 | nan | 28 | 47 | 78 | 14 | 33 | 95.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
91 | 1276 | 4 | 0.000463937 | TGCCTCGTGGTCAAGGGGGGAAGCCAAGGAAATCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV4*00(258.3) | nan | TRAJ42*00(274.2) | TRAC*00(71.5) | 261 | 271 | 297 | 0 | 10 | 50.0 | nan | 28 | 55 | 86 | 15 | 42 | SA30G | 121.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
92 | 1268 | 9 | 0.000461028 | TGTGCAATGAGGCCGAATTCAGGATACAGCACCCTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV14DV4*00(256.1) | nan | TRAJ11*00(292.8) | TRAC*00(65.1) | 273 | 284 | 310 | 0 | 11 | 55.0 | nan | 21 | 49 | 80 | 14 | 42 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
93 | 1260 | 10 | 0.000458119 | TGTGCAACGTTCGAATTCAGGATACAGCACCCTCACCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV14DV4*00(239.4) | nan | TRAJ11*00(292.8) | TRAC*00(63.8) | 273 | 280 | 310 | 0 | 7 | 35.0 | nan | 21 | 49 | 80 | 12 | 40 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
94 | 1250 | 19 | 0.000454483 | TGTGCCGTGAAAGGGGGGAACAACAGACTCGCTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV12-2*00(294.9) | nan | TRAJ7*00(258.7) | TRAC*00(64.5) | 264 | 275 | 297 | 0 | 11 | 55.0 | nan | 27 | 48 | 79 | 15 | 36 | 105.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
95 | 1233 | 4 | 0.000448303 | TGTCTTCTGGGAGATTTCACTGGAGGTAGCAACTATAAACTGACATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV40*00(304) | nan | TRAJ53*00(293.4) | TRAC*00(64.7) | 249 | 267 | 283 | 0 | 18 | SC264T | 76.0 | nan | 27 | 55 | 86 | 20 | 48 | 140.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | |
96 | 1230 | 7 | 0.000447212 | TGCATCGTCAGAGCACGGAATTATGGAGGAAGCCAAGGAAATCTCATCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-1*00(289) | nan | TRAJ42*00(322.6) | TRAC*00(64.8) | 261 | 274 | 297 | 0 | 13 | 65.0 | nan | 21 | 55 | 86 | 17 | 51 | 170.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
97 | 1199 | 7 | 0.000435941 | TGCATCGTCAGAGTGGGGGACATGCGCTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV26-1*00(293.1) | nan | TRAJ43*00(213.9) | TRAC*00(107.1) | 261 | 275 | 297 | 0 | 14 | 70.0 | nan | 31 | 43 | 74 | 18 | 30 | 60.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
98 | 1182 | 17 | 0.00042976 | TGTGCAGGAGTAGGGGGAGGTAGCAACTATAAACTGACATTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV27*00(266.9) | nan | TRAJ53*00(288.4) | TRAC*00(64.8) | 264 | 274 | 294 | 0 | 10 | 50.0 | nan | 28 | 55 | 86 | 15 | 42 | 135.0 | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | nan | ||
99 | 1168 | 12 | 0.000424669 | TGTGCAATGAGCTGGTTTTCTGGTGGCTACAATAAGCTGATTTTT | NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN | TRAV14DV4*00(257.5) | nan | TRAJ4*00(307.7) | TRAC*00(65.1) | 273 | 284 | 310 | 0 | 11 | 55.0 | nan | 21 | 52 | 83 | 14 | 45 | 155.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 milab-human-rna-tcr-umi-multiplex \
{} \
{=s:R1:R2:=} \
{=s:.*/:results/:;s:_R.*::=}"
Under the hood pipeline
Under the hood milab-human-rna-tcr-umi-multiplex
preset 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=false \
--tag-pattern "^N{0:2}tggtatcaacgcagagt(UMI:N{14})N{21}(R1:*) \ ^N{20}(R2:*)" \
--report results/Multi_TRA_FS115_1.report.txt \
--json-report results/Multi_TRA_FS115_1.report.json \
raw/Multi_TRA_FS115_1_R1.fastq.gz \
raw/Multi_TRA_FS115_1_R2.fastq.gz \
results/Multi_TRA_FS115_1.vdjca
Option --report
is specified here explicitly.
--species hsa
- determines the organism species (hsa for Homo Sapiens).
-p generic-amplicon
- a generic preset of MiXCR parameters for amplicon data.
-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 primer sequences were trimmed with
--tagPattern
. -OjParameters.parameters.floatingRightBound=false
- Results in a global alignment algorithm for J gene right bound due to the absence of primer sequences.
-OcParameters.parameters.floatingRightBound=false
- Results in a local alignment algorithm for C gene right bound.We use it because primer sequences were trimmed with
--tagPattern
. --tagPattern "^N{0:2}tggtatcaacgcagagt(UMI:N{14})N{21}(R1:*) \ ^N{20}(R2:*)"
- tag pattern specifies the location of the UMI barcode according to the library structure and trims primers sequences.
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/Multi_TRA_FS115_1.report.txt \
--json-report results/Multi_TRA_FS115_1.report.json \
results/Multi_TRA_FS115_1.vdjca \
results/Multi_TRA_FS115_1.refined.vdjca
assemble
Assembles 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/Multi_TRA_FS115_1.report.txt \
--json-report results/Multi_TRA_FS115_1.report.json \
results/Multi_TRA_FS115_1.refined.vdjca \
results/Multi_TRA_FS115_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/Multi_TRA_FS115_1.clns \
results/Multi_TRA_FS115_1.clones.TRA.txt
-с TRA
- 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
From this plot we can tell that all samples have high alignment rate, more than 90% of reads have been successfully aligned to the reference sequences and CDR3 has been established.
Now we can check chain distribution plot:
mixcr exportQc chainUsage results/*.clns figs/chainUsage.pdf
Here we can see that no cross-contamination has occurred. TCR alpha samples consist only of TRA clones, and TCR betta only of TRB.
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/Multi_TRA_FS115_1.clns \
results/Multi_TRA_FS115_1.report.txt
mixcr exportReports \
--json \
results/Multi_TRA_FS115_1.clns \
results/Multi_TRA_FS115_1.report.json
Show report file
============== Align Report ==============
============== Align Report ==============
Total sequencing reads: 3567679
Successfully aligned reads: 3448515 (96.66%)
Chimeras: 4 (0%)
Paired-end alignment conflicts eliminated: 80984 (2.27%)
Alignment failed, no hits (not TCR/IG?): 2855 (0.08%)
Alignment failed because of absence of V hits: 24752 (0.69%)
Alignment failed because of absence of J hits: 62293 (1.75%)
No target with both V and J alignments: 11677 (0.33%)
Alignment failed because of low total score: 157 (0%)
Absent barcode: 17430 (0.49%)
Overlapped: 3272356 (91.72%)
Overlapped and aligned: 3218859 (90.22%)
Alignment-aided overlaps: 3992 (0.12%)
Overlapped and not aligned: 53497 (1.5%)
No CDR3 parts alignments, percent of successfully aligned: 1211 (0.04%)
Partial aligned reads, percent of successfully aligned: 34409 (1%)
V gene chimeras: 57 (0%)
J gene chimeras: 2004 (0.06%)
TRA chains: 3448272 (99.99%)
TRA non-functional: 355157 (10.3%)
TRB chains: 111 (0%)
TRB non-functional: 20 (18.02%)
TRD chains: 128 (0%)
TRD non-functional: 31 (24.22%)
Realigned with forced non-floating bound: 563770 (15.8%)
Realigned with forced non-floating right bound in left read: 43576 (1.22%)
Realigned with forced non-floating left bound in right read: 43576 (1.22%)
Tag parsing report:
Total reads: 3567679
Matched reads: 3550249 (99.51%)
Projection +R1 +R2: 3550249 (99.51%)
For variant 0:
For projection [1, 2]:
UMI:Left position:
17: + 982994 (27.69%) = 982994 (27.69%)
18: + 1232085 (34.7%) = 2215079 (62.39%)
19: + 1335170 (37.61%) = 3550249 (100%)
UMI:Right position:
31: + 982994 (27.69%) = 982994 (27.69%)
32: + 1232085 (34.7%) = 2215079 (62.39%)
33: + 1335170 (37.61%) = 3550249 (100%)
R1:Left position:
52: + 982994 (27.69%) = 982994 (27.69%)
53: + 1232085 (34.7%) = 2215079 (62.39%)
54: + 1335170 (37.61%) = 3550249 (100%)
R1:Right position:
52~146: + 255675 (7.2%) = 255675 (7.2%)
147: + 659338 (18.57%) = 915013 (25.77%)
148: + 1760294 (49.58%) = 2675307 (75.36%)
149: + 874942 (24.64%) = 3550249 (100%)
R2:Left position: 20
Variants: 0
Cost:
0: + 3433842 (96.72%) = 3433842 (96.72%)
5~10: + 116407 (3.28%) = 3550249 (100%)
UMI length: 14
R1 length:
0~93: + 179299 (5.05%) = 179299 (5.05%)
94: + 1523494 (42.91%) = 1702793 (47.96%)
95: + 951209 (26.79%) = 2654002 (74.76%)
96~97: + 896247 (25.24%) = 3550249 (100%)
R2 length:
15~127: + 291090 (8.2%) = 291090 (8.2%)
128: + 1098902 (30.95%) = 1389992 (39.15%)
129: + 2160257 (60.85%) = 3550249 (100%)
============== RefineTagsAndSort Report ==============
============== RefineTagsAndSort Report ==============
Time spent in correction: 0ns
UMI input diversity: 266703
UMI output diversity: 137060 (51.39%)
UMI input reads: 3448498
UMI output count: 3448497 (100%)
UMI mean reads per tag: 12.93
UMI input core diversity: 266674 (99.99%)
UMI input core reads: 3448469 (100%)
UMI directly corrected diversity: 129642 (48.61%)
UMI directly corrected reads: 144398 (4.19%)
UMI diversity filtered by tag quality: 1 (0%)
UMI reads filtered by tag quality: 1 (0%)
UMI diversity filtered by whitelist: 0 (0%)
UMI recursively corrected: 6799
Number of output records: 3302573 (95.77%)
Filter report:
Number of groups: 137060
Number of groups accepted: 53937 (39.35%)
Total records weight: 3448497.0
Records weight accepted: 3302573.0 (95.77%)
Operator #0:
Effective threshold: 8.0
============== Assemble Report ==============
============== Assemble Report ==============
Number of input groups: 53937
Number of input alignments: 3302573
Number of output pre-clonotypes: 50730
Number of clonotypes per group:
0: + 3239 (6.01%) = 3239 (6.01%)
1: + 50558 (93.83%) = 53797 (99.84%)
2: + 86 (0.16%) = 53883 (100%)
Number of core alignments: 2748073 (83.21%)
Discarded core alignments: 504431 (18.36%)
Empirically assigned alignments: 8264 (0.25%)
Empirical assignment conflicts: 0 (0%)
UMI+VJ-gene empirically assigned alignments: 8264 (0.25%)
VJ-gene empirically assigned alignments: 0 (0%)
UMI empirically assigned alignments: 0 (0%)
Number of ambiguous UMIs: 86
Number of ambiguous V-genes: 3
Number of ambiguous J-genes: 1
Number of ambiguous UMI+V/J-gene combinations: 4
Unassigned alignments: 545221 (16.51%)
Final clonotype count: 27957
Average number of reads per clonotype: 98.38
Reads used in clonotypes, percent of total: 2750375 (77.09%)
Reads used in clonotypes before clustering, percent of total: 2755131 (77.22%)
Number of reads used as a core, percent of used: 2755073 (100%)
Mapped low quality reads, percent of used: 58 (0%)
Reads clustered in PCR error correction, percent of used: 4756 (0.17%)
Reads pre-clustered due to the similar VJC-lists, percent of used: 30 (0%)
Reads dropped due to the lack of a clone sequence, percent of total: 50069 (1.4%)
Reads dropped due to a too short clonal sequence, percent of total: 52 (0%)
Reads dropped due to low quality, percent of total: 0 (0%)
Reads dropped due to failed mapping, percent of total: 1154 (0.03%)
Reads dropped with low quality clones, percent of total: 0 (0%)
Clonotypes eliminated by PCR error correction: 88
Clonotypes dropped as low quality: 0
Clonotypes pre-clustered due to the similar VJC-lists: 1
TRA chains: 27956 (100%)
TRA non-functional: 2833 (10.13%)
TRD chains: 1 (0%)
TRD non-functional: 1 (100%)
{
"align": {
"type": "alignerReport",
"commandLine": "--report results/Multi_TRA_FS115_1.align.report.txt --json-report results/Multi_TRA_FS115_1.align.report.json --preset milab-human-tcr-rna-multiplex-cdr3 /raw/2s3/milab-tcr-mulriplex/Multi_TRA_FS115_1_S149_R1_001.fastq.gz /raw/2s3/milab-tcr-mulriplex/Multi_TRA_FS115_1_S149_R2_001.fastq.gz results/Multi_TRA_FS115_1.vdjca",
"inputFiles": [
"/raw/2s3/milab-tcr-mulriplex/Multi_TRA_FS115_1_S149_R1_001.fastq.gz",
"/raw/2s3/milab-tcr-mulriplex/Multi_TRA_FS115_1_S149_R2_001.fastq.gz"
],
"outputFiles": [
"results/Multi_TRA_FS115_1.vdjca"
],
"version": "4.0.0-432-protocols; built=Tue Oct 11 09:48:35 CEST 2022; rev=2b3b93fc41; lib=repseqio.v2.0",
"trimmingReport": null,
"totalReadsProcessed": 3567679,
"aligned": 3448515,
"notAligned": 119164,
"notAlignedReasons": {
"LowTotalScore": 157,
"NoVHits": 24752,
"NoJHits": 62293,
"NoBarcode": 17430,
"NoCDR3Parts": 0,
"VAndJOnDifferentTargets": 11677,
"NoHits": 2855
},
"chimeras": 4,
"overlapped": 3272356,
"alignmentAidedOverlaps": 3992,
"overlappedAligned": 3218859,
"overlappedNotAligned": 53497,
"pairedEndAlignmentConflicts": 80984,
"vChimeras": 57,
"jChimeras": 2004,
"chainUsage": {
"type": "chainUsage",
"chimeras": 4,
"total": 3448515,
"chains": {
"TRA": {
"total": 3448272,
"nonFunctional": 355157,
"isOOF": 319826,
"hasStops": 35331
},
"TRB": {
"total": 111,
"nonFunctional": 20,
"isOOF": 6,
"hasStops": 14
},
"TRD": {
"total": 128,
"nonFunctional": 31,
"isOOF": 30,
"hasStops": 1
}
}
},
"realignedWithForcedNonFloatingBound": 563770,
"realignedWithForcedNonFloatingRightBoundInLeftRead": 43576,
"realignedWithForcedNonFloatingLeftBoundInRightRead": 43576,
"noCDR3PartsAlignments": 1211,
"partialAlignments": 34409,
"tagParsingReport": {
"total": 3567679,
"matched": 3550249,
"totalBitCost": 595287.5547334387,
"projections": {
"1,2": 3550249
},
"detailedReport": [
{
"nested": {
"1,2": {
"variantId": 0,
"nested": [
[
{
"positionDistributions": {
"UMI:Left": {
"17": 982994,
"18": 1232085,
"19": 1335170
},
"UMI:Right": {
"32": 1232085,
"33": 1335170,
"31": 982994
},
"R1:Left": {
"52": 982994,
"53": 1232085,
"54": 1335170
},
"R1:Right": {
"128": 199,
"129": 217,
"130": 56,
"131": 45,
"132": 304,
"133": 556,
"134": 1335,
"135": 1162,
"136": 1327,
"137": 493,
"138": 640,
"139": 668,
"140": 639,
"141": 1721,
"142": 2101,
"143": 2310,
"144": 1097,
"145": 17799,
"146": 200055,
"147": 659338,
"148": 1760294,
"149": 874942,
"52": 1,
"53": 1,
"54": 4,
"55": 13,
"56": 18,
"58": 13,
"59": 12,
"60": 12,
"66": 5,
"67": 11,
"68": 31,
"69": 31,
"70": 1,
"71": 5,
"72": 6,
"73": 39,
"74": 451,
"75": 685,
"76": 811,
"77": 200,
"78": 11,
"79": 37,
"80": 64,
"81": 45,
"82": 8,
"83": 13,
"84": 27,
"85": 26,
"86": 142,
"87": 185,
"88": 177,
"89": 26,
"90": 26,
"91": 50,
"92": 152,
"93": 222,
"94": 248,
"95": 26,
"96": 4,
"97": 28,
"98": 239,
"99": 1421,
"100": 2175,
"101": 5016,
"102": 4388,
"103": 4156,
"104": 127,
"105": 22,
"106": 37,
"107": 38,
"108": 36,
"109": 20,
"110": 1,
"111": 5,
"113": 1,
"114": 32,
"115": 27,
"116": 35,
"117": 8,
"118": 15,
"119": 224,
"120": 266,
"121": 266,
"122": 35,
"123": 29,
"124": 114,
"125": 104,
"126": 121,
"127": 126
}
},
"lengthDistributions": {
"UMI": {
"14": 3550249
},
"R1": {
"0": 1,
"1": 1,
"2": 32,
"3": 1,
"4": 4,
"5": 4,
"6": 31,
"14": 11,
"15": 66,
"16": 1,
"17": 1,
"18": 2,
"19": 10,
"20": 8,
"21": 119,
"22": 1522,
"23": 528,
"24": 3,
"25": 17,
"26": 34,
"27": 110,
"29": 1,
"30": 44,
"31": 1,
"32": 8,
"33": 65,
"34": 453,
"35": 7,
"36": 9,
"37": 78,
"38": 8,
"39": 36,
"40": 558,
"41": 59,
"42": 2,
"43": 3,
"44": 17,
"45": 87,
"46": 658,
"47": 4287,
"48": 1740,
"49": 10472,
"50": 265,
"51": 28,
"52": 29,
"53": 17,
"54": 58,
"55": 46,
"57": 5,
"58": 1,
"59": 1,
"61": 1,
"62": 86,
"63": 12,
"64": 3,
"65": 6,
"66": 47,
"67": 714,
"68": 10,
"69": 5,
"70": 81,
"71": 19,
"72": 273,
"73": 11,
"74": 60,
"75": 373,
"76": 164,
"77": 14,
"78": 21,
"79": 91,
"80": 1196,
"81": 187,
"82": 2947,
"83": 194,
"84": 543,
"85": 939,
"86": 292,
"87": 719,
"88": 984,
"89": 4240,
"90": 1525,
"91": 14789,
"92": 33248,
"93": 93956,
"94": 1523494,
"95": 951209,
"96": 509864,
"97": 386383
}
},
"costDistributions": {
"0": 3433842,
"5": 109648,
"10": 6759
},
"variantDistribution": {
"0": 3550249
}
}
],
[
{
"positionDistributions": {
"R2:Left": {
"20": 3550249
},
"R2:Right": {
"144": 422,
"145": 3921,
"146": 137570,
"147": 149168,
"35": 1,
"148": 1098902,
"149": 2160257,
"133": 1,
"138": 2,
"139": 2,
"142": 2,
"143": 1
}
},
"lengthDistributions": {
"R2": {
"128": 1098902,
"129": 2160257,
"15": 1,
"113": 1,
"118": 2,
"119": 2,
"122": 2,
"123": 1,
"124": 422,
"125": 3921,
"126": 137570,
"127": 149168
}
},
"costDistributions": {
"0": 3550249
},
"variantDistribution": {
"0": 3550249
}
}
]
],
"costDistributions": {
"0": 3433842,
"5": 109648,
"10": 6759
}
}
}
}
]
}
},
"refineTagsAndSort": {
"type": "refineTagsAndSort",
"commandLine": "--report results/Multi_TRA_FS115_1.refine.report.txt --json-report results/Multi_TRA_FS115_1.refine.report.json results/Multi_TRA_FS115_1.vdjca results/Multi_TRA_FS115_1.refined.vdjca",
"inputFiles": [
"results/Multi_TRA_FS115_1.vdjca"
],
"outputFiles": [
"results/Multi_TRA_FS115_1.refined.vdjca"
],
"version": "4.0.0-432-protocols; built=Tue Oct 11 09:48:35 CEST 2022; rev=2b3b93fc41; lib=repseqio.v2.0",
"correctionReport": {
"inputRecords": 3448498,
"outputRecords": 3302573,
"steps": [
{
"tagName": "UMI",
"inputGroups": 1,
"inputDiversity": 266703,
"inputCount": 3448498,
"coreDiversity": 266674,
"coreCount": 3448469,
"directlyCorrectedDiversity": 129642,
"directlyCorrectedCount": 144398,
"filteredDiversity": 1,
"filteredCount": 1,
"recursivelyCorrected": 6799,
"diversityFilteredByWhitelist": 0,
"outputDiversity": 137060,
"outputCount": 3448497
}
],
"filterReport": {
"type": "filter_groups_report",
"groupingKeys": [
"UMI"
],
"numberOfGroups": 137060,
"numberOfGroupsAccepted": 53937,
"totalWeight": 3448497,
"totalWeightAccepted": 3302573,
"operatorReports": [
{
"operatorReport": {
"type": "generic_hist_report",
"threshold": 8
},
"metricHists": [
{
"bins": [
{
"from": 1,
"to": 1.5848931924611136,
"weight": 58017
},
{
"from": 1.5848931924611136,
"to": 2.51188643150958,
"weight": 9696
},
{
"from": 2.51188643150958,
"to": 3.981071705534973,
"weight": 5234
},
{
"from": 3.981071705534973,
"to": 6.309573444801933,
"weight": 8426
},
{
"from": 6.309573444801933,
"to": 10,
"weight": 4559
},
{
"from": 10,
"to": 15.848931924611142,
"weight": 6605
},
{
"from": 15.848931924611142,
"to": 25.11886431509581,
"weight": 8135
},
{
"from": 25.11886431509581,
"to": 39.810717055349734,
"weight": 8412
},
{
"from": 39.810717055349734,
"to": 63.09573444801933,
"weight": 10494
},
{
"from": 63.09573444801933,
"to": 100,
"weight": 7942
},
{
"from": 100,
"to": 158.48931924611142,
"weight": 6323
},
{
"from": 158.48931924611142,
"to": 251.18864315095823,
"weight": 2185
},
{
"from": 251.18864315095823,
"to": 398.1071705534973,
"weight": 558
},
{
"from": 398.1071705534973,
"to": 630.9573444801937,
"weight": 472
},
{
"from": 630.9573444801937,
"to": 1000,
"weight": 2
}
],
"collectionSpec": {
"log": true,
"binNumber": 0,
"minBinWidth": 0.2
}
}
]
}
]
}
}
},
"assemble": {
"type": "assemblerReport",
"commandLine": "--report results/Multi_TRA_FS115_1.assemble.report.txt --json-report results/Multi_TRA_FS115_1.assemble.report.json results/Multi_TRA_FS115_1.refined.vdjca results/Multi_TRA_FS115_1.clns",
"inputFiles": [
"results/Multi_TRA_FS115_1.refined.vdjca"
],
"outputFiles": [
"results/Multi_TRA_FS115_1.clns"
],
"version": "4.0.0-432-protocols; built=Tue Oct 11 09:48:35 CEST 2022; rev=2b3b93fc41; lib=repseqio.v2.0",
"preCloneAssemblerReport": {
"type": "preCloneAssemblerReport",
"inputGroups": 53937,
"inputAlignments": 3302573,
"clonotypes": 50730,
"clonotypesPerGroup": {
"0": 3239,
"1": 50558,
"2": 86
},
"coreAlignments": 2748073,
"discardedCoreAlignments": 504431,
"empiricallyAssignedAlignments": 8264,
"vjEmpiricallyAssignedAlignments": 0,
"umiEmpiricallyAssignedAlignments": 0,
"gatEmpiricallyAssignedAlignments": 8264,
"empiricalAssignmentConflicts": 0,
"unassignedAlignments": 545221,
"umiConflicts": 86,
"gatConflicts": 4,
"geneConflicts": {
"Variable": 3,
"Joining": 1
},
"coreClonotypesDroppedByTagSuffix": 0,
"coreAlignmentsDroppedByTagSuffix": 0
},
"totalReadsProcessed": 3567679,
"initialClonesCreated": 28046,
"readsDroppedNoTargetSequence": 50069,
"readsDroppedTooShortClonalSequence": 52,
"readsDroppedLowQuality": 0,
"coreReads": 2755073,
"readsDroppedFailedMapping": 1154,
"lowQualityRescued": 58,
"clonesClustered": 88,
"readsClustered": 4756,
"clones": 27957,
"clonesDroppedAsLowQuality": 0,
"clonesPreClustered": 1,
"readsPreClustered": 30,
"readsInClones": 2750375,
"readsInClonesBeforeClustering": 2755131,
"readsDroppedWithLowQualityClones": 0,
"clonalChainUsage": {
"type": "chainUsage",
"chimeras": 0,
"total": 27957,
"chains": {
"TRA": {
"total": 27956,
"nonFunctional": 2833,
"isOOF": 2717,
"hasStops": 116
},
"TRD": {
"total": 1,
"nonFunctional": 1,
"isOOF": 1,
"hasStops": 0
}
}
}
}
}