New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	771005 =	9 (0.480)	27 (1.450)	19/132	0.1	75.0%	intergenic (+394/‑453)	mngB/cydA	alpha‑mannosidase/cytochrome d terminal oxidase, subunit I
	?	NC_000913	= 1050545	NA (NA)					noncoding (768/768 nt)	IS1	repeat region

TCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAA  >  NC_000913/1050493‑1050602                                                     |                                                          tCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCt                                            <  1:412512/68‑1 (MQ=255)      gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTgaga                                       <  1:536211/68‑1 (MQ=255)                         aaCACTATCAATAATGTGGAGTCATTACCGACAATACGTACGCTTGTGATCAGGATTGTCTACTTTGa                    <  2:248118/68‑1 (MQ=37)                         aaCACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGa                    <  1:614212/68‑1 (MQ=255)                          acacTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAg                   >  2:247154/1‑68 (MQ=255)                           cacTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGt                  <  1:244231/68‑1 (MQ=255)                           cacTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGt                  >  2:68407/1‑68 (MQ=255)                             cTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGtat                <  2:407121/68‑1 (MQ=255)                               aTCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATaa              >  1:148830/1‑68 (MQ=255)                                  aaTAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACAt           >  2:407751/1‑68 (MQ=255)                                   aTAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATc          >  2:75229/1‑68 (MQ=255)                                    tAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATcc         >  2:500651/1‑68 (MQ=255)                                      aGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCtt       >  2:430509/1‑68 (MQ=255)                                      aGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCAACTTTGAGTATAACATCCtt       >  2:1630/1‑68 (MQ=255)                                         tGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTtgtg    <  1:525579/68‑1 (MQ=255)                                          tgAGGCATTACCTACAATACGTACGCTTGAGATCAGGATTGTCTTCTTTGTGTATAACATCCTTGTGa   <  1:207165/67‑1 (MQ=255)                                           gAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGaa  <  2:126607/68‑1 (MQ=255)                                                     |                                                          TCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAA  >  NC_000913/1050493‑1050602

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