I need to perform quality control on ancient DNA FASTQ/FASTA files, including read length distribution, GC content, and damage pattern analysis. Please generate Python or bash scripts and Snakemake rules to automate these QC checks and summarize the results.