/gcbyposition

GC content by position in fastq alignment

Primary LanguagePython

#G/C Content by Base Pair

A short python script to calculate the g/c content for each position in an alignment in fastq format.

This script:

  • Reads in a gzipped fq file
  • Returns a dictionary with:
    • key: position in alignment
    • value: proportion of g/c sites at position

I made a change