Perl script to extract sequence by Ids from multifasta file !
#!/usr/bin/perl -w use strict; my $idsfile = "$ARGV[0]"; my $seqfile = "$ARGV[1]"; my %ids = (); open...while() { chomp; $ids{$_} += 1; } close FILE; local $/ = "\n>"; # read...2093 days ago
2072 days ago
2057 days ago
Generates a genome coverage plot with R
library(CoverageView) ##draw a coverage plot for a test case BAM file #get a BAM test file treatBAMfile2052 days ago
Generate simulated polyploid genome !
...sta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheterozygous3percent.fasta #generated 50X of Illumina paired-end reads sim_reads --depth 50 toyheterozygous3percent.fa...2050 days ago
2023 days ago
2001 days ago
1602 days ago
Samtools commands for bioinformatician !
...sej.bam | less ### count reads (-c) by flag specified samt...sults into txt file samtools depth /path/to/sorted_bam.bam > /p...t ### note: returns count of depth at each position. input shoul...one liner to count coverage breadth samtools depth -a sorted_du...1600 days ago
1550 days ago