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...2066 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 ##...1606 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...1562 days ago
1527 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";...428 days ago
302 days ago