pysam | Python module for reading and manipulating SAM | Genomics library

from pysam import FastxFile


fasta_q_file = "INPUT.fasta"
out_filename = "OUTPUT_NAME.fasta"
size_size_trim = 50


with FastxFile(fasta_q_file) as fh, open(out_filename, mode='w') as fout:
    for entry in fh:
        sequence_id = entry.

import os
import sys 
import pysam
from Bio import SeqIO, Seq, SeqRecord

def main(in_file):
    out_file = "%s.fa" % os.path.splitext(in_file)[0]
    with open('test_non_reverse.txt', 'w') as non_reverse:
        with open(out_file, "w")