Comments (7)
Sorry @aallahyar about word "Subreads".I thought long reads could be called that, but I won't be calling them that anymore.
In fact, this lengthy read contained multiple contig, and I wanted to locate multiple contig through this read.So it's a long enzymatic-connection read.Your method is worth trying.I also ran a minimap2 : )
OK, I have a better picture now. Then you are saying that the read itself is concatenation of multiple "fragments" that are produced through enzyme digestion? Then you can scan the read and split it once you found the enzyme recognition sequence. Then map each individual fragments. Due to Nanopore sequencing errors, this is not perfect. But it works much better than mapping the entire read at once.
I hope it helps.
from bwa.
I got the same error.
I wonder if you managed to figure out why this happens?
from bwa.
I got the same error.
I also wonder if you managed to figure out why this happens...
from bwa.
I got the same error.
I also wonder if you managed to figure out why this happens...
Hey @dexon9109
It has been more than 3 years, so I have a vague memory of this. This could have been because I was trying to map a very long sequence (>10kb) and the system specifications (or BWA itself) could not handle it.
Please see if you are doing so and report back. It might be of use to somebody else in the future.
from bwa.
Hello @aallahyar , I did match a longer contig to nanopore subreads.And this error didn't throw a specific sequence or position, so I couldn't really locate the specific sequence.The commands to run are as follows:
bwa bwasw -b 5 -q 2 -r 1 -T 15 -z 10 -t 30 rg.fa nanopore.pass.fq.gz -f test.bam
from bwa.
I am a bit confused as nanopore does not have subreads? PacBio has.
In case the long read is the problem, I can also share that my read had a very low complexity i.e. ratio between frequency of observed A, T, G and C.
Unfortunately, I am not the developer, so I can not help further. My suggestion is to chop up your reads to smaller pieces step by step and see when you can successfully map. I would be surprised if you really need such a lengthy read to successfully find a location in the genome. Of course that is not a perfect solution, sorry about that.
Good luck!
from bwa.
Sorry @aallahyar about word "Subreads".I thought long reads could be called that, but I won't be calling them that anymore.
In fact, this lengthy read contained multiple contig, and I wanted to locate multiple contig through this read.So it's a long enzymatic-connection read.Your method is worth trying.I also ran a minimap2 : )
from bwa.
Related Issues (20)
- how to eliminate the effects of estimated insertion size distribution HOT 1
- bwa mem hangs after a few thousand reads
- Did bwa index use the SA-IS algorithm?
- Fail to read header from fasta input
- When i use bwa for mapping with grch37.p13.fa and hg19.fa,there exists some differences in some regions.
- Build fails with LTO
- [M::bwa_idx_load_from_disk] read 0 ALT contigs
- Indexing large genome failed HOT 3
- no ID within the read group line
- Could anyone tell me the exact meaning of '-5SP'?
- missing `bwa-0.7.18.tar.bz2` for 0.7.18 release HOT 1
- Is there MD tag in output of bwa?
- Set min insert size with -I option. Why do I still get proper pairs with insert size < min insert size?
- bwa alignment of Nanopro sequences to call the transposable elements (TEs)
- How are ambiguity charaters handled in BWA?
- the significance of the -V parameter and XR tag
- BWA and short read mapping
- Continuation after error "paired reads have different names" HOT 1
- How can i obtain all the regions reads hit. (in duplication region) HOT 2
- unexpected high quality scores with chimeric SA alternatives reported HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bwa.