New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1516656 =	90 (0.970)	7 (0.080)	7/270	11.4	7.5%	coding (1079/1425 nt)	patD	gamma‑aminobutyraldehyde dehydrogenase
	?	NC_000913	1516683 =	87 (0.980)					coding (1106/1425 nt)	patD	gamma‑aminobutyraldehyde dehydrogenase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1516810‑1516656 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcatTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAGGGCG  >  NC_000913/1516683‑1516831                                                                                                                                                                                                                                                                                                                  CCATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGG                                                                                                                                                            >  3:257382/1‑150    TACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGC                                                                                                                                                         >  2:252909/1‑150      CCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT                                                                                                                                                       <  1:871245/150‑1      CCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTCTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT                                                                                                                                                       <  1:991283/150‑1                                                                                                                      ATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCACTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACCGCCAGTACGGACT                                       >  1:526793/1‑150                                                                                                                          TCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCG                                                           <  1:379487/126‑1                                                                                                                          TCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCG                                                           >  2:379487/1‑126                                                                                                                            CTGTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATC                                 >  2:10292/1‑150                                                                                                                              GTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGAC                                        <  1:963489/141‑1                                                                                                                              GTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGAC                                        >  2:963489/1‑141                                                                                                                                       CAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGA                      >  2:45957/1‑150                                                                                                                                                       CGCATTACAGGACGATGCTATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAG      >  2:581569/1‑150                                                                                                                                                        GCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAGG     >  3:147862/1‑150                                                                                                                                                        GCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAGG     >  4:93933/1‑150                                                                                                                                                           TTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAGGGCG  >  2:452578/1‑150                                                                                                                                                                                                                                                                                                                  CCATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1516810‑1516656 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcatTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAGGGCG  >  NC_000913/1516683‑1516831

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.