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

What should a standard log of the running process look like? #6

Open
punching-samuel opened this issue Sep 14, 2024 · 2 comments
Open

Comments

@punching-samuel
Copy link

Hello,
I run bamquery

!python xxx/Software/bamquery/BamQuery/BamQuery.py xxx/Software/bamquery/BamQuery try_1 v38_104 --mode normal --t 8 --plots

And my output is as follows,

The output directory already exists in this path. 
BamQuery analysis will continue where it left.
Treatment File : Done!
Reverse Translation : Done!
Alignment : Done!
common_to_modes : Done!

It looks that the process is still running and hasn't finished yet, but my log file hasn't updated at all. My log as follows,

2024-09-14 14:17:39,174 INFO 
2024-09-14 14:17:39,174 INFO 
2024-09-14 14:17:39,174 INFO BamQuery analysis will continue where it left....
2024-09-14 14:17:39,174 INFO Path to input folder : xxx/Software/bamquery/BamQuery/ 
2024-09-14 14:17:39,175 INFO Path to output folder : xxx/Software/bamquery/output/ 
2024-09-14 14:17:39,175 INFO =============== # ===================
2024-09-14 14:17:39,175 INFO =============== Start Parameters ===============
2024-09-14 14:17:39,176 INFO  - BamQuery id : try_1 
2024-09-14 14:17:39,176 INFO  - Mode : normal, Strandedness :  False, Light:  False 
2024-09-14 14:17:39,176 INFO  - Single-Cell experiment (sc) :  False, Count UMIs : False
2024-09-14 14:17:39,176 INFO  - dbSNP :  0, COMMON SNPs : False, Genome Version : v38_104 
2024-09-14 14:17:39,176 INFO  - Plots : True
2024-09-14 14:17:39,176 INFO  - Keep Variant Alignments : False, Keep High Amount Alignments : False
2024-09-14 14:17:39,176 INFO  - Counting overlapping reads : False
2024-09-14 14:17:39,176 INFO  - Mouse Genome : False
2024-09-14 14:17:39,176 INFO  - Threads : 8
2024-09-14 14:17:39,176 INFO =============== End Parameters ===============
2024-09-14 14:17:39,192 INFO Total Bam Files to Query : 8.
2024-09-14 14:17:39,192 INFO Skipping peptide : Peptide because its length. Peptide should be between 8 and 11 aa.
2024-09-14 14:17:39,210 INFO Peptides to evaluate in Peptide Mode : 3617
2024-09-14 14:17:39,211 INFO Peptides to evaluate in Coding Sequence (CS) Mode : 0
2024-09-14 14:17:39,211 INFO Peptides to evaluate in Manual Mode: 0
2024-09-14 14:17:39,211 INFO Total Peptides to evaluate : 3617
2024-09-14 14:17:39,211 INFO ========== Treatment File : Done! ============ 
2024-09-14 14:17:39,214 INFO Fasta file with all the coding sequences already exists in the output path : xxx/Software/bamquery/output/genome_alignments/try_1.fastq --> Skipping this step!
2024-09-14 14:17:39,214 INFO ========== Reverse Translation : Done! ============ 
2024-09-14 14:17:39,401 INFO Alignment file already exists in the output folder : xxx/Software/bamquery/output/genome_alignments//Aligned.out.sam --> Skipping this step!
2024-09-14 14:17:39,401 INFO Alignment information already collected in the output folder : xxx/Software/bamquery/output/alignments//Alignments_information.dic --> Skipping this step!
2024-09-14 14:17:39,403 INFO Total alignments : 1107 
2024-09-14 14:17:39,403 INFO ========== Alignment : Done! ============ 
2024-09-14 14:17:39,403 INFO ========== Common_to_modes : Done! ============ 
2024-09-14 14:17:39,404 INFO Count information already collected in the output folder : xxx/Software/bamquery/output/res/temps_files/try_1_rna_count.csv --> Skipping this step!
2024-09-14 14:17:39,413 INFO ========== Get Count RNA : Done! ============ 

I will appreciate it if you can upload or show some standard log or the information about the intermediate stages.

@punching-samuel
Copy link
Author

punching-samuel commented Sep 14, 2024

And about the '/output/res/temps_files', I checked the try_1_rna_count.csv,

Peptide 2_STAR_sample-1 2_STAR_sample-2 2_STAR_sample-3 2_STAR_sample-4 2_STAR_sample-5 2_STAR_sample-6
1 ASLLDVFVLTR -2 -2 -2 -2 -2 -2
2 ATMELYQISQR -1 -1 -1 -1 -1 -1
3 AVALINAAIQK -1 -1 -1 -1 -1 -1
4 AVIQVSQIVAR -1 -1 -1 -1 -1 -1

All the values in the table are Negative numbers, could you assist me in understanding this situation? I'm looking for clarification and would greatly appreciate your guidance in explaining it further. Thank you in advance for your help.

Samuel.

@punching-samuel
Copy link
Author

I select only 50 peptides in peptides.tsv to run bamquery. The output seems correct.

The output directory already exists in this path. 
BamQuery analysis will continue where it left.
Treatment File : Done!
Reverse Translation : Done!
Alignment : Done!
common_to_modes : Done!
Getting counts for  8  samples
***** WARNING: File xxx/Software/bamquery/lib/genome_versions/genome_v38_104/gencode.v38.primary_assembly.annotation.gtf has inconsistent naming convention for record:
GL000009.2	ENSEMBL	gene	56140	58376	.	-	.	gene_id "ENSG00000278704.1"; gene_type "protein_coding"; gene_name "ENSG00000278704"; level 3;

***** WARNING: File xxx/Software/bamquery/lib/genome_versions/genome_v38_104/gencode.v38.primary_assembly.annotation.gtf has inconsistent naming convention for record:
GL000009.2	ENSEMBL	gene	56140	58376	.	-	.	gene_id "ENSG00000278704.1"; gene_type "protein_coding"; gene_name "ENSG00000278704"; level 3;

========== Genomic and ERE annotation summary by position : Done! ============ 
========== Genomic and ERE annotation by peptide : Done! ============ 
========== Transcription based biotyping by sample: Done! ============ 
========== Transcription based biotyping by group sample: Done! ============ 
========== BamQuery : Done! ============ 

But I still want to ask that if there sets the limited maximun numbers for peptides in peptides.tsv? My total 3617 peptides have total 1107 alignments, which means nearly 2/3 peptides not aligned? How to optimize the situation? Thank you in advance for your help.

With regards.
Samuel.

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

No branches or pull requests

1 participant