History

Our Sponsors



Download BioinformaticsOnline(BOL) Apps in your chrome browser.

Freelancer

AdvertiseWithUs

My most common one-liner in bioinformatics: Revision

script available for doing sequence length disrtibution of fastq file

cat reads.fastq | awk '{if(NR%4==2) print length($1)}' | sort -n | uniq -c > read_length.txt

Plot syntenty

/home/urbe/Tools/Sibelia-3.0.6-Linux/bin/Sibelia -s far --gff -q -o palCh101 ch101_read33155_template_pass_FAH31515.fa

SatsumaSynteny 

/home/urbe/Tools/SATSUMA/satsuma-code-0/SatsumaSynteny -t ../allPalindromeSimulated.fa -q MIRAsimulation_out.padded_CORRECTED.fasta -n 2 -m 8 -o output_directory

MicrosyntenyPlot: Go inside the outfolder of Satsuma synteny and run this 

/home/urbe/Tools/SATSUMA/satsuma-code-0/MicroSyntenyPlot -s 8000 -i xcorr_aligns.seeds.out -o seeds.ps

And to quickly obtain a graph in R:

reads<-read.csv(file="read_length.txt", sep="", header=FALSE)
plot (reads$V2,reads$V1,type="l",xlab="read length",ylab="occurences",col="blue")