SNPs detection by eBWT positional clustering

Marinella Sciortino, Nicola Prezza, Giovanna Rosone, Nadia Pisanti

Risultato della ricerca: Articlepeer review

18 Citazioni (Scopus)


Background: Sequencing technologies keep on turning cheaper and faster, thus putting a growing pressure for data structures designed to efficiently store raw data, and possibly perform analysis therein. In this view, there is a growing interest in alignment-free and reference-free variants calling methods that only make use of (suitably indexed) raw reads data. Results: We develop the positional clustering theory that (i) describes how the extended Burrows-Wheeler Transform (eBWT) of a collection of reads tends to cluster together bases that cover the same genome position (ii) predicts the size of such clusters, and (iii) exhibits an elegant and precise LCP array based procedure to locate such clusters in the eBWT. Based on this theory, we designed and implemented an alignment-free and reference-free SNPs calling method, and we devised a consequent SNPs calling pipeline. Experiments on both synthetic and real data show that SNPs can be detected with a simple scan of the eBWT and LCP arrays as, in accordance with our theoretical framework, they are within clusters in the eBWT of the reads. Finally, our tool intrinsically performs a reference-free evaluation of its accuracy by returning the coverage of each SNP. Conclusions: Based on the results of the experiments on synthetic and real data, we conclude that the positional clustering framework can be effectively used for the problem of identifying SNPs, and it appears to be a promising approach for calling other type of variants directly on raw sequencing data. Availability: The software ebwt2snp is freely available for academic use at:
Lingua originaleEnglish
pagine (da-a)1-13
Numero di pagine13
RivistaAlgorithms for Molecular Biology
Stato di pubblicazionePublished - 2019

All Science Journal Classification (ASJC) codes

  • ???subjectarea.asjc.1300.1315???
  • ???subjectarea.asjc.1300.1312???
  • ???subjectarea.asjc.1700.1703???
  • ???subjectarea.asjc.2600.2604???


Entra nei temi di ricerca di 'SNPs detection by eBWT positional clustering'. Insieme formano una fingerprint unica.

Cita questo