Giter VIP home page Giter VIP logo

Comments (10)

Januaryyiyue avatar Januaryyiyue commented on July 22, 2024

To clarify, alignment worked but the LOG file for the assembly did not work

from mixcr.

mizraelson avatar mizraelson commented on July 22, 2024

Hi, can you share the commands you used and the environment setting (server, SLURM, etc)?

from mixcr.

Januaryyiyue avatar Januaryyiyue commented on July 22, 2024

Sure! I used sbatch on a .sh script with -t 5-00:00:00 --mem=89G. The script is pasted below.
I ran 26 samples; only 2 of them had this problem. The log align and .clns files were successfully generated, but the log assemble and clones files were absent. I am using capTCR-seq data.

#!/bin/bash
#
#$ -cwd
module load java/8
module load mixcr/4.1.2
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar align -p rnaseq-cdr3 --verbose -s hsa -OallowPartialAlignments=true -OvParameters.geneFeatureToAlign=VGeneWithP -r /cluster/path/to/mixcr4/log/LOG_ALIGN_file.txt /cluster/path/to/data/file_R1.fastq /cluster/path/to/data/file_R2.fastq /cluster/path/to/mixcr4/alignments/alignments_file.vdjca
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar assemblePartial /cluster/path/to/mixcr4/alignments/alignments_file.vdjca /cluster/path/to/mixcr4/assemble/rescue1_file.vdjca
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar assemblePartial /cluster/path/to/mixcr4/assemble/rescue1_file.vdjca /cluster/path/to/mixcr4/assemble/rescue2_file.vdjca
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar extend /cluster/path/to/mixcr4/assemble/rescue2_file.vdjca /cluster/path/to/mixcr4/assemble/extended_file.vdjca
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar assemble -r /cluster/path/to/mixcr4/log/LOG_ASSEMBLE_file.txt /cluster/path/to/mixcr4/assemble/extended_file.vdjca /cluster/path/to/mixcr4/assemble/file.clns
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar exportClones --chains TRA /cluster/path/to/mixcr4/assemble/file.clns /cluster/path/to/mixcr4/clones/CLONES_TRAfile.tsv
java -jar /cluster/path/to/mixcr/4.1.2/mixcr.jar exportClones --chains TRB /cluster/path/to/mixcr4/assemble/file.clns /cluster/path/to/mixcr4/clones/CLONES_TRBfile.tsv

Here's the result for the last few lines of the log align file:

Analysis time: 1128.7m
Total sequencing reads: 43390557
Successfully aligned reads: 38195684 (88.03%)
Chimeras: 5 (0%)
Paired-end alignment conflicts eliminated: 93125 (0.21%)
Alignment failed, no hits (not TCR/IG?): 1233379 (2.84%)
Alignment failed because of absence of CDR3 parts: 2129356 (4.91%)
Alignment failed because of low total score: 1832138 (4.22%)
Overlapped: 28612588 (65.94%)
Overlapped and aligned: 26198878 (60.38%)
Alignment-aided overlaps: 1127455 (4.3%)
Overlapped and not aligned: 2413710 (5.56%)
No CDR3 parts alignments, percent of successfully aligned: 44345 (0.12%)
Partial aligned reads, percent of successfully aligned: 21613444 (56.59%)
V gene chimeras: 12995 (0.03%)
J gene chimeras: 429656 (0.99%)
TRA chains: 14828702 (38.82%)
TRA non-functional: 5688903 (38.36%)
TRB chains: 4440587 (11.63%)
TRB non-functional: 570282 (12.84%)
TRD chains: 333584 (0.87%)
TRD non-functional: 59217 (17.75%)
TRG chains: 4814406 (12.6%)
TRG non-functional: 2683093 (55.73%)
TRAD chains: 13778400 (36.07%)
TRAD non-functional: 0 (0%)
Realigned with forced non-floating bound: 31810848 (73.31%)
Realigned with forced non-floating right bound in left read: 167045 (0.38%)
Realigned with forced non-floating left bound in right read: 167045 (0.38%)

from mixcr.

mizraelson avatar mizraelson commented on July 22, 2024

Can you please run it with the latest MiXCR v4.6 and confirm the issue persists?

from mixcr.

Januaryyiyue avatar Januaryyiyue commented on July 22, 2024

Hi, we don't have MiXCR v4.6 on our cluster, and I don't have control over that, sorry about that..

from mixcr.

Januaryyiyue avatar Januaryyiyue commented on July 22, 2024

Could it be related to the input file size? For the ones that had this error the input file_R1.fastq.gz and file_R2.fastq.gz are both >3G each

from mixcr.

mizraelson avatar mizraelson commented on July 22, 2024

It looks like an issue that has been fixed in a later version. Shouldn't have anything to do with the size.

from mixcr.

mizraelson avatar mizraelson commented on July 22, 2024

I wanted to follow up and ask if you could share the problematic sample. That way, I can test whether the latest MiXCR version throws the error or not.

from mixcr.

Januaryyiyue avatar Januaryyiyue commented on July 22, 2024

Hi! Thank you for the follow up. I ran the analysis again with the exact same command, also with v4.1.2, and the error disappeared. This happens sometimes with the licensing error issue as well. Everything is good now, thank you!

from mixcr.

mizraelson avatar mizraelson commented on July 22, 2024

Great!

from mixcr.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.