-
Notifications
You must be signed in to change notification settings - Fork 35
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
no output #484
Comments
Could you please attach the log file the program generated when it ran? |
[akulan@cn2877 V1_outs_090624]$ /home/akulan/Drop-seq_tools-2.5.1/AssignCellsToSamples -m 24g -INPUT_BAM V1_possorted_genome_bam_tagged.bam -CELL_BC_FILE /data/GBMAD/McMahon-Claudia_SC_083024/V1_outs_090624/V1_filtered_feature_bc_matrix/barcodes.tsv.gz -VCF /data/akulan/gsa_may2024/PLINK_130524_0118/GSA_051324_scClaudiaV1labelled_chrAdded_snpsOnly.vcf.gz -OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.txt -VCF_OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.vcf.gz -CELL_BARCODE_TAG CB -LOCUS_FUNCTION_LIST INTRONIC -IGNORED_CHROMOSOMES null -IGNORED_CHROMOSOMES chrX -IGNORED_CHROMOSOMES chrY -IGNORED_CHROMOSOMES chrM INFO 2024-11-14 12:50:13 AssignCellsToSamples Genotype Quality [GQ] not found in header. Disabling GQ_THRESHOLD parameter |
Hi, thank you for submitting the log. You mentioned this is 10x data. In your command, you set the cell barcode tag correctly via It looks like you're very close - just set the MOLECULAR_BARCODE_TAG to UB and try again! |
When I use the --MOLECULAR_BARCODE_TAG UB, I am getting an error [akulan@cn2878 V1_outs_090624]$ /home/akulan/Drop-seq_tools-2.5.1/AssignCellsToSamples -m 24g -INPUT_BAM V1_possorted_genome_bam_tagged.bam -CELL_BC_FILE /data/GBMAD/McMahon-Claudia_SC_083024/V1_outs_090624/V1_filtered_feature_bc_matrix/barcodes.tsv.gz -VCF /data/akulan/gsa_may2024/PLINK_130524_0118/GSA_051324_scClaudiaV1labelled_chrAdded_snpsOnly.vcf.gz -OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.txt -VCF_OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.vcf.gz -CELL_BARCODE_TAG CB -MOLECULAR_BARCODE_TAG UB, -LOCUS_FUNCTION_LIST INTRONIC -IGNORED_CHROMOSOMES null -IGNORED_CHROMOSOMES chrX -IGNORED_CHROMOSOMES chrY -IGNORED_CHROMOSOMES chrM INFO 2024-11-12 20:16:20 AssignCellsToSamples Genotype Quality [GQ] not found in header. Disabling GQ_THRESHOLD parameter |
You set the tag name to be: The warning was because you set a BAM tag name to be 3 characters long, and bam tags are by definition length 2. |
It works! But, I get only one output V1_possorted_genome_bam_tagged.bam.donor_assignments.txt. I don't get the V1_possorted_genome_bam_tagged.bam.donor_assignments.vcf.gz. This is not there. Why is this missing? |
Log file please. |
Sorry, here it is. [akulan@cn2877 V1_outs_090624]$ /home/akulan/Drop-seq_tools-2.5.1/AssignCellsToSamples -m 24g -INPUT_BAM V1_possorted_genome_bam_tagged.bam -CELL_BC_FILE /data/GBMAD/McMahon-Claudia_SC_083024/V1_outs_090624/V1_filtered_feature_bc_matrix/barcodes.tsv.gz -VCF /data/akulan/gsa_may2024/PLINK_130524_0118/GSA_051324_scClaudiaV1labelled_chrAdded_snpsOnly.vcf.gz -OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.txt -VCF_OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.vcf.gz -CELL_BARCODE_TAG CB -MOLECULAR_BARCODE_TAG UB -LOCUS_FUNCTION_LIST INTRONIC -IGNORED_CHROMOSOMES null -IGNORED_CHROMOSOMES chrX -IGNORED_CHROMOSOMES chrY -IGNORED_CHROMOSOMES chrM INFO 2024-11-14 15:25:32 AssignCellsToSamples Genotype Quality [GQ] not found in header. Disabling GQ_THRESHOLD parameter |
Strange, everything in the log looks OK to me at first glance. Could you try using absolute paths instead of relative paths for the outputs? I don't have a strong opinion that will work, but we generally use all of our tools in workflows with absolute paths, so maybe there's a strange bug we never caught? |
Thanks Jim giving the absolute paths worked and did get the vcf output! |
I tried running AssignCellsToSamples with |
Hi,
We have 10x data from a cell village and I am trying to identify donors in the village using donor assignment. I have the 10x bam file, barcodes.tsv, gtf and vcf files.
Step 1: I ran Drop-seq_tools-2.5.1/TagReadWithGeneFunction on the bam file
Step 2: Drop-seq_tools-2.5.1/AssignCellsToSamples - outputs: V1_possorted_genome_bam_tagged.bam.donor_assignments.txt is empty (just header) and the *donor_assignments.vcf is not generated at all.
Here’s the command I am using:
/home/akulan/Drop-seq_tools-2.5.1/AssignCellsToSamples -m 24g -INPUT_BAM V1_possorted_genome_bam_tagged.bam -CELL_BC_FILE /data/GBMAD/McMahon-Claudia_SC_083024/V1_outs_090624/V1_filtered_feature_bc_matrix/barcodes.tsv.gz -VCF /data/akulan/gsa_may2024/PLINK_130524_0118/GSA_051324_scClaudiaV1labelled_chrAdded_snpsOnly.vcf.gz -OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.txt -VCF_OUTPUT V1_possorted_genome_bam_tagged.bam.donor_assignments.vcf.gz -CELL_BARCODE_TAG CB -LOCUS_FUNCTION_LIST INTRONIC -IGNORED_CHROMOSOMES null -IGNORED_CHROMOSOMES chrX -IGNORED_CHROMOSOMES chrY -IGNORED_CHROMOSOMES chrM
Any suggestions will be greatly appreciated.
Thank you very much.
Regards,
Nirmala
The text was updated successfully, but these errors were encountered: