New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4166372	6 (0.170)	6 (0.170)	5/230	2.6	60.8%	intergenic (+87/‑287)	murI/rrsB	glutamate racemase/16S ribosomal RNA of rrnB operon
	?	NC_000913	= 4166399	2 (0.060)					intergenic (+114/‑260)	murI/rrsB	glutamate racemase/16S ribosomal RNA of rrnB operon

CGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACT  >  NC_000913/4166231‑4166515                                                                                                                                              |                                                                                                                                                cGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac                                                                                                                                          <  1:16103/149‑1 (MQ=255)      tACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg                                                                                                                                      <  5:150326/148‑1 (MQ=255)                 ggCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACacggcggcgtgtttgccgttgtt                                                                                                                          <  2:145737‑M1/149‑24 (MQ=255)                           cGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagt                                                                                                                <  1:99753‑M1/149‑34 (MQ=255)                                gAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                                                                                                                                             <  7:112625/115‑1 (MQ=255)                                gAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                                                                                                                                             >  8:112625/1‑115 (MQ=255)                                       tGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtgg                                                                                                              <  5:76145‑M1/139‑36 (MQ=255)                                       tGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtgg                                                                                                              >  6:76145‑M1/1‑104 (MQ=255)                                                   gCTGATTTGGTAGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCgg                                                                                          >  7:93475/1‑147 (MQ=255)                                                      gATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtgg                                                                                                              <  3:36397‑M1/124‑36 (MQ=255)                                                      gATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACccggcggcgtgtttgccgttgttccgtgtcagtgg                                                                                                              >  4:36397‑M1/1‑89 (MQ=255)                                                       aTTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCagag                                                                                     >  3:132290/1‑148 (MQ=40)                                                                              aGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGgtcag                                                                                                                 >  6:31390‑M2/1‑92 (MQ=255)                                                                                                                                       ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgcca              <  5:76145‑M2/139‑105 (MQ=255)                                                                                                                                       ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgcca              >  6:76145‑M2/1‑35 (MQ=255)                                                                                                                                       ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatc                             <  3:36397‑M2/124‑90 (MQ=255)                                                                                                                                       ccacTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatc                             >  4:36397‑M2/1‑35 (MQ=255)                                                                                                                                         acTGACACGGAACAACGGCAAACACGCCGCCGGgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttctgaccgcgcaacattcaaccaaatcagcctaaaactgccagtttttcgagcg  <  1:99753‑M2/149‑117 (MQ=255)                                                                                                                                          cTGACccggcggcgtgtttgccgttgttccgtgtcagtggtggcgcattatagggagttattccggcctgacaagaggaaatttaaaataattttct                                                     >  6:31390‑M1/1‑5 (MQ=255)                                                                                                                                              |                                                                                                                                                CGTTTTACAGCGTTACGGCTTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACT  >  NC_000913/4166231‑4166515

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