Perl script to break the contigs by 'N'
...icating that start and #the end of the NNNNN region in the contig my @regions=(); my @bases=split //,$seq; push(@bases,"E"); #This extra position marks the end of the contig...2055 days ago
Samtools commands for bioinformatician !
...mtools view -q 42 -c sal_sej.bam ### sorting bam file by genome position samtools sort sal_sej.bam >...path/to/coverage_results.txt ### note: returns count of depth at each position. input should be sorted ##...1596 days ago
Bash script to alignment of short reads against reference genome !
...mba markdup SRR1770413.raw.sorted.bam SRR1770413.bam --- this marks reads which appear to be redundant PCR duplicates based on their read mapping position. It uses the same criteria fo...1552 days ago
1517 days ago
Raku script to find palindrome in genomes !
...{ my $substring = $dna.substr($pos, $length); if is-palindrome($substring) { say "Palindrome found at position $pos: $substring";...417 days ago
292 days ago