This Programm verifies which 16S variable region was sequenced.
Input: .fastq file containing 16S sequencing reads or a Directory containing (or subfolders containing) fasta or fastq files Output: print(Probabilities for every region) CSV file containing read_designation, Percentage unaligned reads, Probabilities of every region
Requirements: python 3 bowtie2 samtools bedtools