if you map reads to GRCh38 or hg38, use the following:

ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/000/001/405/GCA_000001405.15_GRCh38/seqs_for_alignment_pipelines.ucsc_ids/GCA_000001405.15_GRCh38_no_alt_analysis_set.fna.gz

There are several other versions of GRCh37/GRCh38. What’s wrong with them? Here are a collection of potential issues:

More at http://lh3.github.io/2017/11/13/which-human-reference-genome-to-use

Address of the bookmark: http://lh3.github.io/2017/11/13/which-human-reference-genome-to-use

Our present objectives is to identify and isolate the substances secreted by Spalax cells, resolve with which components they interact that are active only on cancer cells, in order to unravel the biological mechanisms and pathways that evolved in Spalax cell machinery and ultimately lead to the death of cancer-cells. The study could attest to be a breakthrough in cancer research, using the long lived, hypoxia- and cancer-tolerant Spalax as a significant biological resource for biomedical research that hopefully could open new horizons in treatment and prevention of cancer in humans.

Contact: The applications should be submitted, together with extended CV and bibliography, summary of past accomplishments, and contact information of 3 referees, to Prof of Research Aaron Avivi (aaron@research.haifa.ac.il) AND Dr. Imad Shams (imadshams@gmail.com). (http://bit.ly/1lywShk) aaron@research.haifa.ac.il

More at http://evolution.haifa.ac.il/index.php/29-people/personal-websites/77-personal-site-avivi

Genome assembly projects typically run multiple algorithms in an attempt to find the single best assembly, although those assemblies often have complementary, if untapped, strengths and weaknesses. We present our metassembler algorithm that merges multiple assemblies of a genome into a single superior sequence. 

Address of the bookmark: https://sourceforge.net/projects/metassembler/?source=directory

The Role

We are looking for a Software Engineer - Computational Biologist with 3+
years' experience in a similar role to design and implement a backend system
to support an online individualised genomics interface. This position is a
great opportunity for an ambitious, self-motivated individual to work in a
demanding, challenging and interesting role.

Position Description:
. Participate in planning, design, and implementation of Equinome back
end systems and technologies.
. Implement interfaces and management tools for back end services.
. Manage, analyse, interpret and visualise large genomics data sets.
. Work closely with scientific team to develop new features and
application enhancements
. Design, develop and manage a genomics research database.

Qualification/Experience:
. Minimum MSc in Computer Science, Genetics, Bioinformatics or in a
related field (A Ph.D qualification would be an advantage).
. Proven 3+ years of experience in similar role.
. Highly proficient in Python, SQL, MySQL.
. Excellent knowledge of mammalian genomics, bioinformatics and
statistical/population genetics.
. Hands-on experience working with large data sets.
. Experience with front-end technologies (HTML/CSS/Javascript) an
advantage.
. Experience in rapid web application development: e.g. Django.
. Knowledge or experience of Unix Scripting and R statistical
programming language would be an advantage.
. Ability to work with minimum supervision to deliver high-quality
code on time.
. Fluency in English and good written and communication skills.
. Meticulous attention to detail.

Applications should be submitted before Friday, 26 September 2014 using the
following link:
http://bit.ly/WgbhxS

Note: Full information and application procedure is available at this link:
http://bit.ly/WgbhxS

The tool should be compatible with most UNIX flavors and has been successfully tested on the following operating systems:

• Mac OS X 10.11.1
• Mac OS X 10.10.3
• Ubuntu 14.04 LTS

Address of the bookmark: https://github.com/mculinovic/EAGLER

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

I'm running "OPERA-LG_v2.0.5/bin/preprocess_reads.pl" and have the following error:

fail to open file './temporarySam'

[bwa_aln_core] write to the disk... 0.09 sec
[bwa_aln_core] 70778880 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 161.35 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 70989574 sequences have been processed.
[main] Version: 0.7.15-r1140
[main] CMD: bwa aln -t 30 all_p_ctg.fa -
[main] Real time: 2402.523 sec; CPU: 53429.488 sec
[E::hts_open_format] Failed to open file temporarySam
samtools sort: can't open "temporarySam": No such file or directory
[bwa_aln_core] convert to sequence coordinate... 1.00 sec
[bwa_aln_core] refine gapped alignments... 6.07 sec
[bwa_aln_core] print alignments... PREPROCESS:
Fastq format is recognized
[Thu Jun 14 18:16:47 2018] Building bwa index...
bwa index -p all_p_ctg.fa /home/urbe/Tools/OPERA-LG_v2.0.6/all_p_ctg.fa
[Thu Jun 14 18:18:35 2018] Finding the SA coordinates of the reads using BWA aln...
[Thu Jun 14 18:58:37 2018] Generate alignments of reads using bwa sampe...
bwa samse -n 1 all_p_ctg.fa read.sai - | grep '\(^@\|XT:A:U\)' | /usr/local/bin/samtools view -S -h -b -F 0x4 - | /usr/local/bin/samtools sort -@ 20 -no - temporarySam > FALCON-Unzip-Scaff.bam
Mapping long-reads using blasr...
/home/urbe/Tools/SSpace/SSPACE-LongRead_v1-1/blasr -nproc 40 -m 1 -minMatch 5 -bestn 10 -noSplitSubreads -advanceExactMatches 1 -nCandidates 1 -maxAnchorsPerPosition 1 -sdpTupleSize 7 /media/urbe/MyDDrive/ONTdata/allONT/allONT.fasta /home/urbe/Tools/OPERA-LG_v2.0.6/all_p_ctg.fa | cut -d ' ' -f1-5,7-12 | sed 's/ /\t/g' > FALCON-Unzip-Scaff.map
sh: 1: /home/urbe/Tools/SSpace/SSPACE-LongRead_v1-1/blasr: Permission denied
Sorting mapping results...
sort -k1,1 -k9,9g FALCON-Unzip-Scaff.map > FALCON-Unzip-Scaff.map.sort
Analyzing sorted results...
Extracting linking information...
i3 2000 5000
i2 1000 2000
i4 5000 15000
i0 -200 300
i5 15000 40000
i1 300 1000
Repeat detection...
/home/urbe/Tools/OPERA-LG_v2.0.6/bin//filter_conflicting_edge.pl pairedEdges_i0 contig_length.dat 100 2
Illegal division by zero at /home/urbe/Tools/OPERA-LG_v2.0.6/bin//filter_conflicting_edge.pl line 93.
readline() on closed filehandle FILE at bin/OPERA-long-read.pl line 250.
rm anchor_contig_info.dat contig_length.dat filtered_edges.dat filtered_edges_cov.dat *.sai
rm: cannot remove 'anchor_contig_info.dat': No such file or directory
mv FALCON-Unzip-Scaff.bam FALCON-Unzip-Scaff-with-repeat.bam
/home/urbe/Tools/OPERA-LG_v2.0.6/bin//filter_repeat.pl FALCON-Unzip-Scaff-with-repeat.bam repeat.dat | /usr/local/bin/samtools view - -h -S -b > FALCON-Unzip-Scaff.bam
rm FALCON-Unzip-Scaff-with-repeat.bam
/home/urbe/Tools/OPERA-LG_v2.0.6/bin/OPERA-LG config > log
Analyzing 1 library: FALCON-Unzip-Scaff.bam
min library mean : 0
minimum contig length is 500
Current library: 1 out of 7
Analyzing file: pairedEdges_no_repeat_i0
Analyzing file: pairedEdges_no_repeat_i1
Analyzing file: pairedEdges_no_repeat_i2
Analyzing file: pairedEdges_no_repeat_i3
Analyzing file: pairedEdges_no_repeat_i4
Analyzing file: pairedEdges_no_repeat_i5
ln -s results/scaffoldSeq.fasta scaffoldSeq.fasta

To resolve this, try downloading blasr version 1.3 above and re-run :)

(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

Functional and Comparative Genomics

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

This script has been modified by others to convert InDels as well, e.g. this by David Eccles

./vcf2fq.pl -f <input.fasta> <all-site.vcf> > <output.fastq>

https://github.com/gringer/bioinfscripts/blob/master/vcf2fq.pl

https://github.com/lh3/samtools/blob/master/bcftools/vcfutils.pl