New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1826656	13 (0.700)	3 (0.180)	3/120	1.6	13.1%	coding (268/969 nt)	astE	succinylglutamate desuccinylase
	?	NC_000913	= 1826936	28 (1.660)					coding (1324/1344 nt)	astB	succinylarginine dihydrolase
Rejected: Frequency below/above cutoff threshold.

GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGC  >  NC_000913/1826593‑1826721                                                                |                                                                  gCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCagga                                                               <  1:509817/68‑1 (MQ=255)             cAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCgagggagggggctatg                                                   >  1:529811‑M1/1‑52 (MQ=255)              aaaCATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                  <  1:122750/68‑1 (MQ=255)              aaaCATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGa                                                  <  2:82234/68‑1 (MQ=255)                aCATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATc                                                <  1:332296/68‑1 (MQ=255)                              tCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCa                                  >  1:465405/1‑68 (MQ=255)                                 cTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACg                               <  2:625787/68‑1 (MQ=255)                                       cAATAACGTTTCCCTTGCTTCAGCGCgagggagggggcaatggataattttcttgctctgaccttaac                         >  2:406696‑M1/1‑26 (MQ=255)                                       cAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCGCCAGCAACCGCCAACGTAACgg                         >  1:204145/1‑68 (MQ=255)                                           aaCGTTTCCCTTGCTTCAGCGCgagggagggggcaatggataattttcttgctctgaccttaacgggt                     <  1:406698‑M1/68‑47 (MQ=255)                                             cGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGAtctc                   >  1:310843/1‑68 (MQ=255)                                              gTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGAATCACAAGCAACAGCCAACGTAACGGGAACTCg                  >  2:63721/1‑68 (MQ=255)                                                ttCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGcc                <  1:433949/68‑1 (MQ=255)                                                  cccTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGt              <  2:576784/68‑1 (MQ=255)                                                       gCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGat         <  2:465403/68‑1 (MQ=255)                                                              cgcAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCgcgc  >  1:406823/1‑68 (MQ=255)                                                                |                                                                  GCCAACGACCGCCAAACATTCGATTCATATCGCTATGGCAATAACGTTTCCCTTGCTTCAGCGCAGGAGGATTCCCGAGGATCACCAGCAACCGCCAACGTAACGGGATCTCGCCGTGAGATATCGCGC  >  NC_000913/1826593‑1826721

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