New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3941521	25 (0.180)	10 (0.080)	7/168	NT	29.8%	intergenic (‑194/‑287)	yieP/rrsC	putative transcriptional regulator/16S ribosomal RNA of rrnC operon
	?	NC_000913	= 3941548	NA (NA)					intergenic (‑221/‑260)	yieP/rrsC	putative transcriptional regulator/16S ribosomal RNA of rrnC operon

TGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGA  >  NC_000913/3941518‑3941647                               |                                                                                                    tGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGaaa                                >  2:1275351/1‑100 (MQ=11)  gACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAg                               >  1:405587/1‑100 (MQ=31)  gACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAg                               >  2:2531563/1‑100 (MQ=31)  gACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAg                               <  1:1956857/100‑1 (MQ=31)  gACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAg                               >  2:466791/1‑100 (MQ=31)   acacGGAACAACGGCAAACACGCCGCCGGGACAGCGGGGTTCTCCTGAGAACTCCGAAAAAAAAAAAAAAAATAAATGCTTGACTCTGTAGCGGGAAAg                               <  1:461794/99‑1 (MQ=11)    cacGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCg                             <  1:1058435/100‑1 (MQ=32)      cGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGta                           <  2:3089713/100‑1 (MQ=33)      cGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGta                           >  2:1299061/1‑100 (MQ=33)      cGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGta                           <  1:1350550/100‑1 (MQ=33)      cGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGta                           >  1:104175/1‑100 (MQ=32)       ggAACAACGTCAAACACGCCGCCAGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtat                          <  1:1649140/100‑1 (MQ=12)       ggAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtat                          <  2:3002380/100‑1 (MQ=33)       ggAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtat                          >  2:971630/1‑100 (MQ=33)       ggAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtat                          <  2:1285790/100‑1 (MQ=33)       ggAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGta                           <  1:20791/99‑1 (MQ=33)         aacaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCATAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatta                        <  2:2226718/100‑1 (MQ=18)         aacaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtattn                        <  2:3046807/100‑2 (MQ=18)         aacaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatta                        <  2:585446/100‑1 (MQ=33)         aacaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatta                        >  2:2320385/1‑100 (MQ=33)         aacaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatta                        <  2:978078/100‑1 (MQ=33)          acaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtattat                       <  2:3593547/100‑1 (MQ=33)          acaacGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAACCGtattat                       >  1:1660304/1‑100 (MQ=17)             acGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGca                    >  2:3709112/1‑100 (MQ=33)              cGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcac                   >  1:2527315/1‑100 (MQ=33)              cGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcac                   >  2:2868799/1‑100 (MQ=33)              cGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcac                   >  2:339605/1‑100 (MQ=33)              cGGCAAACACGACGCCGGGTCAGCGGGGTTATCTTGAGAACTACGGCAGAGAAAGCAAAAATAAATGCTTGACTTTGTAGCGGGAAAGCGTATTATGcac                   <  1:3724941/100‑1 (MQ=11)               ggCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcaca                  <  2:2754763/100‑1 (MQ=33)                gCAAACAGGCCGCCTGGTAAGCGCAGTTTTAATGAGAACTCCGCCAGAGAAAGTCGAAATAAATGCTTGACTTTCTAGCGGGAAAGCGTAGTATGCtcac                 >  1:177386/1‑100 (MQ=2)                gCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCacac                 >  1:2000285/1‑100 (MQ=33)                 cAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAcc                >  1:29237/1‑100 (MQ=33)                 cAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAcc                >  2:824276/1‑100 (MQ=33)                        gccgtcGGGCCACCGGCGGTCTCCTGAGAACTCCGGGATGGACAGCAAGAACAAACGCTTGACTCGGCAGCCGACAAGCGAATGATTCACACCCCGcgcc         >  1:1258316/1‑100 (MQ=2)                         ccgccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgccg        >  1:1473322/1‑100 (MQ=32)                         ccgccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGCGAAAGCGTATTATGCACACGCCgcgc          >  1:205529/1‑98 (MQ=12)                          cgccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgccgc       >  2:2873620/1‑100 (MQ=33)                          cgccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgccg        >  1:2144283/1‑99 (MQ=32)                          agccgGGTCAGCGGGGTTCTCCTGAGAACTGCGGCAGAGAAAGCAAAAATAACTGCTTGACTCTGTAGGAGGAAAGCGTATTATGCACACCCCGagcc         >  1:782595/2‑98 (MQ=11)                           gccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCt      <  1:1787918/100‑1 (MQ=32)                           gccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCt      >  2:3412851/1‑100 (MQ=33)                           gccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCt      >  1:3408281/1‑100 (MQ=33)                            ccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTg     >  1:2422810/1‑100 (MQ=33)                            ccgGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTg     <  2:2125030/100‑1 (MQ=33)                               ggTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTgaga  <  1:569590/100‑1 (MQ=33)                               |                                                                                                    TGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGA  >  NC_000913/3941518‑3941647

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