tanghaibao / jcvi

Python library to facilitate genome assembly, annotation, and comparative genomics
BSD 2-Clause "Simplified" License
721 stars 186 forks source link

Why scaffolds with 1 markers could be oriented? #105

Open c2997108 opened 5 years ago

c2997108 commented 5 years ago

Hi,

When I used ALLMAPS with my data, I got the following summary.txt.

*** Summary for consensus map ***
===============================================================================================
                         o               Anchored               Oriented               Unplaced
-----------------------------------------------------------------------------------------------
          Markers (unique)                  1,907                  1,092                     23
            Markers per Mb                    2.7                    2.8                    0.1
             N50 Scaffolds                     33                     20                      1
                 Scaffolds                    168                     62                244,852
   Scaffolds with 1 marker                     28                      2                      5
  Scaffolds with 2 markers                     18                      2                      0
  Scaffolds with 3 markers                     14                      3                      0
Scaffolds with >=4 markers                    108                     55                      2
               Total bases    702,152,522 (80.9%)    390,221,854 (44.9%)    166,293,292 (19.1%)
-----------------------------------------------------------------------------------------------

This result contains 2 oriented scaffolds with 1 marker. I think scaffolds which have only 1 markers could not be oriented. This is bug or not?

Thanks a lot!

Best, Kazutoshi

ptranvan commented 5 years ago

Any answers ?