Skip to content

Commit

Permalink
Merge pull request #676 from uclahs-cds/czhu-fix-rmats
Browse files Browse the repository at this point in the history
fix (moPepGen): rMATs -> rMATS
  • Loading branch information
zhuchcn authored Feb 8, 2023
2 parents 675f269 + d5779e8 commit ae7bf7e
Show file tree
Hide file tree
Showing 9 changed files with 18 additions and 18 deletions.
12 changes: 6 additions & 6 deletions moPepGen/cli/parse_rmats.py
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@

def add_rmats_input_arg(parser:argparse.ArgumentParser, name:str, message:str,
dest:str):
""" add input arg for rMATs """
""" add input arg for rMATS """
message += f" Valid formats: {INPUT_FILE_FORMATS}"
parser.add_argument(
name,
Expand All @@ -48,35 +48,35 @@ def add_subparser_parse_rmats(subparsers:argparse._SubParsersAction):
add_rmats_input_arg(
p, '--se',
message="File path to the SE (skipped exons) junction count file output"
" by rMATs. The file name should look like '*_SE.MATS.JC.txt' or "
" by rMATS. The file name should look like '*_SE.MATS.JC.txt' or "
"'*_SE.MATS.JCEC.txt'.",
dest='skipped_exon'
)
add_rmats_input_arg(
p, '--a5ss',
message="File path to the A5SS (alternative 5' splicint site) junction"
" count file output by rMATs. The file name should look like"
" count file output by rMATS. The file name should look like"
" '_S5SS.MATS.JC.txt' or '*_A5SS.MATS.JCEC.txt'.",
dest='alternative_5_splicing'
)
add_rmats_input_arg(
p, '--a3ss',
message="File path to the A3SS (alternative 3' splicint site) junction"
" count file output by rMATs. The file name should look like"
" count file output by rMATS. The file name should look like"
" '_S3SS.MATS.JC.txt' or '*_A3SS.MATS.JCEC.txt'.",
dest='alternative_3_splicing'
)
add_rmats_input_arg(
p, '--mxe',
message="File path to the MXE (mutually exclusive exons) junction"
" count file output by rMATs. The file name should look like"
" count file output by rMATS. The file name should look like"
" '_MXE.MATS.JC.txt' or '*_MXE.MATS.JCEC.txt'.",
dest='mutually_exclusive_exons'
)
add_rmats_input_arg(
p, '--ri',
message="File path to the RI (retained intron) junction"
" count file output by rMATs. The file name should look like"
" count file output by rMATS. The file name should look like"
" '_RI.MATS.JC.txt' or '*_RI.MATS.JCEC.txt'.",
dest='retained_intron'
)
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
##reference_index=/data/Index/GRCh38-EBI-GENCODE34
##genome_fasta=
##annotation_gtf=
##source=rMATs
##source=rMATS
##CHROM=<Description='{'Description': "'Gene ID'"}'>
##INFO=<ID=TRANSCRIPT_ID,Number=1,Type=String,Description="Transcript ID">
##INFO=<ID=GENE_SYMBOL,Number=1,Type=String,Description="Gene Symbol">
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
##reference_index=/data/Index/GRCh38-EBI-GENCODE34
##genome_fasta=
##annotation_gtf=
##source=rMATs
##source=rMATS
##CHROM=<Description='{'Description': "'Gene ID'"}'>
##INFO=<ID=TRANSCRIPT_ID,Number=1,Type=String,Description="Transcript ID">
##INFO=<ID=GENE_SYMBOL,Number=1,Type=String,Description="Gene Symbol">
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
##reference_index=/data/Index/GRCh38-EBI-GENCODE34
##genome_fasta=
##annotation_gtf=
##source=rMATs
##source=rMATS
##CHROM=<Description='{'Description': "'Gene ID'"}'>
##INFO=<ID=TRANSCRIPT_ID,Number=1,Type=String,Description="Transcript ID">
##INFO=<ID=GENE_SYMBOL,Number=1,Type=String,Description="Gene Symbol">
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
##reference_index=/data/Index/GRCh38-EBI-GENCODE34
##genome_fasta=
##annotation_gtf=
##source=rMATs
##source=rMATS
##CHROM=<Description='{'Description': "'Gene ID'"}'>
##INFO=<ID=TRANSCRIPT_ID,Number=1,Type=String,Description="Transcript ID">
##INFO=<ID=GENE_SYMBOL,Number=1,Type=String,Description="Gene Symbol">
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
##reference_index=/data/Index/GRCh38-EBI-GENCODE34
##genome_fasta=
##annotation_gtf=
##source=rMATs
##source=rMATS
##CHROM=<Description='{'Description': "'Gene ID'"}'>
##INFO=<ID=TRANSCRIPT_ID,Number=1,Type=String,Description="Transcript ID">
##INFO=<ID=GENE_SYMBOL,Number=1,Type=String,Description="Gene Symbol">
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
##reference_index=/data/Index/GRCh38-EBI-GENCODE34
##genome_fasta=
##annotation_gtf=
##source=rMATs
##source=rMATS
##CHROM=<Description='{'Description': "'Gene ID'"}'>
##INFO=<ID=TRANSCRIPT_ID,Number=1,Type=String,Description="Transcript ID">
##INFO=<ID=GENE_SYMBOL,Number=1,Type=String,Description="Gene Symbol">
Expand Down
10 changes: 5 additions & 5 deletions test/integration/test_call_variant_peptides.py
Original file line number Diff line number Diff line change
Expand Up @@ -484,7 +484,7 @@ def test_call_variant_peptide_case25(self):
def test_call_variant_peptide_case26(self):
""" Test case reported in #357 with two alternative splicing events. """
gvf = [
self.data_dir/'comb/CPCG0266_ENST00000381461.6/rMATs.gvf',
self.data_dir/'comb/CPCG0266_ENST00000381461.6/rMATS.gvf',
self.data_dir/'comb/CPCG0266_ENST00000381461.6/gsnp.gvf'
]
expected = self.data_dir/'comb/CPCG0266_ENST00000381461.6_expected.txt'
Expand All @@ -495,7 +495,7 @@ def test_call_variant_peptide_case27(self):
""" Test case reported in #360 with one alternative splicing and an
indel. """
gvf = [
self.data_dir/'comb/CPCG0235_ENST00000525687.5/rMATs.gvf',
self.data_dir/'comb/CPCG0235_ENST00000525687.5/rMATS.gvf',
self.data_dir/'comb/CPCG0235_ENST00000525687.5/gindel.gvf'
]
expected = self.data_dir/'comb/CPCG0235_ENST00000525687.5_expected.txt'
Expand All @@ -508,7 +508,7 @@ def test_call_variant_peptide_case28(self):
collapsing end nodes, the one with lower subgraph level should be kept.
"""
gvf = [
self.data_dir/'comb/CPCG0235_ENST00000590400.1/rMATs.gvf',
self.data_dir/'comb/CPCG0235_ENST00000590400.1/rMATS.gvf',
self.data_dir/'comb/CPCG0235_ENST00000590400.1/gsnp.gvf',
self.data_dir/'comb/CPCG0235_ENST00000590400.1/gindel.gvf'
]
Expand All @@ -522,7 +522,7 @@ def test_call_variant_peptide_case29(self):
aligning nodes in TVG, it failed to limit inside the same subgraph.
"""
gvf = [
self.data_dir/'comb/CPCG0333_ENST00000452737.5/rMATs.gvf',
self.data_dir/'comb/CPCG0333_ENST00000452737.5/rMATS.gvf',
self.data_dir/'comb/CPCG0333_ENST00000452737.5/gsnp.gvf',
self.data_dir/'comb/CPCG0333_ENST00000452737.5/gindel.gvf'
]
Expand All @@ -536,7 +536,7 @@ def test_call_variant_peptide_case30(self):
aligning nodes in TVG, it failed to limit inside the same subgraph.
"""
gvf = [
self.data_dir/'comb/CPCG0462_ENST00000483923.5/rMATs.gvf',
self.data_dir/'comb/CPCG0462_ENST00000483923.5/rMATS.gvf',
self.data_dir/'comb/CPCG0462_ENST00000483923.5/gsnp.gvf'
]
expected = self.data_dir/'comb/CPCG0462_ENST00000483923.5_expected.txt'
Expand Down
2 changes: 1 addition & 1 deletion test/integration/test_parse_rmats.py
Original file line number Diff line number Diff line change
Expand Up @@ -43,7 +43,7 @@ def test_parse_rmats_cli(self):
-o {self.work_dir}/reditools.gvf \\
-g {self.data_dir}/genome.fasta \\
-a {self.data_dir}/annotation.gtf \\
--source rMATs
--source rMATS
"""
res = sp.run(cmd, shell=True, check=False, capture_output=True)
try:
Expand Down

0 comments on commit ae7bf7e

Please sign in to comment.