NBISweden / GenErode

GitHub repository for GenErode, a Snakemake pipeline for the analysis of whole-genome sequencing data from historical and modern samples to study patterns of genome erosion.
GNU General Public License v3.0
21 stars 7 forks source link

filter_mpile.py does not keep first and last nucleotide in a read #68

Open verku opened 6 months ago

verku commented 6 months ago

The filter_mpile.py script removes (prints “N”) the first nucleotide in the “read” because they are accompanied by “^]” (read start^ and a character for mapping quality). The script also always removes (prints “N”) the last nucleotide in the “read” because they are accompanied by “$” (read end$).