New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1516656 =	107 (0.780)	10 (0.080)	7/212	NT	8.6%	coding (1079/1425 nt)	patD	gamma‑aminobutyraldehyde dehydrogenase
	?	NC_000913	1516683 =	113 (0.870)					coding (1106/1425 nt)	patD	gamma‑aminobutyraldehyde dehydrogenase

CCATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1516810‑1516656 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcatTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAG  >  NC_000913/1516683‑1516827                                                                                                                                                                                                                                                                                                              CCATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCG                                                                                                                                                         >  4:1026512/1‑149   ATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGC                                                                                                                                                       <  1:1099732/149‑1    TACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCG                                                                                                                                                      >  3:635606/1‑149     ACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGC                                                                                                                                                     <  1:1327323/149‑1       CGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT                                                                                                                                                   >  4:697110/1‑149                                                                                                                             TGTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGC                                                        >  3:55687/1‑122                                                                                                                             TGTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGC                                                        <  4:55687/122‑1                                                                                                                              GTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTT                                  >  3:315002/1‑143                                                                                                                              GTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTT                                  >  3:82480/1‑143                                                                                                                              GTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTT                                  <  4:315002/143‑1                                                                                                                              GTAATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTT                                  <  4:82480/143‑1                                                                                                                                 ATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCG                                                                                    >  2:901750/1‑90                                                                                                                                 ATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCG                                                                                    <  1:901750/90‑1                                                                                                                                 ATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTT                                                                                                                >  3:996269/1‑62                                                                                                                                 ATGCGCCAGCCAGCAGCGTCGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTT                                                                                                                <  4:996269/62‑1                                                                                                                                                        GCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAG  >  4:481743/1‑149                                                                                                                                                          ATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTT                                  <  2:333987/115‑1                                                                                                                                                          ATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTA                                         <  2:255266/108‑1                                                                                                                                                          ATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTA                                         >  1:255266/1‑108                                                                                                                                                          ATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAAC                                                             >  2:714239/1‑88                                                                                                                                                          ATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAAC                                                             <  1:714239/88‑1                                                                                                                                                          ATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGCGAACTGGGCGAATGACAGCCAGTACGGACTT                                  >  1:333987/1‑115                                                                                                                                                                                                                                                                                                              CCATACCGAAGATGCAAGTCCGTACTGGCTGTCATTCGCCCAGTTCACCACCTGTTCTTCGTTGTCGAAGGGCGTAACACTCACTACTGGACCAAATACCTCTTTTTGCACGATGGCATCGTCCTGTAATGCGCCAGCCAGCAGCGTCGGCGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1516810‑1516656 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcatTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAAAGATGTGGGCAG  >  NC_000913/1516683‑1516827

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ 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.