After meeting with Steven this week, we decided the best course of action was to revise the analysis pipeline and document why we chose the settings we did. The first part of this revision is to go back to
bismark. When I ran my full samples through
bismark, I used the default setting for
-score_min. This lead to a ~20% mapping efficiency. That’s not great.
-score_min option dictates how strict the alignment is. The default option is L,0,-0.2, which is pretty stringent. In this Jupyter notebook, I tested 3 different
-score_min options: L,0,-0.6; L,0,-0.9; L,0,-1.2. Samples were run in the following order:
Therefore, the first mapping efficiency listed belongs to sample 10, second to sample 1, etc. I kept this in mind and generated the following table:
Table 1. Mapping efficiency (%) based on different