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...2255 days ago
2117 days ago
2926 days ago
1256 days ago
2135 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...2670 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...3763 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...3612 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);2452 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(-...2314 days ago
Tags: Bioinformatics, FASTQ, Compression, Decompression, NGS
2926 days ago
Tags: Bioinformatics, FASTQ, NGS, ALE, Framework, Stats, Reads, Evaluation, Quality, Assembly
2915 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....399 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,...399 days ago