Skip to content
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

my proposed changes #185

Open
wants to merge 1 commit into
base: feature/cells_argument
Choose a base branch
from

Conversation

lldelisle
Copy link

No description provided.

@@ -19,11 +19,10 @@ install:
- conda config --add channels defaults
- conda config --add channels conda-forge
- conda config --add channels bioconda
- conda install -c bioconda -c conda-forge snakemake
- conda create -q -n snakemake snakemake>=5.3.1 python=3.6
- conda create -q -n snakemake 'snakemake>=5.3.1' 'python>=3.10'
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Requires python >= 3.10 to get the | None working.

Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I believe this is only for documentation updates, not the code itself.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Aaah. But then, which python is used? It was written python3.6 setup.py install...

Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I updated it on the branch to 3.11

@@ -57,7 +57,7 @@ def parse_barcode_reference(
barcode_pattern = rf"^[ATGC]{{{barcode_length}}}"

header = barcodes_pl.columns
set_dif = set(required_header) - set(header)
set_dif = set([required_header]) - set(header)
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

without the brackets it was splitting every letter of the required_header

Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Keep that one in, I'll merge once the PR is ready

@@ -74,7 +74,7 @@ def parse_barcode_reference(

else:
barcodes_pl = barcodes_pl.with_columns(
reference=pl.col(REFERENCE_COLUMN).str.strip_chars(STRIP_CHARS),
reference=pl.col(required_header).str.strip_chars(STRIP_CHARS),
Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

in barcodes_pl there is not the 'reference' column (at least in the whitelist barcode)

Copy link
Owner

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm changing the behavior here a little bit. I think I'm gonna use "reference" even in the whitelist input. Just makes the logic a bit cleaner

@lldelisle
Copy link
Author

To test I was doing:

$ python cite_seq_count/__main__.py -R1 tests/test_data/fastq/correct_R1_with_cell_barcode_mm.fastq.gz -R2 tests/test_data/fastq/correct_R2.fastq.gz -t tests/test_data/tags/pass/correct_3.csv -cbf 1 -cbl 16 -umif 17 -umil 26 -wl tests/test_data/whitelist.csv
Counting number of reads in file tests/test_data/fastq/correct_R1_with_cell_barcode_mm.fastq.gz
Writing chunks to disk
Reading reads from files: tests/test_data/fastq/correct_R1_with_cell_barcode_mm.fastq.gz, tests/test_data/fastq/correct_R2.fastq.gz
Mapping reads
Mapping done
Correcting barcodes
Barcodes corrected
Traceback (most recent call last):
  File "/home/ldelisle/Documents/mygit/CITE-seq-Count/cite_seq_count/__main__.py", line 208, in <module>
    main()
  File "/home/ldelisle/Documents/mygit/CITE-seq-Count/cite_seq_count/__main__.py", line 101, in main
    final_results=final_results,
UnboundLocalError: local variable 'final_results' referenced before assignment

@Hoohm
Copy link
Owner

Hoohm commented Nov 23, 2023 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants