https://bioinformatics-core-shared-training.github.io/cruk-summer-school-2017/Day1/Session4-seqIntro.html

Address of the bookmark: https://bioinformatics-core-shared-training.github.io/cruk-summer-school-2017/Day1/Session4-seqIntro.html

Eligibility:
M.Sc. in any stream of Life Sciences with minimum 55% marks.

Desirable:
Candidates having experience in Molecular Spectroscopy, Protein Folding and Bioinformatics will be preferred.

Interested candidate may appear in the walk in Interview conducted on September 16, 2014 (Tuesday) 11:00 AM in the Director's Office, Centre for Interdisciplinary Research in Basic Sciences, lamia Millia Islamia, New Delhi-110025.
Candidates are required to bring a set of Xerox copy of their recent CV and qualifying degree (certificate/mark sheet) along with original documents. NoTA/DA will be paid.

For any further information you may e-mail to: mihassan@jmLac.in

Read more at http://jmi.ac.in/upload/advertisement/jobs_cirbs_2014september8.pdf

So far miniasm is in early development stage. It has only been tested on a dozen of PacBio and Oxford Nanopore (ONT) bacterial data sets. Including the mapping step, it takes about 3 minutes to assemble a bacterial genome. Under the default setting, miniasm assembles 9 out of 12 PacBio datasets and 3 out of 4 ONT datasets into a single contig. The 12 PacBio data sets are PacBio E. coli sample, ERS473430, ERS544009, ERS554120, ERS605484, ERS617393, ERS646601, ERS659581, ERS670327, ERS685285, ERS743109 and a deprecated PacBio E. coli data set. ONT data are acquired from the Loman Lab.

For a C. elegans PacBio data set (only 40X are used, not the whole dataset), miniasm finishes the assembly, including reads overlapping, in ~10 minutes with 16 CPUs. The total assembly size is 105Mb; the N50 is 1.94Mb. In comparison, the HGAP3produces a 104Mb assembly with N50 1.61Mb. This dotter plot gives a global view of the miniasm assembly (on the X axis) and the HGAP3 assembly (on Y). They are broadly comparable. Of course, the HGAP3 consensus sequences are much more accurate. In addition, on the whole data set (assembled in ~30 min), the miniasm N50 is reduced to 1.79Mb. Miniasm still needs improvements.

Miniasm confirms that at least for high-coverage bacterial genomes, it is possible to generate long contigs from raw PacBio or ONT reads without error correction. It also shows that minimap can be used as a read overlapper, even though it is probably not as sensitive as the more sophisticated overlapers such as MHAP and DALIGNER. Coupled with long-read error correctors and consensus tools, miniasm may also be useful to produce high-quality assemblies.

Minimap and miniasm are ultrafast tools for (i) mapping and (ii) assembly. Designed for long, noisy reads, they do not have a correction or consensus step, and therefore the resulting assemblies are contiguous (i.e. long) but very noisy (i.e. full of errors)

We start with an all against all comparison:

minimap -Sw5 -L100 -m0 -t8 reads.fq reads.fq | gzip -1 > reads.paf.gz

Then we can assemble

miniasm -f reads.fq reads.paf.gz > reads.gfa

Convert GFA to FASTA:

awk '/^S/{print ">"$2"\n"$3}' reads.gfa | fold > reads.fa

And then count how many contigs:

grep ">" reads.fa | wc -l

# Download sample PacBio from the PBcR website
wget -O- http://www.cbcb.umd.edu/software/PBcR/data/selfSampleData.tar.gz | tar zxf -
ln -s selfSampleData/pacbio_filtered.fastq reads.fq
# Install minimap and miniasm (requiring gcc and zlib)
git clone https://github.com/lh3/minimap && (cd minimap && make)
git clone https://github.com/lh3/miniasm && (cd miniasm && make)
# Overlap
minimap/minimap -Sw5 -L100 -m0 -t8 reads.fq reads.fq | gzip -1 > reads.paf.gz
# Layout
miniasm/miniasm -f reads.fq reads.paf.gz > reads.gfa

Address of the bookmark: https://github.com/lh3/miniasm

(Deemed University)

Hamdard Nagar, New Delhi – 110 062

R E C R U I T M E N T

(Advertisement No. 5/2014)

Applications on prescribed form are invited for filling up the following teaching positions in the Department of Biotechnology, Faculty of Science in the university. Eligible candidates may apply on or before 30.09.2014.

1. Professor/Associate Professor - One in Pay Band of Rs. 37400-67000+ AGP Rs.10000/9000

2. Assistant Professor - Two in Pay Band of Rs. 15600-39100+ AGP Rs. 6000/-

ASSISTANT PROFESSOR – 02 (including 01 SFS)

Specialization : Bioinformatics

Qualification and Experience :

Ph.D. in Biotechnology or an allied discipline with M.Sc. in Biotechnology/ Biochemistry in the First division or equivalent grade from a recognized University/ Institute.

NET in Life Science or allied discipline in addition to the above qualification.

Experience : At least two years of Post-doctoral teaching and/or research experience in Bioinformatics or relevant field in a UGC recognized Institution of repute or international research institute/ University. Proof of research to be evidenced by publications in SCI-indexed journals of high impact factor as the first or corresponding author.

Note : University may consider exempting candidates from NET, who has been awarded Ph.D. degree from ‘A’ Grade accredited University following the procedure as notified by the UGC in its Regulations of 2009 and adopted by Jamia Hamdard.

For more information: http://www.jamiahamdard.ac.in/PDF/Online%20application%20form%20_Teaching_1.pdf
http://www.jamiahamdard.ac.in/PDF/PBAS.pdf

Find the detail of the reads repeats:

fq2fa ONT_A.fastq ONT_A.fasta 

minimap2 -xava-ont ONT_A.fasta ONT_A.fasta -t10 -X > AONT.paf 

awk '{if($1==$6){print}}' AONT.paf > AONTself.paf 

awk '$5=="-"' AONTself.paf | awk '{print $1}'| sort|uniq > invertedrepeat.list

Generated a few palindrome and repeats plots (highlighting only repeats largest than 10, 20 and 30 kb)

minidot -f 5 -m 30000 AONTself.paf > AONTself30000.eps 
sed 's/_template_pass_FAH31515//' AONTself30000.eps > AONTself30000final.eps 

minidot -f 5 -m 20000 AONTself.paf > AONTself20000.eps 
sed 's/_template_pass_FAH31515//' AONTself20000.eps > AONTself20000final.eps 

minidot -f 5 -m 10000 AONTself.paf > AONTself10000.eps 
sed 's/_template_pass_FAH31515//' AONTself10000.eps > AONTself10000final.eps 

Assemble with miniasm:

miniasm -f ONT_A.fasta AONT.paf > AONT.gfa 

grep '^S' AONT.gfa |awk '{print ">"$2"\n"$3}' > AONT_miniasm.fasta 

minimap2 -xasm10 AONT_miniasm.fasta AONT_miniasm.fasta -t1 -X > AONT_miniasm.paf 

awk '{if($1==$6){print}}' AONT_miniasm.paf > AONT_miniasm_self.paf 

minidot -f 5 -m 10000 AONT_miniasm_self.paf > AONT_miniasm_self10000.eps 

Njoy the assembly !

Functional and Comparative Genomics

http://jp.senescence.info/
http://www.senescence.info/
http://www.liv.ac.uk/integrative-biology/staff/joao-de-magalhaes/

Extent of employment: 30 Hours per Week
Duration of employment: 1st of October 2015 to 30th of September 2019
Gross monthly salary and pay grade in terms of collective agreement for university staff (payable 14 times per year): B1, € 1.997,20

Responsibilities
The successful candidate (f/m) will pursue a Ph.D. project related to the interpretation of plant genome and transcriptome sequencing data from next-generation sequencing (NGS) platforms. In particular, the candidate will characterize the unexplored genome of quinoa, a crop plant of long-standing tradition in Latin America. We collaborate with research partners in Austria and abroad, and the candidate’s project will be of central importance in the context of this research network.

Required skills and qualifications
We are looking for a graduate student (f/m) with a Master’s degree in bioinformatics or in a related field, solid programming skills (e.g. developing sequence analysis tools), experience with the analysis of NGS data sets, understanding of lab methods and knowledge of genomics/transcriptomics. The group has successfully performed several projects using NGS technology. We have recently published the reference genome sequence of sugar beet (Dohm et al., Nature, 2014), a crop plant closely related to quinoa (same family, but different genus). Not yet published is a quinoa genome assembly that we have generated, and which will serve as the starting point of the candidate’s project. We are a multidisciplinary team and offer work in a lively and friendly atmosphere, and state-of-the-art computing infrastructure. We are looking forward to expanding our team by a dedicated and strongly motivated person with a distinct interest in the challenges of plant genomics.

Applications can be submitted until: 16th of August 2015

University of Natural Resources and Life Sciences Vienna seeks to increase the number of its female faculty and staff members. Therefore qualified women are strongly encouraged to apply. In case of equal qualification, female candidates will be given preference unless reasons specific to an individual male candidate tilt the balance in his favour.

Please send your job application (incl. letter of motivation, CV, summary of Master’s thesis and contact details for two referees) to Personnel department, University of Natural Resources and Life Sciences, 1190 Vienna, Peter-Jordan-Straße 70; E-Mail: kerstin.buchmueller@boku.ac.at. (Reference code: 59)

We regret that we cannot reimburse applicants travel and lodging expenses incurred as part of the selection and hiring process.

www.boku.ac.at

bwa aln hs37m.fa r1.fq > r1.sai && bwa aln hs37m.fa r2.fq > r2.sai \
&& bwa sampe hs37m r1.sai r2.sai r1.fq r2.fq > r12.bwa.sam

bwa bwasw ../index/bwa/hs37m.fa r12.fq > r12.bwasw.sam

gsnap -A sam -d hs37m r1.fq r2.fq > r12.gsnap.sam

novoalign -r Random -o SAM -f r1.fq r2.fq -i 500 50 -d hs37m-k14s3.novo > r12.novo.sam

smalt map -f samsoft -i 650 -o r12.smalt-k20s13.sam hs37m-k20s13 r1.fq r2.fq

stampy.py -g hs37m -h hs37m -o r12.stampy.sam -M r1.fq,r2.fq

soap -D hs37m.fa.index -a r1.fq -b r2.fq -l 32 -g 3 -u dummy -2 dummy -o r12.soap

• Updated human gene set (GENCODE 21)
• Updated rat gene set including manual annotation from HAVANA
• New species: Vervet-African green monkey
• Imported Transcript Support Levels (TSLs) from UCSC for human and mouse
• Imported APPRIS flag for human and mouse
• Updated Amazon molly gene set

Find more at http://www.ensembl.info/blog/2014/10/02/ensembl-77-has-been-released/