From 5344849b0cf39cea2ad049b5177590473b14e9d7 Mon Sep 17 00:00:00 2001 From: Taylor Reiter Date: Wed, 5 Jun 2024 10:47:37 -0400 Subject: [PATCH 1/5] remove short and longer contigs input --- README.md | 10 +- Snakefile | 59 ++------ config.yml | 31 ++-- demo/README.md | 6 +- demo/config.yml | 32 ++--- ..._than_r2t_minimum_length.fa => contigs.fa} | 135 ++++++++++++++++++ ...contigs_shorter_than_r2t_minimum_length.fa | 135 ------------------ 7 files changed, 185 insertions(+), 223 deletions(-) rename demo/{contigs_longer_than_r2t_minimum_length.fa => contigs.fa} (95%) delete mode 100644 demo/contigs_shorter_than_r2t_minimum_length.fa diff --git a/README.md b/README.md index 98cf989..f3d1714 100644 --- a/README.md +++ b/README.md @@ -34,11 +34,11 @@ snakemake --software-deployment-method conda -j 1 --configfile demo/config.yml ## Input data -The [peptigate pipeline](./Snakefile) requires two pairs of input files (for a total of four input files): -* A transcriptome assembly, split into "long" and "short" contigs. - * Long transcripts/contigs: A transcriptome assembly FASTA file in nucleotide format containing transcripts or contigs longer than X nucleotides (typically, 300-500nt). - * Short transcripts/contigs: contigs that are shorter than X nucleotides (typically, 300-500nt). Some transcriptome assemblers discard short contigs and do not include them in the final assembly. However, some provide them as an intermediate output file. These contigs may contain sORFs and so are included as an input to the peptigate pipeline. If you do not have a file that contains very short contigs, provide a path to an empty file. Short contigs can also be provided as part of the previous file. If that is the case, provide a path to an empty file for this input file. -* Open reading frames predicted from the transcriptome in both amino acid and nucleotide format. The open reading frames in both files should have the same names. Tools like [Transdecoder](https://github.com/TransDecoder/TransDecoder) provide these files in the correct format. +The [peptigate pipeline](./Snakefile) requires three input files: +* A transcriptome assembly: transcriptome assembly FASTA file in nucleotide format containing transcripts or contigs. +* Open reading frames predicted from the transcriptome in both amino acid and nucleotide format. + The open reading frames in both files should have the same names before the first period in the FASTA header name. + Tools like [TransDecoder](https://github.com/TransDecoder/TransDecoder) provide these files in the correct format. * Open reading frames in amino acid format: A FASTA file of predicted open reading frames in amino acid format. * Open reading frames in nucleotide format: A FASTA file of predicted open reading frames in nucleotide format. diff --git a/Snakefile b/Snakefile index 5262839..9470a88 100644 --- a/Snakefile +++ b/Snakefile @@ -21,8 +21,7 @@ OUTPUT_DIR = Path(config["output_dir"]) ORFS_AMINO_ACIDS = Path(config["orfs_amino_acids"]) ORFS_NUCLEOTIDES = Path(config["orfs_nucleotides"]) -CONTIGS_SHORTER = Path(config["contigs_shorter_than_r2t_minimum_length"]) -CONTIGS_LONGER = Path(config["contigs_longer_than_r2t_minimum_length"]) +CONTIGS = Path(config["contigs"]) PLMUTILS_MODEL_DIR = Path(config["plmutils_model_dir"]) ################################################################################ @@ -48,33 +47,6 @@ Note that we follow the conventions of the prokka tool for output file suffixes ## sORF prediction ################################################################################ - -rule combine_contigs: - """ - By default we assume that files provided to this pipeline are reads2transcriptome outputs. - Reads2transcriptome outputs two files that contain contigs. - The first we refer to as contigs_shorter_than_r2t_minimum (or CONTIGS_SHORTER), - which are transcripts output by assemblers that did not meet the r2t runs minimum contig length. - The second we refer to as contigs_longer_than_r2t_minimum (or CONTIGS_LONGER) and are assembled - transcripts that passed the isoform clustering and decontamination steps of r2t. - We no longer need these pools of transcripts differentiated, so we combine them in this rule. - If your input transcriptome only has one file of assembled contigs, sequences only need to be - supplied in contigs_longer_than_r2t_minimum_length. - contigs_shorter_than_r2t_minimum_length can be an empty file. - """ - input: - contigs_shorter=CONTIGS_SHORTER, - contigs_longer=CONTIGS_LONGER, - output: - all_contigs=OUTPUT_DIR / "sORF" / "contigs" / "all_input_contigs.fna", - conda: - "envs/seqkit.yml" - shell: - """ - cat {input.contigs_shorter} {input.contigs_longer} > {output.all_contigs} - """ - - rule get_coding_contig_names: """ Extract amino acid contig names and remove everything after the first period, @@ -95,24 +67,21 @@ rule get_coding_contig_names: rule filter_contigs_to_no_predicted_ORF: """ - The peptigate pipeline takes as input transcripts (provided in two files) and predicted or - annotated coding genes (provided in nucleotide and amino acid formats). - This rule removes contigs with coding genes from the full transcript file. - It assumes that the contig names are the same between to the two files (everything before the - first period) and uses an inverted grep on the sequence names in the coding file to - eliminate transcripts that contain protein-coding genes. - It keeps all other transcripts, regardless of length, to investigate the presence of an sORF - later in the pipeline. - - We expect that read2transcriptome will often be used to used to create input files for - peptigate, or that transdecoder will be used to predict which transcripts contain protein-coding - genes (the r2t pipeline also uses transdecoder to predict open reading frames (ORFs) from - transcripts). By default, only ORFs that are longer than 100 amino acids are kept by - transdecoder. The peptigate pipeline predicts peptides that are 100 amino acids or shorter. - This rule eliminates transcripts that contained a transdecoder-predicted ORF. + The peptigate pipeline takes as input transcripts and predicted or annotated coding genes + (provided in nucleotide and amino acid formats). This rule removes contigs with coding genes + from the transcriptome assembly file (contigs). It assumes that the contig names are the same in + the transcriptome assembly file and the annotated coding genes file (everything before the first + period) and uses an inverted grep on the sequence names in the coding file to eliminate + transcripts that contain protein-coding genes. It keeps all other transcripts, regardless of + length, to investigate the presence of an sORF later in the pipeline. + + We expect that transdecoder will often be used to predict which transcripts contain + protein-coding genes (ORF prediction). By default, only ORFs that are longer than 100 amino + acids are kept by transdecoder. The peptigate pipeline predicts peptides that are 100 amino + acids or shorter. This rule eliminates transcripts that contained a transdecoder-predicted ORF. """ input: - fna=rules.combine_contigs.output.all_contigs, + fna=CONTIGS, names=rules.get_coding_contig_names.output.names, output: fna=OUTPUT_DIR / "sORF" / "contigs" / "contigs_with_no_annotated_orf.fna", diff --git a/config.yml b/config.yml index 3b38af2..7f2124c 100644 --- a/config.yml +++ b/config.yml @@ -14,25 +14,22 @@ input_dir: "path/to/inputdir/" output_dir: "path/to/outputdir" -# All input files are produced by reads2transcriptome. -# While this is not a strict requirement, using these output files gives us the ability to look at very short contiguous sequences. +# Input files and directories +# - contigs: transcriptome assembly contigs. # - orfs_amino_acids: predicted ORFs translated into amino acids. -# Output by transdecoder. -# Used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases. -# - orfs_nucleotides: predicted ORFs as nucleotide sequences. Output by transdecoder. -# Used to compare peptide nucleotide sequences (clustering, dn/ds estimation, etc.). -# - contigs_shorter_than_r2t_minimum_length: contigs that are shorter than X nucleotides (by default 75bp). -# The reads2transcriptome pipeline assembles RNA-seq reads into contigs (transcripts) using multiple assemblers and then merges those assemblies together. -# Before merging, very short contigs are removed (<75bp). -# However, reads2transcriptome outputs a FASTA file containing these transcripts, which is used as input to peptigate here. -# These contigs may contain sORFs and so are included as an input to the peptigate pipeline. -# - contigs_longer_than_r2t_minimum_length: contigs longer than X nucleotides (default is 75bp). -# If the user did not use the reads2transcriptome file to generate their input files, this would be a transcriptome assembly FASTA file in nucleotide format containing transcripts. -# Note that the first step of sORF prediction combines the contigs_shorter_than_r2t_minimum_length and contigs_longer_than_r2t_minimum_length files, so there is no need to perform any pre-processing by length. -# - plmutils_model_dir: path to the directory for the plmutils model that will predict whether sORFs are coding or non-coding. +# ORFs should be predicted from the same transcriptome assembly as the "contigs" input file. +# ORFs should have the same name (before the first period in the name) as the contigs in the +# "contigs" input file. TransDecoder provides files in the proper format. +# Used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases, and to +# remove coding transcripts from the transcriptome assembly before sORF prediction. +# - orfs_nucleotides: predicted ORFs as nucleotide sequences. Should contain the same ORFs as +# "orfs_amino_acids" but in nucleotide format. TransDecoder also provides this file in the proper +# format. If this file contains short ORFs (< 300 nucleotides), they will not be reported as sORFs +# as they are already annotated in the input. +# - plmutils_model_dir: path to the directory for the plmutils model that will predict whether sORFs +# are coding or non-coding. +contigs: "path/to/contigs.fa" orfs_amino_acids: "path/to/orfs.faa" orfs_nucleotides: "path/to/orfs.fa" -contigs_shorter_than_r2t_minimum_length: "path/to/short_contigs.fa" -contigs_longer_than_r2t_minimum_length: "path/to/long_contigs.fa" plmutils_model_dir: "inputs/models/plmutils/" diff --git a/demo/README.md b/demo/README.md index c7e29e7..8580f67 100644 --- a/demo/README.md +++ b/demo/README.md @@ -4,9 +4,9 @@ curl -JLo Amblyomma_americanum_transcriptome_assembly_data.tar.gz https://zenodo tar xf Amblyomma_americanum_transcriptome_assembly_data.tar.gz mv transcriptome_data/* . head -n 200 orthofuser_final_clean.fa.transdecoder.pep > orfs_amino_acids.faa -head -n 200 orthofuser_final_clean.fa.dammit.fasta > all_contigs.fa +head -n 200 orthofuser_final_clean.fa.dammit.fasta > contigs.fa head -n 204 orthofuser_final_clean.fa.transdecoder.cds > orfs_nucleotides.fa ``` -We pulled the "short contigs" file from an internal S3 bucket. -It contains contigs that were filtered from the Amblyomma transcriptome prior to txome merging. +We also pulled short contigs (less than 75 bp) from an internal S3 bucket and added these contigs to the `contigs.fa` file (50 contigs). +These are contigs that were filtered from the *Amblyomma* transcriptome prior to transcriptome merging (mid assembly pipeline). diff --git a/demo/config.yml b/demo/config.yml index bb82214..1fa5d27 100644 --- a/demo/config.yml +++ b/demo/config.yml @@ -14,26 +14,22 @@ input_dir: "inputs/" output_dir: "outputs/demo" -# All input files are produced by reads2transcriptome. -# While this is not a strict requirement, using these output files gives us the ability to look at very short contiguous sequences. -# TER TODO: update names of output files based on Nextflow of reads2transcriptome. +# Input files and directories +# - contigs: transcriptome assembly contigs. # - orfs_amino_acids: predicted ORFs translated into amino acids. -# Output by transdecoder. -# Used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases. -# - orfs_nucleotides: predicted ORFs as nucleotide sequences. Output by transdecoder. -# Used to compare peptide nucleotide sequences (clustering, dn/ds estimation, etc.). -# - contigs_shorter_than_r2t_minimum_length: contigs that are shorter than X nucleotides (by default 75bp). -# The reads2transcriptome pipeline assembles RNA-seq reads into contigs (transcripts) using multiple assemblers and then merges those assemblies together. -# Before merging, very short contigs are removed (<75bp). -# However, reads2transcriptome outputs a FASTA file containing these transcripts, which is used as input to peptigate here. -# These contigs may contain sORFs and so are included as an input to the peptigate pipeline. -# - contigs_longer_than_r2t_minimum_length: contigs longer than X nucleotides (default is 75bp). -# If the user did not use the reads2transcriptome file to generate their input files, this would be a transcriptome assembly FASTA file in nucleotide format containing transcripts. -# Note that the first step of sORF prediction combines the contigs_shorter_than_r2t_minimum_length and contigs_longer_than_r2t_minimum_length files, so there is no need to perform any pre-processing by length. -# - plmutils_model_dir: path to the directory for the plmutils model that will predict whether sORFs are coding or non-coding. +# ORFs should be predicted from the same transcriptome assembly as the "contigs" input file. +# ORFs should have the same name (before the first period in the name) as the contigs in the +# "contigs" input file. TransDecoder provides files in the proper format. +# Used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases, and to +# remove coding transcripts from the transcriptome assembly before sORF prediction. +# - orfs_nucleotides: predicted ORFs as nucleotide sequences. Should contain the same ORFs as +# "orfs_amino_acids" but in nucleotide format. TransDecoder also provides this file in the proper +# format. If this file contains short ORFs (< 300 nucleotides), they will not be reported as sORFs +# as they are already annotated in the input. +# - plmutils_model_dir: path to the directory for the plmutils model that will predict whether sORFs +# are coding or non-coding. +contigs: "demo/contigs.fa" orfs_amino_acids: "demo/orfs_amino_acids.faa" orfs_nucleotides: "demo/orfs_nucleotides.fa" -contigs_shorter_than_r2t_minimum_length: "demo/contigs_shorter_than_r2t_minimum_length.fa" -contigs_longer_than_r2t_minimum_length: "demo/contigs_longer_than_r2t_minimum_length.fa" plmutils_model_dir: "inputs/models/plmutils/" diff --git a/demo/contigs_longer_than_r2t_minimum_length.fa b/demo/contigs.fa similarity index 95% rename from demo/contigs_longer_than_r2t_minimum_length.fa rename to demo/contigs.fa index 4156807..6496905 100644 --- a/demo/contigs_longer_than_r2t_minimum_length.fa +++ b/demo/contigs.fa @@ -198,3 +198,138 @@ GGAAAGCTAGGCCAAGAGCCAAGCCACATCAAAATTCAAAGTGGCTTCTAGAACTCGTTGCCACAGTTGTGATGTAAAGA GTCGGGTTCGAACCCGGGAACTCCGGATCAGTAGTTGGACGCTCCAACCACTGAGCCACCGCGGCGGGGATCGCACGTGGATAATCAATTTCTCCGCCGCGAAGCGCTACCGAAGGTATGCAGACGCACACGTCACTCTTTGGTAGGATGTTGAAGGCTTCTTCTATCTCCCTGGTCCGACCAAAATATGACGGGATGACCACCGTGTTGACAAGCTGCGGAGCGCACTCATGCACCGTCCGATGCCCCGCTAATTGACCGCTGATAGCGCACTCCAGGGATCCCCCGTGCTCTCGAAGCCGTTCATTCAAACAGCGCTCCGACCAATGTAGGTGCGTCCGCATGACAAAGGTATTCCACAAACCCCGCTTTTCACACATTTCACGCCCTCTCCGGTGTTCCCTCCTAAACTGATTCTTTT >Transcript_99 len=421 CDS=1-421 exon=0-421 gene=0-421 mRNA=0-421 CAAGGACAAGGATGGCAAGAACTTGCTGAAGGTGGTTATGAGGACATGGCTTCCTGCTGGAGATGCACTCTTCGAAATGATCACAATCCACCTTCCATCTCCTGTGACGGCACAACGCTACCGAATGGAAATCCTCTATGAGGGTCCCCTTGATGATGAGGCTGCTGTTGCTGTAAAGGCGTGCGACCCAGAAGGCCCACTGATGATGTACGTCTCAAAGATGGTACCGACATCGGACAAGGGCCGCTTTTACGCATTTGGACGAGTCTTCTCTGGTGTTGTGTCATCGGGACAAAAGGTCCGCATCATGGGCCCAAACTACACTCCTGGCAAGAAAGAGGACTTGGCTGAGAAGGCCATTCAGAGGACTGTTCTGATGATGGGTCGGTACGTGGAGCCCATTGAAGATGTGCCCTGCGGT +>peaae2cellline_NODE_191792_length_57_cov_5.882979_g176311_i0 +AAAATACGAAATGGTGCTTTTTGGTTTGCAATAAACACAGTTAGGTGTGCCTGTTTC +>peaae2cellline_NODE_191793_length_57_cov_2.702128_g176312_i0 +GGATTCGAAGTCGGGCAAAGAGCAGGGAACATTCCACGCTCTCAGTTATGGTGTCAT +>peaae2cellline_NODE_191799_length_47_cov_1.500000_g176318_i0 +AATAGCTGTCGAATTCTGCCACTTTCATTATGCCATCTTCGACGGGA +>peaae2cellline_NODE_191803_length_46_cov_1.710843_g176322_i0 +TACAGCACCACAGCACTTTTGCACCGTGGGTTTGCGCCTTGCCTGC +>peaae2cellline_NODE_191804_length_46_cov_1.228916_g176323_i0 +CTTTGAGAAGTAATTAAATTACATTACAAATTACTCTTCTCAAAAA +>peaae2cellline_NODE_191809_length_38_cov_1.253333_g176328_i0 +CCTTGTCTTGAAAATGAAAACCGGTTTTGGGGAAAAGA +>peaae1cellline_NODE_200257_length_72_cov_1.972477_g183439_i0 +TGTTCCTGTGGTGCTGGGTCTTTGCAACTGAAATAAATAAACAAGCCTAGACTGACCCAG +CACCACAGGAAC +>peaae1cellline_NODE_200261_length_59_cov_1.343750_g183443_i0 +AGGCAGTTACTTACTGCTTTGGCAGGACGGTCGAAGGGGCAAGGTAGACAATGTAGTGG +>peaae1cellline_NODE_200262_length_58_cov_2.178947_g183444_i0 +CTGTACAACGCACTGTTTTCCTCTGTTGCTAACTATTGTATTTTGGTCTGGGGTACAA +>peaae1cellline_NODE_200264_length_56_cov_1.892473_g183446_i0 +GAGAGACCACACTTCTCATTCATGTGGCACCAGCACAGAAAAACGTGCAAGATCTA +>peaae1cellline_NODE_200266_length_52_cov_1.393258_g183448_i0 +TCTAATTCTGTGCTTCGAGGACGATTTCTTGTCTAAGATGCACAGAATTAGG +>peaae1cellline_NODE_200272_length_44_cov_5.555556_g183454_i0 +GGCCATGAACAACCAGTCCTGTCCTGTCTTGTTTCTTCCTGGTA +>peaae1cellline_NODE_200273_length_42_cov_44.962025_g183455_i0 +TACCGAGCTGCGAGGCGGCCTGTATAAAAACGCTGTCAGTGT +>petxwholemale_NODE_12966_length_61_cov_3.076923_g12631_i0 +CCTACACTTTCAAACCTGAAAGGTTTCTGCCTGAGAGCAAAGACCTCCTGAAGCCGTTCT +C +>petxwholefemale_NODE_12365_length_47_cov_1.677778_g12066_i0 +ATATTTCGAATCAGATAACGCGATGACGGCCGTCAGCTTCCATCAGC +>petx24wholefemale_NODE_12274_length_66_cov_2.137615_g11799_i0 +CTACCGAGAGTCATTTTCTTATGCCGTTCTAATACTTATAGTAGGCTATAGAAAAAATAT +CTGTTC +>petx0wholemale_NODE_1133_length_67_cov_2.293103_g1132_i0 +GCAAGGTGTCTTACGGGTGGTGATGCCCACCAAGACCCTGCGCATCACCACCCGTAAGAC +ACCTTGT +>petx0wholefemale_NODE_618915_length_74_cov_1.682927_g519641_i0 +AGGCTGAACCAAGACACCATGACGACTTGCCATCTGCAACATGGCTGAATGTAAAATGTA +AATTGGTTCAGCCC +>petx0wholefemale_NODE_618917_length_73_cov_2.401639_g519643_i0 +TAGTATAAAAACATGTCAAAGTGGCAAGGTAAACTGTGTAAACATGGCAGTAGCACTCAA +GGTTGGCATTTAT +>petx0wholefemale_NODE_618918_length_72_cov_3.909091_g519644_i0 +TTCGCCTAACAGCATTAGCAGCAGCCAAGAAGCAAATTGCCAGCAACAAGGTTGGTAAAA +GCAGGGGGCGAC +>petx0wholefemale_NODE_618919_length_72_cov_3.148760_g519645_i0 +TAATTGCGCGACGTTTGTCTCTCGACGGGATTGCACCAACTTTCACTGAAGGGCTTTCAC +GTGGCGTACTTC +>petx0wholefemale_NODE_618920_length_72_cov_3.090909_g519646_i0 +TTTAGCTTGTTTGGTTAGTTGTGTCTTTTGGTCGTTTGGCTCGCAAAGTCAGAACTAACC +AAACAAGCGAAC +>petx0wholefemale_NODE_618921_length_72_cov_2.272727_g519647_i0 +CCCAACCAGGGGAAATCAGCAGTTGCTTTTTCCTGTCTCCCCCTCCATCTTTCACTTTCA +TATCTCTCCTTC +>petx0wholefemale_NODE_618923_length_71_cov_12.450000_g519649_i0 +TATTCTGCTGACATAATAAAAAACTGTCAAATATGGGCCCTCGCACAGTTTTTTATTATG +ACAGCAGAATC +>petx0wholefemale_NODE_618925_length_71_cov_3.791667_g519651_i0 +TCGTGTTTTCAGTTCATAACCCACCTGTGGCAAGCATGGTGTTTATGATGGGTTATGAAC +TGAATACACGT +>petx0wholefemale_NODE_618926_length_70_cov_30.613445_g519652_i0 +TTAACATATTATTGAAACGGGACACTGATAACGCAGCCTTGATCTGTGTCCCGTTTCAAT +AATATGTTAT +>petx0wholefemale_NODE_618927_length_70_cov_1.722689_g519653_i0 +AACAGTTTTGACAAAGCTGCTCGCAAAACCAGCTGATGGCGATATTACCTGCATTTTATT +TTTTTTTACC +>petx0wholefemale_NODE_618928_length_70_cov_1.512605_g519654_i0 +ACGCACCTTTGGTCAATTTCATCCTATCGCTGTCACACGACGACACGACAGTATTTCTAG +AATTCTGCAG +>petx0wholefemale_NODE_618929_length_69_cov_4.567797_g519655_i0 +TAGAATGTTATTCAGGTCCTCCACTACATCCCCCAGCCGGATGTAGTGGAGGACCTGAAT +AACATTCTT +>petx0wholefemale_NODE_618931_length_68_cov_8.094017_g519657_i0 +GAACGAGTGCAATTCGCTTGTCCCATTCGATGGGGCTTGGAAATGGGCCAAGCGAATTGC +ACTCGTTA +>petx0wholefemale_NODE_618932_length_68_cov_3.153846_g519658_i0 +GTGGCAGTTGCCGCCATGCTCGGCTTTTTTTTTATCTTTGTTTAACTGTGTGGTAAATTA +TTAAAATT +>petx0wholefemale_NODE_618933_length_68_cov_2.717949_g519659_i0 +GCGGTCTGGCCGGTCCGTATTTATTTTTCATCTTTGGAGCCGTTATCGTTTTTTTTTTTT +CACCCGGT +>petx0wholefemale_NODE_618934_length_68_cov_2.282051_g519660_i0 +ATAAAACCGGGGCGAGGCATGACCCAATGCATATCTCCTTCTTTTGCAGATAAAACTCCT +TGTCAAAG +>petx0wholefemale_NODE_618935_length_68_cov_1.145299_g519661_i0 +CAAGTTACCACAATACAATGCTCAAAAAATAAAGGGTTGCATTTTGAGCATTGTATTGGG +GTAACGTT +>petx0wholefemale_NODE_618937_length_67_cov_10.827586_g519663_i0 +AGACATCTGCTCTATGTGTGCTGTCATACTGACTATTCATATGGCAGCACACATAGAGCC +GAGGTCG +>petx0wholefemale_NODE_618939_length_67_cov_1.206897_g519665_i0 +CCTATACTGCTTAAGAGCAGTCAGGGCTGGATATCTCACTGTTGGGCCTGCTCTTAAGCA +GTATAGA +>petx0wholefemale_NODE_618940_length_67_cov_1.172414_g519666_i0 +CTTGCCTTTCTCGGGTAGAGGTATCGTCCCTGGTTATAAAGCCAGCAACAGCAAAAGAAA +TGCAGTA +>petx0wholefemale_NODE_618941_length_66_cov_1.008696_g519667_i0 +CCTTTAATTTTAAGACTCACTTCCGCTCCCTCGAGCACAGCAGCAGGTGAGTCTTAAAAT +TAAAGT +>petx0wholefemale_NODE_618944_length_65_cov_5.570175_g519670_i0 +GCCGTCCATCTCTATACACGACACAGACATGTGTGTGGCCGCGTTTCTGTCGCTGAGTTG +CTGAG +>petx0wholefemale_NODE_618945_length_65_cov_4.210526_g519671_i0 +GAAGGACGCATTAACGGAGCTCTCTTGGATTCACTGTCCGCTGGTGGCATTGATTACAGC +GTCCA +>petx0wholefemale_NODE_618947_length_65_cov_2.035088_g519673_i0 +GCTTCATGCAGAATGAAGAGAAACCGGAACGGATGGCAGCAGGTAGCATTTTGTCTGCAT +GAAGG +>petx0wholefemale_NODE_618948_length_65_cov_2.035088_g519674_i0 +TTCTGTCTCGCTTACGCCATTCTCGCCACGTTGGCAGCTGTCAAGACTCTCGCTTCGCCA +GTTTT +>petx0wholefemale_NODE_618949_length_63_cov_1.366071_g519675_i0 +CAAAAGCGTACAACTGTCAGCGACGAGCCACGAGCGCCCTGCGCTGACAGTTGTACGCTT +TTT +>petx0wholefemale_NODE_618951_length_61_cov_38.827273_g519677_i0 +TTTATTCCACATATGTAAAATCAAGATGAACAAAAGTGTGGTGCTTTTCAGCTGGAAATA +T +>petx0wholefemale_NODE_618952_length_61_cov_1.945455_g519678_i0 +GCCATCTTAACGAGCTATTATTATTTTAATATTATAAATAATAATCGCGCGATAAGATGG +G +>petx0wholefemale_NODE_618953_length_61_cov_1.845455_g519679_i0 +GCTGAAAGGTCAGCAGCCGAGCGCTATAACCATTGCGCCACCACGGCGCCAACCTTTCAG +A +>petx0wholefemale_NODE_618954_length_61_cov_1.200000_g519680_i0 +AAAGATCCCCATCAAAGGTGGTTTTTGTCTCTCTTAAAAAGCCGCTTTGTTGGGGATGTT +G +>petx0wholefemale_NODE_618956_length_60_cov_1.394495_g519682_i0 +CTTCAAACAATAACTGCAGTTCGCGATTCACTGGTAACTCCCATGTAAACGTGAGGCCCT +>petx0wholefemale_NODE_618957_length_60_cov_1.229358_g519683_i0 +AGCTTTATTAGCCGTGCATAAAACGGAAGTCTTTCCGCATCCAGCAGGGCCAGTGAAGAA +>petx12wholefemaleecoli_NODE_281510_length_74_cov_17.659864_g256556_i0 +CTTTCACGTCCCTGTACATAGCCCTTCCTTTGTTTATGGCTTGTTAAATAAAAGCTGGTT +TACCTTGTTGAAAA diff --git a/demo/contigs_shorter_than_r2t_minimum_length.fa b/demo/contigs_shorter_than_r2t_minimum_length.fa deleted file mode 100644 index da0f38d..0000000 --- a/demo/contigs_shorter_than_r2t_minimum_length.fa +++ /dev/null @@ -1,135 +0,0 @@ ->peaae2cellline_NODE_191792_length_57_cov_5.882979_g176311_i0 -AAAATACGAAATGGTGCTTTTTGGTTTGCAATAAACACAGTTAGGTGTGCCTGTTTC ->peaae2cellline_NODE_191793_length_57_cov_2.702128_g176312_i0 -GGATTCGAAGTCGGGCAAAGAGCAGGGAACATTCCACGCTCTCAGTTATGGTGTCAT ->peaae2cellline_NODE_191799_length_47_cov_1.500000_g176318_i0 -AATAGCTGTCGAATTCTGCCACTTTCATTATGCCATCTTCGACGGGA ->peaae2cellline_NODE_191803_length_46_cov_1.710843_g176322_i0 -TACAGCACCACAGCACTTTTGCACCGTGGGTTTGCGCCTTGCCTGC ->peaae2cellline_NODE_191804_length_46_cov_1.228916_g176323_i0 -CTTTGAGAAGTAATTAAATTACATTACAAATTACTCTTCTCAAAAA ->peaae2cellline_NODE_191809_length_38_cov_1.253333_g176328_i0 -CCTTGTCTTGAAAATGAAAACCGGTTTTGGGGAAAAGA ->peaae1cellline_NODE_200257_length_72_cov_1.972477_g183439_i0 -TGTTCCTGTGGTGCTGGGTCTTTGCAACTGAAATAAATAAACAAGCCTAGACTGACCCAG -CACCACAGGAAC ->peaae1cellline_NODE_200261_length_59_cov_1.343750_g183443_i0 -AGGCAGTTACTTACTGCTTTGGCAGGACGGTCGAAGGGGCAAGGTAGACAATGTAGTGG ->peaae1cellline_NODE_200262_length_58_cov_2.178947_g183444_i0 -CTGTACAACGCACTGTTTTCCTCTGTTGCTAACTATTGTATTTTGGTCTGGGGTACAA ->peaae1cellline_NODE_200264_length_56_cov_1.892473_g183446_i0 -GAGAGACCACACTTCTCATTCATGTGGCACCAGCACAGAAAAACGTGCAAGATCTA ->peaae1cellline_NODE_200266_length_52_cov_1.393258_g183448_i0 -TCTAATTCTGTGCTTCGAGGACGATTTCTTGTCTAAGATGCACAGAATTAGG ->peaae1cellline_NODE_200272_length_44_cov_5.555556_g183454_i0 -GGCCATGAACAACCAGTCCTGTCCTGTCTTGTTTCTTCCTGGTA ->peaae1cellline_NODE_200273_length_42_cov_44.962025_g183455_i0 -TACCGAGCTGCGAGGCGGCCTGTATAAAAACGCTGTCAGTGT ->petxwholemale_NODE_12966_length_61_cov_3.076923_g12631_i0 -CCTACACTTTCAAACCTGAAAGGTTTCTGCCTGAGAGCAAAGACCTCCTGAAGCCGTTCT -C ->petxwholefemale_NODE_12365_length_47_cov_1.677778_g12066_i0 -ATATTTCGAATCAGATAACGCGATGACGGCCGTCAGCTTCCATCAGC ->petx24wholefemale_NODE_12274_length_66_cov_2.137615_g11799_i0 -CTACCGAGAGTCATTTTCTTATGCCGTTCTAATACTTATAGTAGGCTATAGAAAAAATAT -CTGTTC ->petx0wholemale_NODE_1133_length_67_cov_2.293103_g1132_i0 -GCAAGGTGTCTTACGGGTGGTGATGCCCACCAAGACCCTGCGCATCACCACCCGTAAGAC -ACCTTGT ->petx0wholefemale_NODE_618915_length_74_cov_1.682927_g519641_i0 -AGGCTGAACCAAGACACCATGACGACTTGCCATCTGCAACATGGCTGAATGTAAAATGTA -AATTGGTTCAGCCC ->petx0wholefemale_NODE_618917_length_73_cov_2.401639_g519643_i0 -TAGTATAAAAACATGTCAAAGTGGCAAGGTAAACTGTGTAAACATGGCAGTAGCACTCAA -GGTTGGCATTTAT ->petx0wholefemale_NODE_618918_length_72_cov_3.909091_g519644_i0 -TTCGCCTAACAGCATTAGCAGCAGCCAAGAAGCAAATTGCCAGCAACAAGGTTGGTAAAA -GCAGGGGGCGAC ->petx0wholefemale_NODE_618919_length_72_cov_3.148760_g519645_i0 -TAATTGCGCGACGTTTGTCTCTCGACGGGATTGCACCAACTTTCACTGAAGGGCTTTCAC -GTGGCGTACTTC ->petx0wholefemale_NODE_618920_length_72_cov_3.090909_g519646_i0 -TTTAGCTTGTTTGGTTAGTTGTGTCTTTTGGTCGTTTGGCTCGCAAAGTCAGAACTAACC -AAACAAGCGAAC ->petx0wholefemale_NODE_618921_length_72_cov_2.272727_g519647_i0 -CCCAACCAGGGGAAATCAGCAGTTGCTTTTTCCTGTCTCCCCCTCCATCTTTCACTTTCA -TATCTCTCCTTC ->petx0wholefemale_NODE_618923_length_71_cov_12.450000_g519649_i0 -TATTCTGCTGACATAATAAAAAACTGTCAAATATGGGCCCTCGCACAGTTTTTTATTATG -ACAGCAGAATC ->petx0wholefemale_NODE_618925_length_71_cov_3.791667_g519651_i0 -TCGTGTTTTCAGTTCATAACCCACCTGTGGCAAGCATGGTGTTTATGATGGGTTATGAAC -TGAATACACGT ->petx0wholefemale_NODE_618926_length_70_cov_30.613445_g519652_i0 -TTAACATATTATTGAAACGGGACACTGATAACGCAGCCTTGATCTGTGTCCCGTTTCAAT -AATATGTTAT ->petx0wholefemale_NODE_618927_length_70_cov_1.722689_g519653_i0 -AACAGTTTTGACAAAGCTGCTCGCAAAACCAGCTGATGGCGATATTACCTGCATTTTATT -TTTTTTTACC ->petx0wholefemale_NODE_618928_length_70_cov_1.512605_g519654_i0 -ACGCACCTTTGGTCAATTTCATCCTATCGCTGTCACACGACGACACGACAGTATTTCTAG -AATTCTGCAG ->petx0wholefemale_NODE_618929_length_69_cov_4.567797_g519655_i0 -TAGAATGTTATTCAGGTCCTCCACTACATCCCCCAGCCGGATGTAGTGGAGGACCTGAAT -AACATTCTT ->petx0wholefemale_NODE_618931_length_68_cov_8.094017_g519657_i0 -GAACGAGTGCAATTCGCTTGTCCCATTCGATGGGGCTTGGAAATGGGCCAAGCGAATTGC -ACTCGTTA ->petx0wholefemale_NODE_618932_length_68_cov_3.153846_g519658_i0 -GTGGCAGTTGCCGCCATGCTCGGCTTTTTTTTTATCTTTGTTTAACTGTGTGGTAAATTA -TTAAAATT ->petx0wholefemale_NODE_618933_length_68_cov_2.717949_g519659_i0 -GCGGTCTGGCCGGTCCGTATTTATTTTTCATCTTTGGAGCCGTTATCGTTTTTTTTTTTT -CACCCGGT ->petx0wholefemale_NODE_618934_length_68_cov_2.282051_g519660_i0 -ATAAAACCGGGGCGAGGCATGACCCAATGCATATCTCCTTCTTTTGCAGATAAAACTCCT -TGTCAAAG ->petx0wholefemale_NODE_618935_length_68_cov_1.145299_g519661_i0 -CAAGTTACCACAATACAATGCTCAAAAAATAAAGGGTTGCATTTTGAGCATTGTATTGGG -GTAACGTT ->petx0wholefemale_NODE_618937_length_67_cov_10.827586_g519663_i0 -AGACATCTGCTCTATGTGTGCTGTCATACTGACTATTCATATGGCAGCACACATAGAGCC -GAGGTCG ->petx0wholefemale_NODE_618939_length_67_cov_1.206897_g519665_i0 -CCTATACTGCTTAAGAGCAGTCAGGGCTGGATATCTCACTGTTGGGCCTGCTCTTAAGCA -GTATAGA ->petx0wholefemale_NODE_618940_length_67_cov_1.172414_g519666_i0 -CTTGCCTTTCTCGGGTAGAGGTATCGTCCCTGGTTATAAAGCCAGCAACAGCAAAAGAAA -TGCAGTA ->petx0wholefemale_NODE_618941_length_66_cov_1.008696_g519667_i0 -CCTTTAATTTTAAGACTCACTTCCGCTCCCTCGAGCACAGCAGCAGGTGAGTCTTAAAAT -TAAAGT ->petx0wholefemale_NODE_618944_length_65_cov_5.570175_g519670_i0 -GCCGTCCATCTCTATACACGACACAGACATGTGTGTGGCCGCGTTTCTGTCGCTGAGTTG -CTGAG ->petx0wholefemale_NODE_618945_length_65_cov_4.210526_g519671_i0 -GAAGGACGCATTAACGGAGCTCTCTTGGATTCACTGTCCGCTGGTGGCATTGATTACAGC -GTCCA ->petx0wholefemale_NODE_618947_length_65_cov_2.035088_g519673_i0 -GCTTCATGCAGAATGAAGAGAAACCGGAACGGATGGCAGCAGGTAGCATTTTGTCTGCAT -GAAGG ->petx0wholefemale_NODE_618948_length_65_cov_2.035088_g519674_i0 -TTCTGTCTCGCTTACGCCATTCTCGCCACGTTGGCAGCTGTCAAGACTCTCGCTTCGCCA -GTTTT ->petx0wholefemale_NODE_618949_length_63_cov_1.366071_g519675_i0 -CAAAAGCGTACAACTGTCAGCGACGAGCCACGAGCGCCCTGCGCTGACAGTTGTACGCTT -TTT ->petx0wholefemale_NODE_618951_length_61_cov_38.827273_g519677_i0 -TTTATTCCACATATGTAAAATCAAGATGAACAAAAGTGTGGTGCTTTTCAGCTGGAAATA -T ->petx0wholefemale_NODE_618952_length_61_cov_1.945455_g519678_i0 -GCCATCTTAACGAGCTATTATTATTTTAATATTATAAATAATAATCGCGCGATAAGATGG -G ->petx0wholefemale_NODE_618953_length_61_cov_1.845455_g519679_i0 -GCTGAAAGGTCAGCAGCCGAGCGCTATAACCATTGCGCCACCACGGCGCCAACCTTTCAG -A ->petx0wholefemale_NODE_618954_length_61_cov_1.200000_g519680_i0 -AAAGATCCCCATCAAAGGTGGTTTTTGTCTCTCTTAAAAAGCCGCTTTGTTGGGGATGTT -G ->petx0wholefemale_NODE_618956_length_60_cov_1.394495_g519682_i0 -CTTCAAACAATAACTGCAGTTCGCGATTCACTGGTAACTCCCATGTAAACGTGAGGCCCT ->petx0wholefemale_NODE_618957_length_60_cov_1.229358_g519683_i0 -AGCTTTATTAGCCGTGCATAAAACGGAAGTCTTTCCGCATCCAGCAGGGCCAGTGAAGAA ->petx12wholefemaleecoli_NODE_281510_length_74_cov_17.659864_g256556_i0 -CTTTCACGTCCCTGTACATAGCCCTTCCTTTGTTTATGGCTTGTTAAATAAAAGCTGGTT -TACCTTGTTGAAAA From 74fc391128584d264843275acb1b324ac5f71289 Mon Sep 17 00:00:00 2001 From: Taylor Reiter Date: Wed, 5 Jun 2024 10:48:05 -0400 Subject: [PATCH 2/5] lint and format --- Snakefile | 1 + 1 file changed, 1 insertion(+) diff --git a/Snakefile b/Snakefile index 9470a88..1bf1288 100644 --- a/Snakefile +++ b/Snakefile @@ -47,6 +47,7 @@ Note that we follow the conventions of the prokka tool for output file suffixes ## sORF prediction ################################################################################ + rule get_coding_contig_names: """ Extract amino acid contig names and remove everything after the first period, From fe21372ecb440c566e5f7ff9536e729737e25a8c Mon Sep 17 00:00:00 2001 From: Taylor Reiter Date: Wed, 5 Jun 2024 11:25:42 -0400 Subject: [PATCH 3/5] rev order of demo data cat --- demo/contigs.fa | 270 ++++++++++++++++++++++++------------------------ 1 file changed, 135 insertions(+), 135 deletions(-) diff --git a/demo/contigs.fa b/demo/contigs.fa index 6496905..994b6ac 100644 --- a/demo/contigs.fa +++ b/demo/contigs.fa @@ -1,3 +1,138 @@ +>peaae2cellline_NODE_191792_length_57_cov_5.882979_g176311_i0 +AAAATACGAAATGGTGCTTTTTGGTTTGCAATAAACACAGTTAGGTGTGCCTGTTTC +>peaae2cellline_NODE_191793_length_57_cov_2.702128_g176312_i0 +GGATTCGAAGTCGGGCAAAGAGCAGGGAACATTCCACGCTCTCAGTTATGGTGTCAT +>peaae2cellline_NODE_191799_length_47_cov_1.500000_g176318_i0 +AATAGCTGTCGAATTCTGCCACTTTCATTATGCCATCTTCGACGGGA +>peaae2cellline_NODE_191803_length_46_cov_1.710843_g176322_i0 +TACAGCACCACAGCACTTTTGCACCGTGGGTTTGCGCCTTGCCTGC +>peaae2cellline_NODE_191804_length_46_cov_1.228916_g176323_i0 +CTTTGAGAAGTAATTAAATTACATTACAAATTACTCTTCTCAAAAA +>peaae2cellline_NODE_191809_length_38_cov_1.253333_g176328_i0 +CCTTGTCTTGAAAATGAAAACCGGTTTTGGGGAAAAGA +>peaae1cellline_NODE_200257_length_72_cov_1.972477_g183439_i0 +TGTTCCTGTGGTGCTGGGTCTTTGCAACTGAAATAAATAAACAAGCCTAGACTGACCCAG +CACCACAGGAAC +>peaae1cellline_NODE_200261_length_59_cov_1.343750_g183443_i0 +AGGCAGTTACTTACTGCTTTGGCAGGACGGTCGAAGGGGCAAGGTAGACAATGTAGTGG +>peaae1cellline_NODE_200262_length_58_cov_2.178947_g183444_i0 +CTGTACAACGCACTGTTTTCCTCTGTTGCTAACTATTGTATTTTGGTCTGGGGTACAA +>peaae1cellline_NODE_200264_length_56_cov_1.892473_g183446_i0 +GAGAGACCACACTTCTCATTCATGTGGCACCAGCACAGAAAAACGTGCAAGATCTA +>peaae1cellline_NODE_200266_length_52_cov_1.393258_g183448_i0 +TCTAATTCTGTGCTTCGAGGACGATTTCTTGTCTAAGATGCACAGAATTAGG +>peaae1cellline_NODE_200272_length_44_cov_5.555556_g183454_i0 +GGCCATGAACAACCAGTCCTGTCCTGTCTTGTTTCTTCCTGGTA +>peaae1cellline_NODE_200273_length_42_cov_44.962025_g183455_i0 +TACCGAGCTGCGAGGCGGCCTGTATAAAAACGCTGTCAGTGT +>petxwholemale_NODE_12966_length_61_cov_3.076923_g12631_i0 +CCTACACTTTCAAACCTGAAAGGTTTCTGCCTGAGAGCAAAGACCTCCTGAAGCCGTTCT +C +>petxwholefemale_NODE_12365_length_47_cov_1.677778_g12066_i0 +ATATTTCGAATCAGATAACGCGATGACGGCCGTCAGCTTCCATCAGC +>petx24wholefemale_NODE_12274_length_66_cov_2.137615_g11799_i0 +CTACCGAGAGTCATTTTCTTATGCCGTTCTAATACTTATAGTAGGCTATAGAAAAAATAT +CTGTTC +>petx0wholemale_NODE_1133_length_67_cov_2.293103_g1132_i0 +GCAAGGTGTCTTACGGGTGGTGATGCCCACCAAGACCCTGCGCATCACCACCCGTAAGAC +ACCTTGT +>petx0wholefemale_NODE_618915_length_74_cov_1.682927_g519641_i0 +AGGCTGAACCAAGACACCATGACGACTTGCCATCTGCAACATGGCTGAATGTAAAATGTA +AATTGGTTCAGCCC +>petx0wholefemale_NODE_618917_length_73_cov_2.401639_g519643_i0 +TAGTATAAAAACATGTCAAAGTGGCAAGGTAAACTGTGTAAACATGGCAGTAGCACTCAA +GGTTGGCATTTAT +>petx0wholefemale_NODE_618918_length_72_cov_3.909091_g519644_i0 +TTCGCCTAACAGCATTAGCAGCAGCCAAGAAGCAAATTGCCAGCAACAAGGTTGGTAAAA +GCAGGGGGCGAC +>petx0wholefemale_NODE_618919_length_72_cov_3.148760_g519645_i0 +TAATTGCGCGACGTTTGTCTCTCGACGGGATTGCACCAACTTTCACTGAAGGGCTTTCAC +GTGGCGTACTTC +>petx0wholefemale_NODE_618920_length_72_cov_3.090909_g519646_i0 +TTTAGCTTGTTTGGTTAGTTGTGTCTTTTGGTCGTTTGGCTCGCAAAGTCAGAACTAACC +AAACAAGCGAAC +>petx0wholefemale_NODE_618921_length_72_cov_2.272727_g519647_i0 +CCCAACCAGGGGAAATCAGCAGTTGCTTTTTCCTGTCTCCCCCTCCATCTTTCACTTTCA +TATCTCTCCTTC +>petx0wholefemale_NODE_618923_length_71_cov_12.450000_g519649_i0 +TATTCTGCTGACATAATAAAAAACTGTCAAATATGGGCCCTCGCACAGTTTTTTATTATG +ACAGCAGAATC +>petx0wholefemale_NODE_618925_length_71_cov_3.791667_g519651_i0 +TCGTGTTTTCAGTTCATAACCCACCTGTGGCAAGCATGGTGTTTATGATGGGTTATGAAC +TGAATACACGT +>petx0wholefemale_NODE_618926_length_70_cov_30.613445_g519652_i0 +TTAACATATTATTGAAACGGGACACTGATAACGCAGCCTTGATCTGTGTCCCGTTTCAAT +AATATGTTAT +>petx0wholefemale_NODE_618927_length_70_cov_1.722689_g519653_i0 +AACAGTTTTGACAAAGCTGCTCGCAAAACCAGCTGATGGCGATATTACCTGCATTTTATT +TTTTTTTACC +>petx0wholefemale_NODE_618928_length_70_cov_1.512605_g519654_i0 +ACGCACCTTTGGTCAATTTCATCCTATCGCTGTCACACGACGACACGACAGTATTTCTAG +AATTCTGCAG +>petx0wholefemale_NODE_618929_length_69_cov_4.567797_g519655_i0 +TAGAATGTTATTCAGGTCCTCCACTACATCCCCCAGCCGGATGTAGTGGAGGACCTGAAT +AACATTCTT +>petx0wholefemale_NODE_618931_length_68_cov_8.094017_g519657_i0 +GAACGAGTGCAATTCGCTTGTCCCATTCGATGGGGCTTGGAAATGGGCCAAGCGAATTGC +ACTCGTTA +>petx0wholefemale_NODE_618932_length_68_cov_3.153846_g519658_i0 +GTGGCAGTTGCCGCCATGCTCGGCTTTTTTTTTATCTTTGTTTAACTGTGTGGTAAATTA +TTAAAATT +>petx0wholefemale_NODE_618933_length_68_cov_2.717949_g519659_i0 +GCGGTCTGGCCGGTCCGTATTTATTTTTCATCTTTGGAGCCGTTATCGTTTTTTTTTTTT +CACCCGGT +>petx0wholefemale_NODE_618934_length_68_cov_2.282051_g519660_i0 +ATAAAACCGGGGCGAGGCATGACCCAATGCATATCTCCTTCTTTTGCAGATAAAACTCCT +TGTCAAAG +>petx0wholefemale_NODE_618935_length_68_cov_1.145299_g519661_i0 +CAAGTTACCACAATACAATGCTCAAAAAATAAAGGGTTGCATTTTGAGCATTGTATTGGG +GTAACGTT +>petx0wholefemale_NODE_618937_length_67_cov_10.827586_g519663_i0 +AGACATCTGCTCTATGTGTGCTGTCATACTGACTATTCATATGGCAGCACACATAGAGCC +GAGGTCG +>petx0wholefemale_NODE_618939_length_67_cov_1.206897_g519665_i0 +CCTATACTGCTTAAGAGCAGTCAGGGCTGGATATCTCACTGTTGGGCCTGCTCTTAAGCA +GTATAGA +>petx0wholefemale_NODE_618940_length_67_cov_1.172414_g519666_i0 +CTTGCCTTTCTCGGGTAGAGGTATCGTCCCTGGTTATAAAGCCAGCAACAGCAAAAGAAA +TGCAGTA +>petx0wholefemale_NODE_618941_length_66_cov_1.008696_g519667_i0 +CCTTTAATTTTAAGACTCACTTCCGCTCCCTCGAGCACAGCAGCAGGTGAGTCTTAAAAT +TAAAGT +>petx0wholefemale_NODE_618944_length_65_cov_5.570175_g519670_i0 +GCCGTCCATCTCTATACACGACACAGACATGTGTGTGGCCGCGTTTCTGTCGCTGAGTTG +CTGAG +>petx0wholefemale_NODE_618945_length_65_cov_4.210526_g519671_i0 +GAAGGACGCATTAACGGAGCTCTCTTGGATTCACTGTCCGCTGGTGGCATTGATTACAGC +GTCCA +>petx0wholefemale_NODE_618947_length_65_cov_2.035088_g519673_i0 +GCTTCATGCAGAATGAAGAGAAACCGGAACGGATGGCAGCAGGTAGCATTTTGTCTGCAT +GAAGG +>petx0wholefemale_NODE_618948_length_65_cov_2.035088_g519674_i0 +TTCTGTCTCGCTTACGCCATTCTCGCCACGTTGGCAGCTGTCAAGACTCTCGCTTCGCCA +GTTTT +>petx0wholefemale_NODE_618949_length_63_cov_1.366071_g519675_i0 +CAAAAGCGTACAACTGTCAGCGACGAGCCACGAGCGCCCTGCGCTGACAGTTGTACGCTT +TTT +>petx0wholefemale_NODE_618951_length_61_cov_38.827273_g519677_i0 +TTTATTCCACATATGTAAAATCAAGATGAACAAAAGTGTGGTGCTTTTCAGCTGGAAATA +T +>petx0wholefemale_NODE_618952_length_61_cov_1.945455_g519678_i0 +GCCATCTTAACGAGCTATTATTATTTTAATATTATAAATAATAATCGCGCGATAAGATGG +G +>petx0wholefemale_NODE_618953_length_61_cov_1.845455_g519679_i0 +GCTGAAAGGTCAGCAGCCGAGCGCTATAACCATTGCGCCACCACGGCGCCAACCTTTCAG +A +>petx0wholefemale_NODE_618954_length_61_cov_1.200000_g519680_i0 +AAAGATCCCCATCAAAGGTGGTTTTTGTCTCTCTTAAAAAGCCGCTTTGTTGGGGATGTT +G +>petx0wholefemale_NODE_618956_length_60_cov_1.394495_g519682_i0 +CTTCAAACAATAACTGCAGTTCGCGATTCACTGGTAACTCCCATGTAAACGTGAGGCCCT +>petx0wholefemale_NODE_618957_length_60_cov_1.229358_g519683_i0 +AGCTTTATTAGCCGTGCATAAAACGGAAGTCTTTCCGCATCCAGCAGGGCCAGTGAAGAA +>petx12wholefemaleecoli_NODE_281510_length_74_cov_17.659864_g256556_i0 +CTTTCACGTCCCTGTACATAGCCCTTCCTTTGTTTATGGCTTGTTAAATAAAAGCTGGTT +TACCTTGTTGAAAA >Transcript_0 len=4692 CDS=4152-4692 exon=0-4692 five_prime_UTR=0-4152 gene=0-4692 mRNA=0-4692 TCGATCCCGGCCGCGGTGGTTGAATTTCGATGGAGGCGAAATTCTAGAGGCCTGTGTGCTGTACGATATCAGTGCACGGTAAAGAACCCCAGATGGTTGAAATTCTTGAAGCCTCCCACTACAACGTCCCTCATAGCTTGAGTCACTTTGGGACATTAAACCATCATACACCATAACTGCCATCTCCCACGGTGGTCAGTTTGCTCACAATATAGTGGGCAAGCTGTGTTGACGCCACGAGACATGAACTAGGTGGCGCGATGCCCTTCTCCGCACAAGACAAGATTCTTGCTAGCGGTTGGTACATATGGAGTGGTGCTTCTGCACAAAAACTGGAGCACACACTTAAGCTCCGCCTAAAGGGTATGACCCGACAGTGTTGATGGGCTAATGCCCACGTATGGGAAAGTGGTCATGCTCTACTTTATATTTATACATCCCTGGGAGTCCTCGTACCTCTCCTGGTGCAGTGGTTAAACAATGTGACACTGCTCTGCGATGGCAGCTGTCGTCACTGGTGGGGCGTGTGTGACCTTGGTTTGTATTCCCAAGCAATCTCTCCTAACCAATCATTCATTAAATTGCCACCTGTAGCGGGGCAGTTTGCTCGCGATCCGGTCGGCAGGTTATGATGGTGTTACTCACGTGAGTTAGGTTGCTTCTCTGCGGTTTTTTCATCACTGTTTTCACTTGCAACCAGTGACGCCAGCAACGACGCTGGCTTTTCTGTGATATGGTGCTTTTAACACCGTCGCTTTCAAGGCGGATTCACATTGGCGAATCGCGGAGGTCGCGTGACCGCATTCAATCAAATTGCTATGCTCCTTCCCTGCTTTCCTAGTGCTTCGGAAACCGTATCTAATTACTTGAAGTCGGTCGCGTGATGTTTGTGACTCACCAGTGTGAAACCACTTTAAAACTGACCTGTTGCAGCTGCTGTCTTGTGAAGTCTGCAGTTCCACATTTTCATGCATTCAGGGATGCAAAGCATCCATATTTCATGATATAAGATGCAGTGATGTTTGATTTGATTTTTTCTTGCTTGTTTGTAATATTTATATGAGCTTCTGTATGTCGTTTTTAGTTTAAAATGTTTGTGGACAACTTGAGTCGGTATTCTACACTTTATGAGTAGGACTCAGAGCAGTAAATAGCTATATTTTTGTTATTGTTCCTACAAGAAGTTTGGCCTTTTTTCTTTTTCTCTTGCAACTAGTCAGCATTTCTAGCGTGGCTCAGCCATGCTTTCTCAACATTTCTTTTAACATTTTCAAAGAGACTAACTGCTGCCACTTGCCACCAGTGGTGTTGTGCTGCTGTTCGCAGAAGCGGTTGCTGGTATCTATGATAAATATTGTCATGTAGTGGTGATGGCAGTGTGGTGAACAAGAAGACAGCAAAAGGCTTCCAGAATGACTGCAACTCTTTGGCGAAATCAATTTAAACCTTTGTAAGGAACCATTGAGATAAGGGACGAAAAGAAACTATGACAATCTGGAAAAAGGTGGAAGCAGTTCCGCGTGGATGAATGGAGATATGTCTGGTCGCATCTCGCATCACAGACCAAATGATGAAGCCATATAGCAGTGGCTTTGAGCATCAACAAACAGTACAAAGATTGGCAGCAATATCAGCAGTAATAATTAAGTTATATAGTATTTTTAAGACCGCCTTTGTGTTTTCACCTGTAGTTTTGATAAAACTTGTGCATCTCATACAATGCCAAGTTCATTACTGCCTCTGTTTGTTGCCTGGGTGGCCCGGTAGTCCTGTCCATAGAACTTCATAGAGTTGCTATACTTGTAAAGTGTGTCCTGAGTGTAAGTGTTTGAAGCCACTGTCTGCCATGAGAAAACAAGACTACGATATCAAAATATTCGATTTAAAATGCTTCTGCTACGTATCTCAGTGCCACAGCGTAGATTAGCTTATTTGGACAGAGCTTTCGATTGCTGCAAAAAAAAAAAATGAATAAAAGTCAGCAGACTGTTCCTTAACCAAACTACTACAGTTTGCAGCCTTGTAGTTTTTAATCTTATCTACTTCTCTGTTTCATTGTGCAGGTATCGCAACAAAGATTGGCAACCAATGAGTACACCAGCTGCCCGAAGGACCATTTTCCAAACCCACTGCTTCCAGTATTTGCAGGCTAAACTGCTTGATGTTCATAACGTGAAATGGCATAGCCAGTGATATCAGAATGGTGACGGTGCTCAAAATGAGAGATATGAGAGCCTGCCTTGTTAATGCATAGACCACAACCATCATTTCAACAAGCCACAGTCTCTGTCATTGCAGTGCTCAGCTGCTGACCTGAAAGATGTGGGTTCAGTCCTGGCCGAGGCAAAATGTTAGGAGGCCCATGTACTTTGCGATGCCAGTGCACGTTAAAGAACTGCAGGTGGTCCTAATTACCGGAACCCCTCCCCTATGGCATCCCTCATAGCCTGAGTCGCTGTGGGATGTTAAACCCAACAAAACAATCAATCATCACTGCGACATGAACAGTGGATGTTGTGAAACATTCTGTATGTGAAATTGTTCATGAATGCAGAGGAAACCTTGTGATCTAGCACATCGTGCTACAAGTTGCACTTTTTCTCTATATTCTGTGCAGGTTTTTCTGAAAAAAGTACACAAGCTGGAGAGTACACAGAAAAAAAAGAACACCTGCAGTAAAAAAGTGCCTGTATTTTGCGGCTGTGCCTATTGTGCATCATAATTGGTCTTGCTTTTCATGCTCTAGTGTCCCCAATTTAGGTCTTCCATTTCTATTGTGGTTCTACCACCCTTTGAGTTTCCGTGTCTCATTATTTTATGGCTTATTGTTTTCTGTATTCTTCTGTATACGACACCCTTATAAAGTGAGGTCTGAAAAGCAGTAGTAATATCTTTTTGGTATGCTGTCGGATGTTCTGCAGTGAACTATTTTATCTGATGTCTACATTAGCTTGTATTGCATGTGCTTCTGATTCTGTTTGCTGTTAGCAAATATTTTCGAGTCCGGTTAATTGTGAGATGATATTGATGATTCGTGCTTATTTTAGTGCTGTGTATTACAGTGGTGACACTGCCGATAGAATTGGTTCTGTCCTTGTCGTTTGCGGCATGGCTGCTTGTGTTCTGAATTTATTTTTTGAGATCATCAGGCCTGATTATGCTAGATTGAAAAGAAAAGTTTTGTCTTCGTAATTGCTTTTTATTTTGCTCTCTTAACATTCAGAAATTCTGTGCATGTTTTAGCTATTCATTCATTGCAATTCAGTTTGGTAAAAGTATCAATCTGCACCTATTTTAAGTCATGAAGCCTTGCTTTATTTTTTTTGAGGCTGCAGGAAGAAAATTGCTTCCTCATGATTCTGGCAGAGCAGATGGAGGTTCGTGTTCCGGCACAATAAACTGAACACGGTAGCTGCCAGACGTGCACACTGTGAAATGGTGCCCGAAGCAAATGATAATTACTGCGTTGTTTTAGGTCATGACAGGTCAGCTCAGATATATATATTTTAGCATTTGGCTACATTTAGCTTAAGAGAACTAGTAATCATTGAGCAGTTTTTCTCCAGCAATGTGCTAGAGGTTTTTTCTGCAACATGCAGAGTAATTAACCTTGAAGCTTTTATGGACAAAAGGAATGCACTCATGAGCAGTGTACTCTACCATGAGTCAAAGCTCGCAGACATTGTTATTCGTGAAAGACTTGCTCAGAATTGCAAGATCACTCCGATGAAGCCAAGCATTATATGAGGACTTTGCTGTCCTGGTGCGTAAGTTGGCTAAAGACTGCTCATCGCTCTCAGTAATTGTTCGAAAGCTGAGACATTATCCCTCCAACTCCATGCTAATTATAAATATCAGAAAATTTGGTTTAAAAACTTTCTTTAGAGGCACAACCCATTGACCATGGAATGAACTCTTAGCCTAAATGTCGTGTAATATTGCCATTTTCATTTGCAGCTTACTGGCCGTGCTACCAAATGGGCTCTAATACGCTGAAAATACTCTCCTTTTCTTGGTGCAGGCCTGGCCACGGCTGATCCCAAGCTGAGAGTGTCGCAAGGTGAAAACAGGACTAAGCAGCTATACCAGTCCTCCCTCTGCAGTCGGAGTTTTATTAAGAAAAGCCCTCTAGAGCCGCACCTTCTCACCCAAATGAGTGATCATCCTTTCCAGTGTGGCCACTGTGGAAAGAGCTTTCCACAAAAGTATGACTCAGTGCAACATCTCTGTACCCACACATATGAGCGTCCATACAAGTGTAGCCAGTGCAGGAGCAGCTTTACTCAAAAGGCCCACCTGGTGCAACACCTTCGTACCCACATGGATGAGCATCTGTACAGGTGTGACCACTGTGAGAGCAGCTTTGCTCTAAAGGTTCACCTGGAGCGACACCTTCTTACTCACACAGGTGGGTGTCCATGCAAGTGTGCCCACTGTGACAGCAGCTTTGTTCAAAAGAGCCATCTGAAACATCTTTGTACCCACACGGGTGAGCGTCCATACAGGTGCAACCACTGTGACAGCAGCTTTGCTCAAAAGGGCACCCTGGAGCAACATCTTCGTACCCACACGGGTGAGTGTCCATACAAGTGTGACCATTGTGACAGCAGCTTTTTTGAAAAAGGCAACCTGAAGCAACACCCTTGTACCCACACGGGTGAGCGTCCATACAAGTGTGACCACTGTGACAGC >Transcript_1 len=4685 CDS=3755-4031,4322-4685 exon=0-4685,0-4685 five_prime_UTR=4031-4685 gene=0-4685,0-4685 mRNA=0-4685,0-4685 three_prime_UTR=0-3755,0-4322 @@ -198,138 +333,3 @@ GGAAAGCTAGGCCAAGAGCCAAGCCACATCAAAATTCAAAGTGGCTTCTAGAACTCGTTGCCACAGTTGTGATGTAAAGA GTCGGGTTCGAACCCGGGAACTCCGGATCAGTAGTTGGACGCTCCAACCACTGAGCCACCGCGGCGGGGATCGCACGTGGATAATCAATTTCTCCGCCGCGAAGCGCTACCGAAGGTATGCAGACGCACACGTCACTCTTTGGTAGGATGTTGAAGGCTTCTTCTATCTCCCTGGTCCGACCAAAATATGACGGGATGACCACCGTGTTGACAAGCTGCGGAGCGCACTCATGCACCGTCCGATGCCCCGCTAATTGACCGCTGATAGCGCACTCCAGGGATCCCCCGTGCTCTCGAAGCCGTTCATTCAAACAGCGCTCCGACCAATGTAGGTGCGTCCGCATGACAAAGGTATTCCACAAACCCCGCTTTTCACACATTTCACGCCCTCTCCGGTGTTCCCTCCTAAACTGATTCTTTT >Transcript_99 len=421 CDS=1-421 exon=0-421 gene=0-421 mRNA=0-421 CAAGGACAAGGATGGCAAGAACTTGCTGAAGGTGGTTATGAGGACATGGCTTCCTGCTGGAGATGCACTCTTCGAAATGATCACAATCCACCTTCCATCTCCTGTGACGGCACAACGCTACCGAATGGAAATCCTCTATGAGGGTCCCCTTGATGATGAGGCTGCTGTTGCTGTAAAGGCGTGCGACCCAGAAGGCCCACTGATGATGTACGTCTCAAAGATGGTACCGACATCGGACAAGGGCCGCTTTTACGCATTTGGACGAGTCTTCTCTGGTGTTGTGTCATCGGGACAAAAGGTCCGCATCATGGGCCCAAACTACACTCCTGGCAAGAAAGAGGACTTGGCTGAGAAGGCCATTCAGAGGACTGTTCTGATGATGGGTCGGTACGTGGAGCCCATTGAAGATGTGCCCTGCGGT ->peaae2cellline_NODE_191792_length_57_cov_5.882979_g176311_i0 -AAAATACGAAATGGTGCTTTTTGGTTTGCAATAAACACAGTTAGGTGTGCCTGTTTC ->peaae2cellline_NODE_191793_length_57_cov_2.702128_g176312_i0 -GGATTCGAAGTCGGGCAAAGAGCAGGGAACATTCCACGCTCTCAGTTATGGTGTCAT ->peaae2cellline_NODE_191799_length_47_cov_1.500000_g176318_i0 -AATAGCTGTCGAATTCTGCCACTTTCATTATGCCATCTTCGACGGGA ->peaae2cellline_NODE_191803_length_46_cov_1.710843_g176322_i0 -TACAGCACCACAGCACTTTTGCACCGTGGGTTTGCGCCTTGCCTGC ->peaae2cellline_NODE_191804_length_46_cov_1.228916_g176323_i0 -CTTTGAGAAGTAATTAAATTACATTACAAATTACTCTTCTCAAAAA ->peaae2cellline_NODE_191809_length_38_cov_1.253333_g176328_i0 -CCTTGTCTTGAAAATGAAAACCGGTTTTGGGGAAAAGA ->peaae1cellline_NODE_200257_length_72_cov_1.972477_g183439_i0 -TGTTCCTGTGGTGCTGGGTCTTTGCAACTGAAATAAATAAACAAGCCTAGACTGACCCAG -CACCACAGGAAC ->peaae1cellline_NODE_200261_length_59_cov_1.343750_g183443_i0 -AGGCAGTTACTTACTGCTTTGGCAGGACGGTCGAAGGGGCAAGGTAGACAATGTAGTGG ->peaae1cellline_NODE_200262_length_58_cov_2.178947_g183444_i0 -CTGTACAACGCACTGTTTTCCTCTGTTGCTAACTATTGTATTTTGGTCTGGGGTACAA ->peaae1cellline_NODE_200264_length_56_cov_1.892473_g183446_i0 -GAGAGACCACACTTCTCATTCATGTGGCACCAGCACAGAAAAACGTGCAAGATCTA ->peaae1cellline_NODE_200266_length_52_cov_1.393258_g183448_i0 -TCTAATTCTGTGCTTCGAGGACGATTTCTTGTCTAAGATGCACAGAATTAGG ->peaae1cellline_NODE_200272_length_44_cov_5.555556_g183454_i0 -GGCCATGAACAACCAGTCCTGTCCTGTCTTGTTTCTTCCTGGTA ->peaae1cellline_NODE_200273_length_42_cov_44.962025_g183455_i0 -TACCGAGCTGCGAGGCGGCCTGTATAAAAACGCTGTCAGTGT ->petxwholemale_NODE_12966_length_61_cov_3.076923_g12631_i0 -CCTACACTTTCAAACCTGAAAGGTTTCTGCCTGAGAGCAAAGACCTCCTGAAGCCGTTCT -C ->petxwholefemale_NODE_12365_length_47_cov_1.677778_g12066_i0 -ATATTTCGAATCAGATAACGCGATGACGGCCGTCAGCTTCCATCAGC ->petx24wholefemale_NODE_12274_length_66_cov_2.137615_g11799_i0 -CTACCGAGAGTCATTTTCTTATGCCGTTCTAATACTTATAGTAGGCTATAGAAAAAATAT -CTGTTC ->petx0wholemale_NODE_1133_length_67_cov_2.293103_g1132_i0 -GCAAGGTGTCTTACGGGTGGTGATGCCCACCAAGACCCTGCGCATCACCACCCGTAAGAC -ACCTTGT ->petx0wholefemale_NODE_618915_length_74_cov_1.682927_g519641_i0 -AGGCTGAACCAAGACACCATGACGACTTGCCATCTGCAACATGGCTGAATGTAAAATGTA -AATTGGTTCAGCCC ->petx0wholefemale_NODE_618917_length_73_cov_2.401639_g519643_i0 -TAGTATAAAAACATGTCAAAGTGGCAAGGTAAACTGTGTAAACATGGCAGTAGCACTCAA -GGTTGGCATTTAT ->petx0wholefemale_NODE_618918_length_72_cov_3.909091_g519644_i0 -TTCGCCTAACAGCATTAGCAGCAGCCAAGAAGCAAATTGCCAGCAACAAGGTTGGTAAAA -GCAGGGGGCGAC ->petx0wholefemale_NODE_618919_length_72_cov_3.148760_g519645_i0 -TAATTGCGCGACGTTTGTCTCTCGACGGGATTGCACCAACTTTCACTGAAGGGCTTTCAC -GTGGCGTACTTC ->petx0wholefemale_NODE_618920_length_72_cov_3.090909_g519646_i0 -TTTAGCTTGTTTGGTTAGTTGTGTCTTTTGGTCGTTTGGCTCGCAAAGTCAGAACTAACC -AAACAAGCGAAC ->petx0wholefemale_NODE_618921_length_72_cov_2.272727_g519647_i0 -CCCAACCAGGGGAAATCAGCAGTTGCTTTTTCCTGTCTCCCCCTCCATCTTTCACTTTCA -TATCTCTCCTTC ->petx0wholefemale_NODE_618923_length_71_cov_12.450000_g519649_i0 -TATTCTGCTGACATAATAAAAAACTGTCAAATATGGGCCCTCGCACAGTTTTTTATTATG -ACAGCAGAATC ->petx0wholefemale_NODE_618925_length_71_cov_3.791667_g519651_i0 -TCGTGTTTTCAGTTCATAACCCACCTGTGGCAAGCATGGTGTTTATGATGGGTTATGAAC -TGAATACACGT ->petx0wholefemale_NODE_618926_length_70_cov_30.613445_g519652_i0 -TTAACATATTATTGAAACGGGACACTGATAACGCAGCCTTGATCTGTGTCCCGTTTCAAT -AATATGTTAT ->petx0wholefemale_NODE_618927_length_70_cov_1.722689_g519653_i0 -AACAGTTTTGACAAAGCTGCTCGCAAAACCAGCTGATGGCGATATTACCTGCATTTTATT -TTTTTTTACC ->petx0wholefemale_NODE_618928_length_70_cov_1.512605_g519654_i0 -ACGCACCTTTGGTCAATTTCATCCTATCGCTGTCACACGACGACACGACAGTATTTCTAG -AATTCTGCAG ->petx0wholefemale_NODE_618929_length_69_cov_4.567797_g519655_i0 -TAGAATGTTATTCAGGTCCTCCACTACATCCCCCAGCCGGATGTAGTGGAGGACCTGAAT -AACATTCTT ->petx0wholefemale_NODE_618931_length_68_cov_8.094017_g519657_i0 -GAACGAGTGCAATTCGCTTGTCCCATTCGATGGGGCTTGGAAATGGGCCAAGCGAATTGC -ACTCGTTA ->petx0wholefemale_NODE_618932_length_68_cov_3.153846_g519658_i0 -GTGGCAGTTGCCGCCATGCTCGGCTTTTTTTTTATCTTTGTTTAACTGTGTGGTAAATTA -TTAAAATT ->petx0wholefemale_NODE_618933_length_68_cov_2.717949_g519659_i0 -GCGGTCTGGCCGGTCCGTATTTATTTTTCATCTTTGGAGCCGTTATCGTTTTTTTTTTTT -CACCCGGT ->petx0wholefemale_NODE_618934_length_68_cov_2.282051_g519660_i0 -ATAAAACCGGGGCGAGGCATGACCCAATGCATATCTCCTTCTTTTGCAGATAAAACTCCT -TGTCAAAG ->petx0wholefemale_NODE_618935_length_68_cov_1.145299_g519661_i0 -CAAGTTACCACAATACAATGCTCAAAAAATAAAGGGTTGCATTTTGAGCATTGTATTGGG -GTAACGTT ->petx0wholefemale_NODE_618937_length_67_cov_10.827586_g519663_i0 -AGACATCTGCTCTATGTGTGCTGTCATACTGACTATTCATATGGCAGCACACATAGAGCC -GAGGTCG ->petx0wholefemale_NODE_618939_length_67_cov_1.206897_g519665_i0 -CCTATACTGCTTAAGAGCAGTCAGGGCTGGATATCTCACTGTTGGGCCTGCTCTTAAGCA -GTATAGA ->petx0wholefemale_NODE_618940_length_67_cov_1.172414_g519666_i0 -CTTGCCTTTCTCGGGTAGAGGTATCGTCCCTGGTTATAAAGCCAGCAACAGCAAAAGAAA -TGCAGTA ->petx0wholefemale_NODE_618941_length_66_cov_1.008696_g519667_i0 -CCTTTAATTTTAAGACTCACTTCCGCTCCCTCGAGCACAGCAGCAGGTGAGTCTTAAAAT -TAAAGT ->petx0wholefemale_NODE_618944_length_65_cov_5.570175_g519670_i0 -GCCGTCCATCTCTATACACGACACAGACATGTGTGTGGCCGCGTTTCTGTCGCTGAGTTG -CTGAG ->petx0wholefemale_NODE_618945_length_65_cov_4.210526_g519671_i0 -GAAGGACGCATTAACGGAGCTCTCTTGGATTCACTGTCCGCTGGTGGCATTGATTACAGC -GTCCA ->petx0wholefemale_NODE_618947_length_65_cov_2.035088_g519673_i0 -GCTTCATGCAGAATGAAGAGAAACCGGAACGGATGGCAGCAGGTAGCATTTTGTCTGCAT -GAAGG ->petx0wholefemale_NODE_618948_length_65_cov_2.035088_g519674_i0 -TTCTGTCTCGCTTACGCCATTCTCGCCACGTTGGCAGCTGTCAAGACTCTCGCTTCGCCA -GTTTT ->petx0wholefemale_NODE_618949_length_63_cov_1.366071_g519675_i0 -CAAAAGCGTACAACTGTCAGCGACGAGCCACGAGCGCCCTGCGCTGACAGTTGTACGCTT -TTT ->petx0wholefemale_NODE_618951_length_61_cov_38.827273_g519677_i0 -TTTATTCCACATATGTAAAATCAAGATGAACAAAAGTGTGGTGCTTTTCAGCTGGAAATA -T ->petx0wholefemale_NODE_618952_length_61_cov_1.945455_g519678_i0 -GCCATCTTAACGAGCTATTATTATTTTAATATTATAAATAATAATCGCGCGATAAGATGG -G ->petx0wholefemale_NODE_618953_length_61_cov_1.845455_g519679_i0 -GCTGAAAGGTCAGCAGCCGAGCGCTATAACCATTGCGCCACCACGGCGCCAACCTTTCAG -A ->petx0wholefemale_NODE_618954_length_61_cov_1.200000_g519680_i0 -AAAGATCCCCATCAAAGGTGGTTTTTGTCTCTCTTAAAAAGCCGCTTTGTTGGGGATGTT -G ->petx0wholefemale_NODE_618956_length_60_cov_1.394495_g519682_i0 -CTTCAAACAATAACTGCAGTTCGCGATTCACTGGTAACTCCCATGTAAACGTGAGGCCCT ->petx0wholefemale_NODE_618957_length_60_cov_1.229358_g519683_i0 -AGCTTTATTAGCCGTGCATAAAACGGAAGTCTTTCCGCATCCAGCAGGGCCAGTGAAGAA ->petx12wholefemaleecoli_NODE_281510_length_74_cov_17.659864_g256556_i0 -CTTTCACGTCCCTGTACATAGCCCTTCCTTTGTTTATGGCTTGTTAAATAAAAGCTGGTT -TACCTTGTTGAAAA From 6db82c31ffbb6a70f09e5779d0b27dc27dc79bc0 Mon Sep 17 00:00:00 2001 From: Taylor Reiter Date: Wed, 5 Jun 2024 18:27:15 -0400 Subject: [PATCH 4/5] Apply suggestions from code review Co-authored-by: Keith Cheveralls Signed-off-by: Taylor Reiter --- README.md | 2 +- config.yml | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/README.md b/README.md index f3d1714..7357f35 100644 --- a/README.md +++ b/README.md @@ -35,7 +35,7 @@ snakemake --software-deployment-method conda -j 1 --configfile demo/config.yml ## Input data The [peptigate pipeline](./Snakefile) requires three input files: -* A transcriptome assembly: transcriptome assembly FASTA file in nucleotide format containing transcripts or contigs. +* A transcriptome assembly as a FASTA file in nucleotide format containing transcripts or contigs. * Open reading frames predicted from the transcriptome in both amino acid and nucleotide format. The open reading frames in both files should have the same names before the first period in the FASTA header name. Tools like [TransDecoder](https://github.com/TransDecoder/TransDecoder) provide these files in the correct format. diff --git a/config.yml b/config.yml index 7f2124c..94d0eaa 100644 --- a/config.yml +++ b/config.yml @@ -20,7 +20,7 @@ output_dir: "path/to/outputdir" # ORFs should be predicted from the same transcriptome assembly as the "contigs" input file. # ORFs should have the same name (before the first period in the name) as the contigs in the # "contigs" input file. TransDecoder provides files in the proper format. -# Used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases, and to +# This file is used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases, and to # remove coding transcripts from the transcriptome assembly before sORF prediction. # - orfs_nucleotides: predicted ORFs as nucleotide sequences. Should contain the same ORFs as # "orfs_amino_acids" but in nucleotide format. TransDecoder also provides this file in the proper From ec0bf2e9bffb8c5eac88948761d835159b14bc66 Mon Sep 17 00:00:00 2001 From: Taylor Reiter Date: Wed, 5 Jun 2024 18:31:23 -0400 Subject: [PATCH 5/5] suggestions from code review --- demo/README.md | 2 +- demo/config.yml | 29 ----------------------------- 2 files changed, 1 insertion(+), 30 deletions(-) diff --git a/demo/README.md b/demo/README.md index 8580f67..5397eed 100644 --- a/demo/README.md +++ b/demo/README.md @@ -9,4 +9,4 @@ head -n 204 orthofuser_final_clean.fa.transdecoder.cds > orfs_nucleotides.fa ``` We also pulled short contigs (less than 75 bp) from an internal S3 bucket and added these contigs to the `contigs.fa` file (50 contigs). -These are contigs that were filtered from the *Amblyomma* transcriptome prior to transcriptome merging (mid assembly pipeline). +These are contigs that were filtered from the *Amblyomma* transcriptome prior to transcriptome merging during the [transcriptome assembly pipeline](https://github.com/Arcadia-Science/2023-amblyomma-americanum-txome-assembly/). diff --git a/demo/config.yml b/demo/config.yml index 1fa5d27..248fe4d 100644 --- a/demo/config.yml +++ b/demo/config.yml @@ -1,34 +1,5 @@ -################################################################################ -## Input file descriptions -################################################################################ - -# Config parameters in this file are used as defaults by the pipeline (Snakefile). -# To override the defaults, create a copy of this file and pass your new file to snakemake using the --configfile flag. -# Any parameters in the new config file will overwrite the defaults listed here. - -########### -# File IO # -########### - -# Specify the input and output directories input_dir: "inputs/" output_dir: "outputs/demo" - -# Input files and directories -# - contigs: transcriptome assembly contigs. -# - orfs_amino_acids: predicted ORFs translated into amino acids. -# ORFs should be predicted from the same transcriptome assembly as the "contigs" input file. -# ORFs should have the same name (before the first period in the name) as the contigs in the -# "contigs" input file. TransDecoder provides files in the proper format. -# Used for cleavage peptide prediction and annotation of nonribosomal peptide synthetases, and to -# remove coding transcripts from the transcriptome assembly before sORF prediction. -# - orfs_nucleotides: predicted ORFs as nucleotide sequences. Should contain the same ORFs as -# "orfs_amino_acids" but in nucleotide format. TransDecoder also provides this file in the proper -# format. If this file contains short ORFs (< 300 nucleotides), they will not be reported as sORFs -# as they are already annotated in the input. -# - plmutils_model_dir: path to the directory for the plmutils model that will predict whether sORFs -# are coding or non-coding. - contigs: "demo/contigs.fa" orfs_amino_acids: "demo/orfs_amino_acids.faa" orfs_nucleotides: "demo/orfs_nucleotides.fa"