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...2265 days ago
2128 days ago
2937 days ago
1267 days ago
2146 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...2680 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...3774 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...3623 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);2463 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(-...2325 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....410 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,...410 days ago