2069 days ago
Generate simulated polyploid genome !
...t.fasta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheterozygous3percent.fasta #generated 50X of Illumina paired-end reads sim_reads --depth 50 toyheterozygous3pe...2047 days ago
2028 days ago
1828 days ago
1719 days ago
Samtools commands for bioinformatician !
...le into fasta samtools fasta reads_mapped.sam > reads.fasta #..._sej.bam | less ### sort reads by flag specified and show th...sej.bam | less ### count reads (-c) by flag specified samto...le has two columns: QC-passed reads and QC-failed reads # and ro...1597 days ago
Bash script to download SRA file !
...for the first and second mate in each pair. We'll use them in this format when aligning. wget ftp://ftp-trace.ncbi.nlm.nih.gov/sra/sra instant/reads/ByRun/sra/SRR/SRR177/SRR17704...1553 days ago
Bash script to alignment of short reads against reference genome !
...by line: #alignment with bwa: bwa mem -t $threads -R '@RG\tID:K12\tSM:K12' --- this says "align using so many threads" and also "give the reads the...conversion to BAM: samtools view -b - --- this reads SAM from stdin (the - specifi...1553 days ago
To convert just one specific read group to fastq
..., | grep SRR | head -$N > selected.txt # Store the data in the reads folder. mkdir -p reads # Download the SRR data fo...GROUP-SRR1972919 all.bam # Reverting the process is to extract reads, tagged with readgroups to pa...1536 days ago
Perl script to reads and extract webpage contents !
use 5.010; use strict; use warnings; use WWW::Mechanize; my ($url) = @ARGV; die "Usage: $0 URL\n" if not $url; my $w = WWW::Mechanize->new; $w->get($url); say $w->content; #Run #jit@jit-HP-Pro-3335-MT:~/Downloads/testDock$ perl web.pl https://bioinformaticsonline.com1507 days ago