Perl script to run SATSUMA in loop !
#!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage perl 1by1.pl for SATSUMA analysis # User need to set the reference multifasta file name here...2154 days ago
2141 days ago
2140 days ago
2137 days ago
Installing Porechop on Ubuntu !
...default: 1) -t THREADS, --threads THREADS Number of threads to use for adapter...he trimming of adapters from read ends --end_size END_SIZE The number of base pairs at each end of...2120 days ago
2089 days ago
Perl script to delete the adjacent repeats !
.../g #Alternatively, if you want to ignore case when comparing, you could use: $string =~ s/(.)\1/$1/gi #And finally, if you wanted to trim any number of duplicates down to a singl...1564 days ago
To convert just one specific read group to fastq
.... set -uex # The SRR BioProject number for the sequencing data. PROJECT=PRJNA257197 # The number of datasets to subselect from...rst N elements. Keep only valid SRR numbers. cat runinfo.txt | cut -f 1...s view -H all.bam echo "" echo "Number of alignments with read group...1556 days ago
1229 days ago
1456 days ago