BLASTPLOT: a PERL module to plot next generation sequencing NCBI-BLAST results

BACKGROUND: The development of Next Generation Sequencing (NGS) during the last decade has created an unprecedented amount of sequencing data, as well as the ability to rapidly sequence specimens of interest. Read-based BLAST analysis of NGS data is a common procedure especially in the case of metagenomic samples. However, coverage is usually not enough to allow for de novo assembly. This type of read-based analysis often creates the question of how the reads that align to the same sequence are distributed. The same question applies to preparation of primers or probes for microarray experiments. Although there are several packages that allow the visualization of DNA segments in relation to a reference, in most cases they require the visualization of one reference at a time and the capture of screen shots for each segment. Such a procedure could be tedious and time consuming. The field is in need of a solution that automates the capture of coverage plots for all the segments of interest. RESULTS: We have created BLASTPLOT, a PERL module to quickly plot the BLAST results from short sequences (primers, probes, reads) against reference targets. CONCLUSIONS: BLASTPLOT is a simple to use PERL module that allows the generation of PNG graphs for all the reference sequences associated with a BLAST result set

Finishing genomes with limited resources: lessons from an ensemble of microbial genomes

While new sequencing technologies have ushered in an era where microbial genomes can be easily sequenced, the goal of routinely producing high-quality draft and finished genomes in a cost-effective fashion has still remained elusive. Due to shorter read lengths and limitations in library construction protocols, shotgun sequencing and assembly based on these technologies often results in fragmented assemblies. Correspondingly, while draft assemblies can be obtained in days, finishing can take many months and hence the time and effort can only be justified for high-priority genomes and in large sequencing centers. In this work, we revisit this issue in light of our own experience in producing finished and nearly-finished genomes for a range of microbial species in a small-lab setting. These genomes were finished with surprisingly little investments in terms of time, computational effort and lab work, suggesting that the increased access to sequencing might also eventually lead to a greater proportion of finished genomes from small labs and genomics cores

Characterization of pPCP1 Plasmids in Yersinia pestis Strains Isolated from the Former Soviet Union

Complete sequences of 9.5-kb pPCP1 plasmids in three Yersinia pestis strains from the former Soviet Union (FSU) were determined and compared with those of pPCP1 plasmids in three well-characterized, non-FSU Y. pestis strains (KIM, CO92, and 91001). Two of the FSU plasmids were from strains C2614 and C2944, isolated from plague foci in Russia, and one plasmid was from strain C790 from Kyrgyzstan. Sequence analyses identified four sequence types among the six plasmids. The pPCP1 plasmids in the FSU strains were most genetically related to the pPCP1 plasmid in the KIM strain and least related to the pPCP1 plasmid in Y. pestis 91001. The FSU strains generally had larger pPCP1 plasmid copy numbers compared to strain CO92. Expression of the plasmid's pla gene was significantly (P ≤ .05) higher in strain C2944 than in strain CO92. Given pla's role in Y. pestis virulence, this difference may have important implications for the strain's virulence

Field-induced polarisation of Dirac valleys in bismuth

Electrons are offered a valley degree of freedom in presence of particular lattice structures. Manipulating valley degeneracy is the subject matter of an emerging field of investigation, mostly focused on charge transport in graphene. In bulk bismuth, electrons are known to present a threefold valley degeneracy and a Dirac dispersion in each valley. Here we show that because of their huge in-plane mass anisotropy, a flow of Dirac electrons along the trigonal axis is extremely sensitive to the orientation of in-plane magnetic field. Thus, a rotatable magnetic field can be used as a valley valve to tune the contribution of each valley to the total conductivity. According to our measurements, charge conductivity by carriers of a single valley can exceed four-fifth of the total conductivity in a wide range of temperature and magnetic field. At high temperature and low magnetic field, the three valleys are interchangeable and the three-fold symmetry of the underlying lattice is respected. As the temperature lowers and/or the magnetic field increases, this symmetry is spontaneously lost. The latter may be an experimental manifestation of the recently proposed valley-nematic Fermi liquid state.Comment: 14 pages + 5 pages of supplementary information; a slightly modified version will appear as an article in Nature physic