Estimate Genome Size with Jellyfish and R
jellyfish count -t 8 -C -m 19 -s 5G -o 19mer_out --min-qual-char=? /common/Tutorial/Genome_estimation/sample_read_1.fastq /common/Tutorial/Genome_estimation/sample_read_2.fast...2248 days ago
2212 days ago
2212 days ago
Perl script to find the distance beetween all the contigs and scaffolds
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; $| = 1; #Script...e distance beetween all the contigs and scaffolds #Usage: perl clustalReads.pl genome.fa > HammingDist.txt #Depend...2142 days ago
Perl script to run SATSUMA in loop !
#!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage per...e reference multifasta file name here my $seqfile=""; my $queryfile = "genome.fasta"; # Ur query genome my $tarfile = "renamedAdinet...2142 days ago
Perl script to convert GFF 2 FASTA !
...my $outfile_pep = Bio::SeqIO->new( -format => 'fasta', -file =>...y $outfile_gene = Bio::SeqIO->new( -format => 'fasta', -file =>...oter) ### First, index the genome my $file_fasta = $ARGV[0]; my $db = Bio::DB::Fasta->new($file_fasta); print ("Genome fasta parsed\n"); ###...2141 days ago
2139 days ago
Perl subroutine to read genome/reads fasta file !
sub readSeqFromFasta{ my $file = $_[0]; my (%map,$id,$seq); open(IN,$file) or die("Cannot open file for reading $file:$!\n"); while(){ chomp; if($_ =~ /^>...2135 days ago
2128 days ago
2077 days ago