Comments (5)
Hi! Apologies for the delayed answer. Indeed the FOF strategy doesn't work as is.
You'll need to
- if applicable, for each insert size, concatenate all the left (resp. right) files from that insert size into a single left (resp. right) file
- pass as argument each insert size library, e.g.
-1 left_insertA.fq.gz -2 right_insertA.gz --mp-1 left_insertB.fq.gz --mp-2 right_insertB.fq.gz
in order from smaller insert size to larger (ie in my example, insert size A is smaller than B)
Alternatively:
- run gatb-pipeline a FOF of all libraries in any order. This will produce contigs but no scaffold. Note that these will be the same contigs as if you had specified that librairies were paired/mate-pairs, as Minia does not care about pairing when making contigs. Then you can run BESST stand-alone (using the gatb-pipeline script by using the
-c
argument, or just the BESST program, or any other scaffolder) manually using the contigs produced by gatb-pipeline.
from gatb-minia-pipeline.
Thank you for the suggestion and apologies for the delayed response!
I did use a mixture of the two options though. I concatenated the smaller insert fastqs and used the mate-pair separately.
from gatb-minia-pipeline.
sounds good, did it work?
from gatb-minia-pipeline.
Yup, it did! I got a decent contiguity as well.
Got the 2.5Gb genome in 9876 (no joke) scaffolds over 1kb in length with an N50 of 1.83Mb. Having a 8Kb and 20Kb insert MP libraries did help a lot combined with dual scaffolding.
from gatb-minia-pipeline.
very nice!
from gatb-minia-pipeline.
Related Issues (20)
- make test failed and number of cores not limitable? HOT 5
- gatb-pipeline requires a Linux 64 bits system HOT 2
- Problem with BESST in the pipeline HOT 12
- Execution of 'perl tools/sga-deinterleave.pl' HOT 4
- Problems with assembling the data beyond kmer 127 and other queries. HOT 4
- hdf5 issue - unable to open file HOT 5
- Add version info on prerequisites HOT 1
- Error with `make test` HOT 3
- 920: syntax error in expression (error token is "920") HOT 1
- GFA output HOT 1
- BESST: error: unrecognized arguments: fr HOT 3
- Flag to remove intermediate files HOT 4
- Release package?
- Question about input file type HOT 2
- ZeroDivisionError while running BESST HOT 1
- BESST: error: argument -orientation is required HOT 3
- OSError: [Errno 17] File exists HOT 1
- test if bwa version is high enough
- Error with samtools in `make test` HOT 1
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 gatb-minia-pipeline.