BBTools for bioinformatician !
...h means N (in either fasta or fastq format). The BBMap package au...mat.sh in=reads.sam out=reads.fastq Verify pairing and optio...0 adderrors=false out=perfect.fastq reads=1000 minlength=18 maxle...omma-delimited (e.g. in=file1.fastq,file2.fastq,file3.fastq). And...2233 days ago
2096 days ago
2905 days ago
1235 days ago
2114 days ago
Genome Assembly Tools and Software - PART2 !!
...alability in mind. It can assemble large and high-coverage genomes from fastq files in a short time and pro...e-assembled contig sequences (FASTA) and NGS paired-read data (FASTA or FASTQ). The final scaffolds are pro...2648 days ago
Useful Publications and Websites for Deep Sequencing Data Analysis
...Res 21(5):734-40, 2011. Full Text Compression of DNA sequence reads in FASTQ format. Deorowicz & Grabo...equencing errors. Kelley et al, Genome Biol 11(11):R116, 2010. PubMed FastQC: a quality control tool for...3742 days ago
Perl one-liner for bioinformatician !!!
...encing data most of the bioinformaticians mung and wrangle around massive amounts of genomics text. There are several "standardized" file formats (FASTQ, SAM, VCF, etc.) and some too...3591 days ago
Extract the fastq sequence with range in Perl
use Bio::DB::Fasta; open(POSITIONS,"positions.txt"); while(){ chomp; my ($seqName,$begin,$end) = split(/\s/); my $db = Bio::DB::Fasta->new('allGenomeContacted.fa'); my $seq = $db->seq("$seqName", $begin => $end); print "$seq\n"; } close(POSITIONS);2431 days ago
Convert fastq to fasta in Perl
use Bio::SeqIO; #convert .fastq.gz to .fasta open my $zcat, 'zcat seq.fastq.gz |' or die $!; my $in=Bio::SeqIO->new(-fh=>$zcat, -format=>'fastq'); my $out=Bio::SeqIO->new(-...2293 days ago
Tags: Bioinformatics, FASTQ, Compression, Decompression, NGS
2905 days ago
Tags: Bioinformatics, FASTQ, NGS, ALE, Framework, Stats, Reads, Evaluation, Quality, Assembly
2894 days ago
Comment on "A quick guide to Phred scaling"
...quality scores are typically encoded in the quality score field of the FASTQ file format, which is a commo...h base call and are typically encoded in the quality score field of the FASTQ file format....378 days ago
Comment on "Short-read assembly using Spades !"
.... This step can be performed using various tools such as Trimmomatic or FastQC. Running SPAdes: After...is usually as follows: spades.py -o output_directory -1 forward_reads.fastq -2 reverse_reads.fastq Here,...378 days ago