WrapFasta.pl - Wrapping sequence is recommended by NCBI for fasta format. Some programs will return errors if this is not done properly. Here are a couple examples of this type of error: "Each line of the fasta entry must be the same length except the last." or "Each line of the qual file must be less than 65,536 characters."
If you are receiving a similar error from a program you are trying to run then WrapFasta.pl will create a reformatted fasta file from your improperly wrapped file.
DEPENDENCIES:
WrapFasta.pl requires the following BioPerl libraries: Bio::SeqIO and Bio::Seq.
USAGE: WrapFasta.pl -f [FASTA filename]
fasta2fastq.pl - Script converts FASTA and QUAL files to FASTQ format. It is a slight modification of a script posted to SeqAnwsers http://seqanswers.com/forums/showthread.php?t=2775
USAGE: fasta2fastq -b [base filename e.g. filename minus .fasta] -d [directory]
filter_by_length.pl -Script output fasta records less then the set minimum length.
USAGE: filter_by_length.pl [FASTA filename or filenames]
pre_post_cleaning_metrics.pl - This script summarizes average read lengths etc before and after cleaning for multiple single end or paired end files. Output log files from prinseq (v prinseq-lite-0.20.3) are taken as the input and cleaning metrics are written to the file pre_post_clean_reads.csv.
USAGE: perl pre_post_cleaning_metrics.pl [.log filename/s]
run_broken_pair_removal.sh - Runs the broken pair removal script from the velvet package on multiple files. This script can convert illumina headers so that they end in "/1" or "/2" if your reads are not illumina cassava (post v1.8) you should comment out lines 28 and 29. This script uses a script called stats.sh to monitor and report resource usage on a SGE cluster like Beocat at KSU. But this section (lines 36-42) also can be commented out.
USAGE: run_broken_pair_removal.sh [filenames_minus_1.fastq or _2.fastq]
run_prinseq.sh - paired end read cleaning with prinseq. This script can convert illumina headers so that they end in "/1" or "/2" if your reads are not illumina cassava (post v1.8) you should comment out lines 28 and 29. This script uses a script called stats.sh to monitor and report resource usage on a SGE cluster like Beocat at KSU. But this section (lines 36-42) also can be commented out.
USAGE: bash run_prinseq.sh [filenames_minus_R1_001_h.fastq or _R2_001_h.fastq]
sanger2fastq.pl - Script converts older paired sanger reads (with an ANC, FASTA, and QUAL file) to FASTQ format. The headers and files will be similar to Illumina mate pair output (i.e. RF orientation). You will need to customize line 14 to 64 to import your data because the sanger ANC files seem to have a variable order to the columns. The conversion from FASTA to FASTQ is a modification of a script posted to SeqAnwsers http://seqanswers.com/forums/showthread.php?t=2775 Uses the '.anc' file to split reads into forward and reverse files and makes the header pairs identical except for a trailing /1 or /2
USAGE: sanger2fasta_qual.pl
sanger2fastq_brokenpair_rm.pl - Script output unbroken pairs and singletons, sequences with more than two template_ids in the anc file are removed
USAGE: sanger2fastq_brokenpair_rm.pl
make_contigs_from_fasta.pl Script breaks scaffolds at Ns and outputs a new multi fasta file with one entry pre contig
USAGE: perl make_contigs_from_fasta.pl [fasta]
stats.sh - see prinseq.sh