Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [long-reads]

Questions relating to long reads, such as pacbio or nanopore reads.

3
votes
1answer
41 views

Error correction within the long read

I am going to get some data from plasmid sequencing to identify SNPs on the plasmids. What is done in the lab is the following: The plasmids are purified by size. We amplify the plasmids using the ...
1
vote
1answer
33 views

Using ONT MinION, why is that we cannot get a full length DNA read?

Here, https://www.youtube.com/watch?v=E9-Rm5AoZGw, ONT claims they can read non-fragmented DNA and that with the hairpin like structure the double strands can be read at once. Then why is that MinION ...
3
votes
1answer
62 views

How to simulate nanopore reads?

I have looked already here: Tools for simulating Oxford Nanopore reads . This doesn't answer my question, because it lists a few Nanopore read simulators, but I have specific problems with each of ...
3
votes
1answer
63 views

How to map short sequences to long reads, recovering all multiply-mapped high-quality matches

The dilemma: I have a problem where I have around one million short sequences (21 bp to several 100s of basepairs) for which I need to identify all occurrences of in 20-30x coverage noisy long reads (...
2
votes
1answer
181 views

Can I pilon-polish long reads with Illumina short reads to improve structural variant detection?

I have pacbio Sequel data at 50 x coverage for a strain of animal. I would like to find structural variants compared to the reference genome sequence. At the moment, I align my reads to the reference ...
0
votes
0answers
32 views

Canu failed with 'didn't find any mers?'

I have installed canu using conda in an environment containing Python 3 on a Mac OS X. I checked that canu is properly installed by pulling out the help and it works. I launched a canu assembly job by ...
3
votes
2answers
240 views

Can I stop my nanopore sequencing run if there are no more reads being produced?

I am sequencing a whole genome on MinION. I have used this flow cell for 6 hours for the phage lambda control. That gave me around 300,000 reads. I started a 48h sequencing run on the same flow cell ...
2
votes
1answer
60 views

How to interpret the SMRT Link base modification algorithm output?

I have got PacBio data from a Sequel machine at a coverage > 50X. I have run it in SMRT Link with their Base Modification Detection algorithm. I am now trying to figure out how to interpret the output ...
3
votes
1answer
331 views

How many reads has my sequencing run produced on minion?

I am running a 48 h sequencing protocol on a FLO-MIN107, on MinION, with DNA library-prepped with SQK-RAD004. MinKNOW is installed on my Mac and controlling the MinION. I have found the following ...
1
vote
1answer
123 views

What to do with the configuration test cell after configuring the MinION with it?

What to do with the configuration test cell after configuring the MinION with it? I received a MinION with a configuration test cell in it. After running a configuration in the MinKNOW software while ...
3
votes
1answer
250 views

Does MinKNOW work with Mac OSX high sierra 10.13.1?

I have done a compatibility check as described here on my computer to know if I can use MinKNOW on my mac laptop. It says that my OS is incompatible and that the OSX must be Yosemite or El Capitan. ...
2
votes
1answer
62 views

Scaffolding a genome with hybrid data

I am assembling a ~500MB genome, and have ~150x long reads and ~200x 150bp PE short reads, with a ~400bp insert size. I've done a lot of work with minimap+miniasm, and have what I think is a good set ...
5
votes
1answer
233 views

Coverage calculation: long reads (RNA-seq)

Say your aim is to calculate the coverage of an RNA-seq experiment generated with long-read sequencing (so, uneven read length). Up to now, I relied on the Lander/Waterman equation: $$C = L*N / G$$...
2
votes
1answer
100 views

Mark BWA-SW split alignments in output for long reads

Is there a way to remove split alignments of single-end long reads from BWA-SW output? BWA-MEM has an option to include or flag split alignments in the result, but it seems BWA-SW method include them ...
6
votes
1answer
554 views

How can I use Nanopore reads to close gaps or resolve repeats in a short-read assembly?

Low coverage MinION reads should be useful to close gaps and resolve repeats left by short-read assemblers. However, I haven't had any success with the software I know about. I'm aware of the ...
4
votes
1answer
115 views

Detecting structural variants with MinION data

Working on various cancers I have an interest in detecting structural variation (SV) in human, we've successfully used various tools like Pindel, SVDetect, Manta, and LUMPY, to name a few for ...
12
votes
2answers
79 views

Finding the location and unit length of repetitive sequences within a long read

After discovering a few difficulties with genome assembly, I've taken an interest in finding and categorising repetitive DNA sequences, such as this one from Nippostrongylus brasiliensis [each base is ...
5
votes
2answers
112 views

Structural variant calling for low-coverage PacBio data

PacBio is selling ~10x PacBio SEQUEL long reads as an upgrade to Illumina data for SV discovery. In a clinical setting, the main requirements are proper sensitivity and specificity but also the ...
4
votes
1answer
112 views

Rapid metagenomics classifiers on long read data [closed]

I recently used the minION (Nanopore, 9.4 flow cell, RAD001 kit) to generate a metagenome out of environmental samples. Passed reads weren't brilliant (196, average 1,594bp lenght), but working with ...
15
votes
1answer
313 views

How can I improve a long-read assembly with a repetitive genome?

I'm currently trying to assembly a genome from a rodent parasite, Nippostrongylus brasiliensis. This genome does have an existing reference genome, but it is highly fragmented. Here are some ...
14
votes
2answers
294 views

How to deal with heterozygosity during polishing of genome assembly based on long reads?

All the long-read sequencing platforms are based on single-molecule sequencing which causes higher per-base error rates. For this reason a polishing step was added to genome assembly pipelines - ...