New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3423376 =	16 (0.170)	11 (0.120)	5/196	NT	48.2%	intergenic (+182/+47)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
	?	NC_000913	= 4213162	8 (0.090)					intergenic (+3/‑72)	rrfE/yjaA	5S ribosomal RNA of rrnE operon/stress‑induced protein

AGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTG  >  NC_000913/3423264‑3423413                                                                                                                 |                                      agtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt    >  3:47975‑M1/113‑148 (MQ=255)        cggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   <  4:1034728‑M1/5‑1 (MQ=255)        cggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   >  3:1034728‑M1/106‑110 (MQ=255)           tggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       >  1:343197‑M1/103‑103 (MQ=255)           tggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       <  2:343197‑M1/1‑1 (MQ=255)             gtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   <  1:265316‑M1/5‑1 (MQ=255)             gtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   >  2:265316‑M1/101‑105 (MQ=255)               cccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacttcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc       >  4:1038098‑M1/99‑131 (MQ=255)               cccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaCTAAAACAAAAGGCCCAGTCTTCCGACTGAGcc       <  3:1038098‑M1/33‑1 (MQ=255)                      gaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   >  1:119114‑M1/92‑96 (MQ=255)                      gaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   <  2:119114‑M1/5‑1 (MQ=255)                       cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                    >  3:661421/1‑110 (MQ=40)                       cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCAGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                    <  4:661421/110‑1 (MQ=255)                         ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                   >  2:635460/1‑93 (MQ=39)                         ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                   <  1:635460/93‑1 (MQ=39)                          ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       <  1:70762‑M1/1‑1 (MQ=255)                          ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       >  2:70762‑M1/88‑88 (MQ=255)                           catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGa          <  2:491071‑M1/30‑1 (MQ=255)                           catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGa          >  1:491071‑M1/87‑116 (MQ=255)                           catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       <  1:1103330‑M1/1‑1 (MQ=255)                           catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       >  2:1103330‑M1/87‑87 (MQ=255)                              cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg         <  4:71976/114‑1 (MQ=40)                              cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg         >  3:71976/1‑114 (MQ=40)                                  gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              <  4:713140/105‑1 (MQ=255)                                  gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              >  3:713140/1‑105 (MQ=40)                                        gaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTg           <  4:881736‑M1/29‑1 (MQ=255)                                        gaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTg           >  3:881736‑M1/74‑102 (MQ=255)                                               aacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       <  4:965349‑M1/1‑1 (MQ=255)                                               aacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaa                                       >  3:965349‑M1/67‑67 (MQ=255)                                                 cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   <  2:905020‑M1/5‑1 (MQ=255)                                                 cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATaaa                                   >  1:905020‑M1/65‑69 (MQ=255)                                                 cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt      >  1:358880‑M1/65‑98 (MQ=255)                                                 cgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt      <  2:358880‑M1/34‑1 (MQ=255)                                                      tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                        <  3:629843/75‑1 (MQ=255)                                                      tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGccc                        >  4:629843/1‑75 (MQ=255)                                                         aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      >  1:691611/1‑90 (MQ=255)                                                         aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      <  2:791787/90‑1 (MQ=255)                                                         aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      <  2:691611/90‑1 (MQ=255)                                                         aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt      >  1:791787/1‑90 (MQ=255)                                                                 cccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg  <  3:149945/86‑1 (MQ=255)                                                                 cccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg  >  4:149945/1‑86 (MQ=255)                                                                               cctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc       >  3:926911‑M1/35‑67 (MQ=255)                                                                               cctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc       <  4:926911‑M1/33‑1 (MQ=255)                                                                                       ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              >  3:474855/1‑52 (MQ=255)                                                                                       ggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa              <  4:474855/52‑1 (MQ=255)                                                                                                                 |                                      AGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTG  >  NC_000913/3423264‑3423413

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