-------------------------------------------------------------------------------- [log] Date: 2020-08-01 [log] Sample: sample_1_fished [log] Input file(s): sample_1_fished.fastq sample_2_fished.fastq [log] Reference: 7ca733b18edab1fc837a34e3764b7d55e8dc93b0 -------------------------------------------------------------------------------- [alignment] Analyzing read length [alignment] Pseudoaligning with Kallisto: kallisto pseudo -i /projects/rwarren_prj2/bin/arcasHLA/scripts/../dat/ref/hla.idx -t 8 -o /tmp/arcas_d6b8f8dd-df3b-4ffc-87dc-4c802b929d2b/ sample_1_fished.fastq sample_2_fished.fastq [quant] fragment length distribution will be estimated from the data [index] k-mer length: 31 [index] number of targets: 11,226 [index] number of k-mers: 264,111 [index] number of equivalence classes: 25,101 [quant] running in paired-end mode [quant] will process pair 1: sample_1_fished.fastq sample_2_fished.fastq [quant] finding pseudoalignments for the reads ... done [quant] processed 273 reads, 23 reads pseudoaligned [alignment] Processing pseudoalignment [alignment] Processed 573 reads, 23 pseudoaligned to HLA reference [alignment] 23 reads mapped to a single HLA gene [alignment] Observed HLA genes: gene abundance read count classes HLA-A 11.70% 3 3 HLA-B 15.73% 4 4 HLA-E 3.95% 1 1 HLA-G 16.85% 4 4 HLA-J 23.47% 6 3 HLA-L 11.72% 3 2 HLA-V 16.57% 2 1 -------------------------------------------------------------------------------- [genotype] Genotyping parameters: population: prior minimum count: 75 max iterations: 1000 tolerance: 1e-06 drop iterations: 20 drop threshold: 0.1 zygosity threshold: 0.15 -------------------------------------------------------------------------------- [genotype] Genotyping HLA-A [genotype] Not enough reads aligned to HLA-A to genotype. -------------------------------------------------------------------------------- [genotype] Genotyping HLA-B [genotype] Not enough reads aligned to HLA-B to genotype. -------------------------------------------------------------------------------- [genotype] Genotyping HLA-C [genotype] Not enough reads aligned to HLA-C to genotype. -------------------------------------------------------------------------------- [genotype] Genotyping HLA-DPA1 [genotype] Not enough reads aligned to HLA-DPA1 to genotype. -------------------------------------------------------------------------------- [genotype] Genotyping HLA-DPB1 [genotype] Not enough reads aligned to HLA-DPB1 to genotype. -------------------------------------------------------------------------------- [genotype] Genotyping HLA-DQA1 [genotype] Not enough reads aligned to HLA-DQA1 to genotype. -------------------------------------------------------------------------------- [genotype] Genotyping HLA-DQB1 [genotype] Not enough reads aligned to HLA-DQB1 to genotype. --------------------------------------------------------------------------------