New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 223839	22 (0.330)	129 (2.420) +ATGAAGGGAGCTTGCTCCTGGATTCAGC	17/242	4.3	86.4%	noncoding (69/1542 nt)	rrsH	16S ribosomal RNA of rrnH operon
	?	NC_000913	3941911 =	28 (0.430)					noncoding (104/1542 nt)	rrsC	16S ribosomal RNA of rrnC operon
Rejected: Coverage evenness skew score above cutoff.

GCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTA  >  NC_000913/223689‑223889                                                                                                                                                       |                                                   gCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACgg                                                    >  1:27492/1‑151 (MQ=31)             aTGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGgaagaa                                        <  1:432712/151‑1 (MQ=31)             aTGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGgaaga                                         >  2:291760/1‑150 (MQ=31)              tGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGgaaga                                         >  1:61623/1‑149 (MQ=31)                 aTACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgctt                                    >  1:438654/1‑151 (MQ=31)                   aCCAAGCCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTACATTGAAGAGTTTGATCATGTCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttg                                   >  1:605635/1‑150 (MQ=9)                    ccAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttgct                                 >  1:231998/1‑151 (MQ=31)                    ccAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttgct                                 >  2:327238/1‑151 (MQ=31)                            cAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTg                           <  1:417946/149‑1 (MQ=31)                            cAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGc                          <  1:895323/150‑1 (MQ=31)                             aaGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTg                        <  2:904441/151‑1 (MQ=31)                                      cacGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc               >  2:415246/1‑151 (MQ=31)                                       acGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCg              >  2:725466/1‑151 (MQ=31)                                        cGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGCAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCgg             <  2:814470/151‑1 (MQ=14)                                          tAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGa            >  1:148963/1‑150 (MQ=31)                                          tAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAc           >  2:1018384/1‑151 (MQ=31)                                           aaTTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAc           <  2:319612/150‑1 (MQ=31)                                               cATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGt       <  1:971145/150‑1 (MQ=31)                                                aTTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTg      >  1:483094/1‑150 (MQ=31)                                                aTTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGa     <  2:137789/151‑1 (MQ=31)                                                   aCAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGt   <  2:945088/150‑1 (MQ=31)                                                    cAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTa  >  1:143158/1‑150 (MQ=31)                                                     aaaGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTC‑AACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAAGGGTGAGTa  >  1:906447/1‑148 (MQ=2)                                                                                                                                                       |                                                   GCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACAAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTA  >  NC_000913/223689‑223889

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