New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 664687	1 (0.030)	23 (0.670)	19/144	0.9	94.0%	intergenic (‑78/‑48)	gltL/insJ‑2	ATP‑binding protein of glutamate/aspartate transport system; b0652_1/IS150 hypothetical protein
	?	REL606	666134 =	2 (0.060)					pseudogene (354/387 nt)	gltL	ATP‑binding protein of glutamate/aspartate transport system; b0652_1

TCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAT  >  REL606/666062‑666205                                                                         |                                                                        ctgctgaccgccggaaagctgcgccggaaacttattggcgtgagcagaaagcccgacacgctccagcagtttCAgg                                                                      <  1:420868‑M2/4‑1 (MQ=255)                   ctgcgccggaaacttattggcgtgagcagaaagcccgacacgctccagaagtttCAGGGCTTTTTCACGAGCCGgc                                                    >  1:671953‑M2/55‑76 (MQ=255)                       gccggaaacttattggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCtt                                                <  2:855939‑M2/26‑1 (MQ=255)                        ccggaaacttattggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCttt                                               <  2:768293‑M2/27‑1 (MQ=255)                        ccggaaacttattggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCttt                                               <  1:369159‑M2/27‑1 (MQ=255)                        ccTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt                                               <  1:428265/76‑1 (MQ=255)                          tataTGCCTCGTGTTTCACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTAt                                             >  2:389235/1‑76 (MQ=255)                                  attggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTCATCGCGTTTa                                     <  1:159939‑M2/37‑1 (MQ=255)                                  attggcgtgagcagaaagcccgacacgctccagcagtgtCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTa                                     <  1:804860‑M2/37‑1 (MQ=255)                                   ttggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTaa                                    >  2:764794‑M2/39‑76 (MQ=255)                                    tggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTGAg                                   >  1:549596‑M2/38‑76 (MQ=255)                                     ggcgtgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGc                                  <  1:692573‑M2/40‑1 (MQ=255)                                         tgagcagaaagcccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt                              >  2:310510‑M2/33‑76 (MQ=255)                                                    cccgacacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTACACCTGCGcc                   <  1:569089‑M2/55‑1 (MQ=255)                                                     ccgacacgctccagtagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCa                  <  2:749299‑M2/56‑1 (MQ=255)                                                        acacgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAggg               >  2:509742‑M2/18‑76 (MQ=255)                                                           cgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGTTTTATCGCGTTTAAGCACTTTCACCTGCTCCAGGGTCa            <  1:806370‑M2/62‑1 (MQ=255)                                                           cgctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCa            <  1:80832‑M2/62‑1 (MQ=255)                                                            gctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAg           >  1:379162‑M2/14‑76 (MQ=255)                                                            gctccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAg           >  2:56497‑M2/14‑76 (MQ=255)                                                              tccagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGt         >  1:1045741‑M2/12‑76 (MQ=255)                                                                cagcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCCTTTAANCACTTTNACCTGCGCCAGGGTCATgttn       >  1:972415‑M2/10‑75 (MQ=255)                                                                 agcagtttCAGGGTTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGtttt      <  2:367744‑M2/68‑1 (MQ=255)                                                                  gcagtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTc     <  2:784371‑M2/69‑1 (MQ=255)                                                                     gtttCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAt  <  1:330008‑M2/72‑1 (MQ=255)                                                                         |                                                                        TCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAT  >  REL606/666062‑666205

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