exon-benchmarks 0.3.4-beta.5

# Script that given a path to a FASTA will count the number
# of sequences that do not start with M

import gzip
import sys

from Bio import SeqIO

def main():
    fasta_path = sys.argv[1]
    count = 0

    if fasta_path.endswith('.gz'):
        fasta_path = gzip.open(fasta_path, 'rt')
    else:
        fasta_path = open(fasta_path, 'rt')

    for record in SeqIO.parse(fasta_path, 'fasta'):
        if record.seq[0] != 'M':
            count += 1
    print(count)

main()