New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4596277	198 (1.320)	10 (0.070)	8/168	NT	4.9%	coding (288/1023 nt)	lgoD	L‑galactonate oxidoreductase; L‑gulonate oxidoreductase
	?	NC_000913	= 4596288	198 (1.400)					coding (299/1023 nt)	lgoD	L‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Rejected: Frequency below cutoff threshold.

ATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4596174‑4596277 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcccgCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAAT  <  NC_000913/4596288‑4596191                                                                                                                                                                                                            ATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATG                                                                                                        <  1:790805/100‑1                       TGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCAT                                                                                  <  1:2497515/100‑1                       TGGGCTGGGTAAAAATATTCCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCAT                                                                                  <  1:669232/100‑1                         GGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTG                                                                                <  1:3348675/100‑1                         GGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTG                                                                                <  1:2034229/100‑1                          GCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTGC                                                                               <  1:507478/100‑1                           CTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTGCT                                                                              <  1:1669804/100‑1                            TGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTGCTG                                                                             <  1:3120456/100‑1                              GGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTGCTGAC                                                                           >  1:172235/1‑100                               GTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTGCTGACA                                                                          <  2:172232/100‑1                                TAAAAATATTGCTTATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCGCTTTTACACGCCGGGCATTGCTGACAG                                                                         <  1:3762459/100‑1                                                                                                    GCCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAAC     <  1:1795225/100‑1                                                                                                    GCCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAAC     >  2:1362688/1‑100                                                                                                     CCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACA    >  2:371777/1‑100                                                                                                     CCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACA    >  2:1628235/1‑100                                                                                                     CCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACA    >  2:1551835/1‑100                                                                                                     CCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACA    >  1:1290667/1‑100                                                                                                     CCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACA    >  1:3944051/1‑100                                                                                                     CCCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACA    <  1:1875724/100‑1                                                                                                      CCGCTTTTACACGCCGGGCATTGCTTACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAA   <  2:3379703/100‑1                                                                                                      CCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAA   <  1:3938662/100‑1                                                                                                      CCGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAA   <  1:1842945/100‑1                                                                                                       CGCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAAT  >  2:2452389/1‑100                                                                                                                                                                                                            ATGAAATATGTGGGGAGATTGTTGGGCTGGGTAAAAATATTGCTGATCTTAAAAATGGTCAGCAAGTTGCTGTGATCCCTTATGTTGCCTGTCAGCAATGCCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4596174‑4596277 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcccgCTTTTACACGCCGGGCATTGCTGACAGGCAACATAAGGGATCACAGCAACTTGCTGACCATTTTTAAGATCAGCAATATTTTTACCCAGCCCAACAAT  <  NC_000913/4596288‑4596191

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