hsinnan75 / MapCaller

MapCaller – An efficient and versatile approach for short-read alignment and variant detection in high-throughput sequenced genomes
MIT License
30 stars 5 forks source link

How do i detect fragment size automatically? #10

Closed tseemann closed 4 years ago

tseemann commented 4 years ago
-size Sequencing fragment size [default: 500, 
MapCaller can predict the fragment size automatically]

Default is 500 but i want to use auto-detection. How?

hsinnan75 commented 4 years ago

If the fragment size of the sequencing library is known, you may use the option "-size" to indicate the fragment size. Otherwise, MapCaller will predict the fragment size according to the paired alignments, but it needs at least 1000 paired alignments to infer the fragment size.

tseemann commented 4 years ago

It says the default is -size 500. So that means it is choosing 500 for me. How do I disable that so it does automatic -size ?

hsinnan75 commented 4 years ago

The default size is only used for the initial fragment size for selecting paired alignments when there are multiple pairs. Thus the fragment size will be refined periodically during read mapping. The default size or designated size given by the user work the same way. In fact, it does not matter what the value is given. MapCaller will always infer the average fragment size during the read mapping.

tseemann commented 4 years ago

Ok. Thank you for making that clear.