New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 734441	57 (0.650)	4 (0.050)	4/266	NT	6.8%	pseudogene (222/1434 nt)	ybfO	conserved protein, rhs‑like YbfO
	?	NC_000913	= 3622200	NA (NA)					coding (3009/4236 nt)	rhsB	rhs element protein RhsB

GGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAG  >  NC_000913/3622152‑3622346                                                 |                                                                                                                                                   ggTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGAt                                               >  1:985834/1‑150 (MQ=32)  gTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATa                                              >  2:8706/1‑150 (MQ=32)   tCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAc                                             <  1:1535224/150‑1 (MQ=32)        aGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAAc                                        <  2:505149/150‑1 (MQ=32)          tCGCTATCTTTACGAACCGCTGGGCCGCAGGGTGGCAAAACGGGAGTGGCGGCGTGAACGGGACCTGACGGGCTGGATTTCGCTGTCACGGACACCGCAAGTGCCCTGGAACGGCTGGGACGGCGACCGGCTGACCACGAAACAGAACGa                                      >  2:1434430/1‑150 (MQ=11)           cGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGAc                                     <  1:630361/150‑1 (MQ=32)            gCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACa                                    <  1:1296362/150‑1 (MQ=32)            gCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACa                                    >  2:1483887/1‑150 (MQ=32)            gCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACa                                    >  1:319853/1‑150 (MQ=32)                 cTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGc                               <  2:369778/150‑1 (MQ=35)                  tttACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGcc                              <  2:479902/150‑1 (MQ=35)                  tttACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACAGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGcc                              >  2:692389/1‑150 (MQ=25)                    tACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGc                            <  2:837734/150‑1 (MQ=37)                      cGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCAt                          >  2:1347934/1‑150 (MQ=37)                       gACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATc                         >  1:368961/1‑150 (MQ=37)                         cccGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCa                       <  1:1220211/150‑1 (MQ=37)                          ccGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAg                      >  2:129205/1‑150 (MQ=37)                               gggCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGa                                 >  2:204415/1‑134 (MQ=33)                               gggCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGa                                 <  1:204415/134‑1 (MQ=33)                                ggCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGAtt                <  1:901687/150‑1 (MQ=37)                                 gCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGAttt               >  2:922645/1‑150 (MQ=37)                                 gCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGAttt               >  2:43415/1‑150 (MQ=37)                                   cGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACTGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTAt             <  1:248055/150‑1 (MQ=37)                                   cGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTAt             <  1:210845/150‑1 (MQ=37)                                    gCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATc            <  1:189100/150‑1 (MQ=37)                                      aGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAg          <  2:915993/150‑1 (MQ=37)                                       gggTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGc         >  1:500418/1‑150 (MQ=37)                                         gTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGCCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCg       >  1:749959/1‑150 (MQ=37)                                           ggCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCggg     <  2:983050/150‑1 (MQ=37)                                           ggCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCggg     >  2:915617/1‑150 (MQ=37)                                           ggCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCggg     <  1:1483887/150‑1 (MQ=37)                                           ggCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCggg     >  2:1299069/1‑150 (MQ=37)                                           cgCAAAACGGATGTGGCGGCGTGAACAGGACCTCACCGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCggg     >  2:1195682/2‑150 (MQ=255)                                             cAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGa   >  1:95868/1‑150 (MQ=37)                                              aaaaCGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAg  <  2:716300/150‑1 (MQ=37)                                                 |                                                                                                                                                   GGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAG  >  NC_000913/3622152‑3622346

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