Samtools commands for bioinformatician !
...total+=1}END{print (total/c)*100}' ### flagstat # output file has two columns: QC-passed reads and QC-failed reads # and rows: total reads, duplicates, mapped. samtools flagstat s...1600 days ago
Bash script to alignment of short reads against reference genome !
...w.bam --- sort the BAM file, writing it to .sorted.bam. #marking PCR duplicates: sambamba markdup SRR1770413....SRR1770413.bam --- this marks reads which appear to be redundant PCR duplicates based on their read mapping p...1557 days ago
Perl script to delete the adjacent repeats !
...tively, 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 single letter, this...1548 days ago