New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4040480	NA (NA)	12 (0.170)	5/282	10.3	30.1%	intergenic (+58/‑36)	rrlA/rrfA	23S ribosomal RNA of rrnA operon/5S ribosomal RNA of rrnA operon
	?	NC_000913	= 4040504	29 (0.390)					intergenic (+82/‑12)	rrlA/rrfA	23S ribosomal RNA of rrnA operon/5S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.

TCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGC  >  NC_000913/4040462‑4040652                                           |                                                                                                                                                     tCAGCCTGATACAGATGAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATGTGCCTGGCGGCAGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGTCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCAt                                            <  1:592597/149‑1 (MQ=255)     ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGCAAAACCATAGCGCCAATGGTAGTGTGGGGTCTCCTCATGCga                                        >  1:1154656/1‑149 (MQ=11)     ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCga                                        <  2:1234890/149‑1 (MQ=32)     ccTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCCGGCGGCAGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACGCAGAAGTTAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCga                                        <  2:294485/149‑1 (MQ=255)       tGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCgaga                                      <  1:1135181/149‑1 (MQ=32)          tACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGCGTa                                   >  1:78114/1‑149 (MQ=17)          tACAGATTAAATCAGAACGCAGAAGCGGTCTGACAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTa                                   <  2:1165382/149‑1 (MQ=17)               aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAAGTCAGAAGTGAAAAGCCGTAGCGCCGATGGTAGTGGGGGGTCTCCTCAGGCCAGGCTAGGGaa                              >  2:856336/1‑149 (MQ=11)               aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGaa                              <  1:621184/149‑1 (MQ=32)               aTTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCAAGAGTAGGGaa                              >  1:283066/1‑149 (MQ=17)                tgcaagcAGAACGCAGAAGTGGTCTGATAAAACAGAATTGGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGGCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAAc                             <  1:1053357/143‑1 (MQ=11)                 tAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACt                            >  1:970338/1‑149 (MQ=32)                     tCAGAACGCAGAAGCCGTCTGATAAAACAGAATTTGCCCGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCa                        <  1:122556/149‑1 (MQ=17)                       aGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAgg                      <  2:853448/149‑1 (MQ=32)                       aGAACGCAGAAGCGGTCTGATAAAACAGAATCTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAgg                      <  2:1219487/149‑1 (MQ=17)                           cGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATc                  <  1:445425/149‑1 (MQ=32)                              aGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGACCTCAGAAGTGAAACGCCGTAGCGCCGATGGAAGTGTGGGGTCTCCTCATGAGAGAGTAGGGAACGTCCAGCCACCaaa               >  1:283621/1‑149 (MQ=11)                              aGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGGGAAGGGAACTGCCAGGCATCaaa               >  2:1027562/1‑149 (MQ=17)                              aGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCaaa               <  1:400335/149‑1 (MQ=32)                                  gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaa           <  1:1153554/149‑1 (MQ=32)                                   cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCCTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCAAATAATGTGAGAGTAGGGAACTGCCAGGCATCAAATaaaa          <  2:269975/149‑1 (MQ=11)                                   cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaaa          <  1:237136/149‑1 (MQ=32)                                    ggTCTGATAAAGCAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAAc         <  1:672716/149‑1 (MQ=17)                                     gTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGGGCAGGGAACTGCCAGGCCTGAAATAAAACg        >  2:1047643/1‑149 (MQ=255)                                      tCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGa       <  2:795013/149‑1 (MQ=32)                                         gATAAAACAGAATTTGCCCGGCGGCAGCAGCGCGGTGGACCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg    <  2:1035503/149‑1 (MQ=12)                                          aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAgg   >  2:233771/1‑149 (MQ=32)                                           tAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGc  >  1:421332/1‑149 (MQ=32)                                           tAAAACAGAACTCGCCTGGCGGCAGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGc  <  1:17328/149‑1 (MQ=12)                                           |                                                                                                                                                     TCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGC  >  NC_000913/4040462‑4040652

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