@@ -392,6 +392,7 @@ def remap_symbols(row, gene_positiondf):
392392 row ["Hugo_Symbol" ] = symbol
393393 return row
394394
395+
395396class bed (FileTypeFormat ):
396397 """GENIE bed format"""
397398
@@ -457,10 +458,10 @@ def create_gene_panel(self, beddf, seq_assay_id, gene_panel_path, parentid):
457458 if not beddf .empty :
458459 # exonsdf = beddf[beddf["Feature_Type"] == "exon"]
459460 # Only include genes that should be included in the panels
460- #include_exonsdf = exonsdf[exonsdf["includeInPanel"]]
461-
461+ # include_exonsdf = exonsdf[exonsdf["includeInPanel"]]
462+
462463 # Include all genes in the panels
463- include_exonsdf = beddf [beddf ["Feature_Type" ] == "exon" ]
464+ include_exonsdf = beddf [beddf ["Feature_Type" ] == "exon" ]
464465 # Write gene panel
465466 null_genes = include_exonsdf ["Hugo_Symbol" ].isnull ()
466467 unique_genes = set (include_exonsdf ["Hugo_Symbol" ][~ null_genes ])
@@ -515,10 +516,10 @@ def _process(self, beddf, seq_assay_id, newpath, parentid, create_panel=True):
515516 seq_assay_id = seq_assay_id .upper ()
516517 seq_assay_id = seq_assay_id .replace ("_" , "-" )
517518
518- # add in placeholder includeInPanel column and value if it doesn't exits
519+ # add in placeholder includeInPanel column and value if it doesn't exist
519520 if len (beddf .columns ) == 4 :
520521 beddf [4 ] = True
521-
522+
522523 # Add in 6th column which is the clinicalReported
523524 if len (beddf .columns ) > 5 :
524525 if all (beddf [5 ].apply (lambda x : x in [True , False ])):
@@ -535,7 +536,7 @@ def _process(self, beddf, seq_assay_id, newpath, parentid, create_panel=True):
535536 exon_gtf_path , gene_gtf_path = create_gtf (process_functions .SCRIPT_DIR )
536537 LOGGER .info ("REMAPPING {}" .format (seq_assay_id ))
537538 # bedname = seq_assay_id + ".bed"
538-
539+
539540 beddf .columns = [
540541 "Chromosome" ,
541542 "Start_Position" ,
@@ -638,7 +639,7 @@ def _validate(self, beddf):
638639 "Start_Position" ,
639640 "End_Position" ,
640641 "Hugo_Symbol" ,
641- #"includeInPanel",
642+ # "includeInPanel",
642643 ]
643644 if len (beddf .columns ) < len (newcols ):
644645 total_error += (
@@ -682,15 +683,15 @@ def _validate(self, beddf):
682683 "Hugo_Symbol column, not the strand column\n "
683684 )
684685
685- #warn, error = process_functions.check_col_and_values(
686+ # warn, error = process_functions.check_col_and_values(
686687 # beddf,
687688 # "includeInPanel",
688689 # [True, False],
689690 # filename="BED file",
690691 # required=True,
691- #)
692- #warning += warn
693- #total_error += error
692+ # )
693+ # warning += warn
694+ # total_error += error
694695
695696 if to_validate_symbol :
696697 gene_position_table = self .syn .tableQuery ("SELECT * FROM syn11806563" )
0 commit comments