New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3941549 =	58 (0.390)	15 (0.110)	3/168	NT	21.5%	intergenic (‑222/‑259)	yieP/rrsC	putative transcriptional regulator/16S ribosomal RNA of rrnC operon
	?	NC_000913	4035245 =	NA (NA)					intergenic (+92/‑286)	hemG/rrsA	protoporphyrin oxidase, flavoprotein/16S ribosomal RNA of rrnA operon

TGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACA  >  NC_000913/4035146‑4035262                                                                                                    |                  tGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGaca                   >  1:641330/1‑100 (MQ=255) tGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTACTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCAGTGaca                   >  2:397269/1‑100 (MQ=255)   aaaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg                 <  1:13513/100‑1 (MQ=255)   aaaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACg                 <  1:2947017/100‑1 (MQ=255)    aaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                <  1:1374130/100‑1 (MQ=255)    aaTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACgg                <  1:671666/100‑1 (MQ=255)     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa               <  1:1532360/100‑1 (MQ=255)     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa               >  2:3735664/1‑100 (MQ=255)     aTAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGa               <  2:2961880/100‑1 (MQ=255)      tAAGCATAAAGAATATAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaa              <  1:4118897/100‑1 (MQ=255)      ataGCAGAAAGAATACAAAATGCGCGGGCAGAAAACTATTTCAAGTTTCCTCTTGTCAGGCCGGAAGAACTCCCTATAATGCGCCACCACTACCACGGa               <  2:2572576/97‑1 (MQ=14)       aagaaaaaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac             <  2:3441103/94‑1 (MQ=255)       aaGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaac             >  1:2588685/1‑100 (MQ=255)        aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaaca            >  2:1641318/1‑100 (MQ=255)        aGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATCCGCCACCACTGACACGGaaca            >  2:1659101/1‑100 (MQ=255)          cATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGaacaac          <  1:1956036/100‑1 (MQ=40)           agaaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACTACTGACACGGAACAACg         <  2:809104/98‑1 (MQ=255)           aTAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACg         >  1:395546/1‑100 (MQ=40)             aaaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGc       >  1:51585/1‑100 (MQ=38)              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa      <  1:2175458/100‑1 (MQ=35)              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa      <  2:2705357/100‑1 (MQ=35)              aaGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCa      >  2:3302785/1‑100 (MQ=35)               aGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACGCGGAACAACGGCaa     >  2:3893002/1‑100 (MQ=25)                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac   <  2:1077742/100‑1 (MQ=32)                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac   <  2:3413538/100‑1 (MQ=32)                 aaTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAac   <  1:431831/100‑1 (MQ=32)                  aTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAaca  >  1:439864/1‑100 (MQ=31)                  aTAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAaca  >  1:3625248/1‑100 (MQ=31)                                                                                                    |                  TGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACA  >  NC_000913/4035146‑4035262

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