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...2314 days ago
2177 days ago
2986 days ago
1316 days ago
2195 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...2729 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...3823 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...3672 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);2512 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(-...2374 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....459 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,...459 days ago