2113 days ago
Perl script to find palindromic pattern !
#!/usr/bin/perl use strict; use warnings; use strict; my %invert; @invert{ qw[ A C G T ] } = qw[ T G C A ]; m...s), substr( $in, $p1-$pals, ($pals+1)*2 ), $p1-$pals; } } __DATA__...2113 days ago
2089 days ago
2089 days ago
2048 days ago
2018 days ago
CollectGcBiasMetrics.jar will generate a GC bias plot for each contig
samtools index aln-pe.mapped.sorted.bam for i in $(samtools view -H aln-pe.mapped.sorted.bam | awk -F"\t" '/@SQ...n-pe.mapped.sorted.$i.bam; java -Xmx2g -jar $(which CollectGcBiasMetrics.jar) R=data/C...2018 days ago
1567 days ago
Perl script to remove duplicated lines !
#!/usr/bin/perl use strict; use warnings; { $_ = ; my $next_line; while( $next_line = ) {..._line; } continue { $_ = $next_line; } print $_ if eof; } __DATA__...1574 days ago
Bash script to download SRA file !
#We can use the sratoolkit to directly pull the sequence data (in paired FASTQ format) from the archive. fast...the SRA toolkit. It allows directly downloading data from a particular sequencing...tq-dump --split-files SRR1770413 #SRA stores data i...1573 days ago