New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1683535 =	20 (1.000)	4 (0.230)	3/114	1.7	24.1%	coding (654/1302 nt)	rstB	sensory histidine kinase of RstAB two‑component system
	?	NC_000913	1683597 =	8 (0.470)					coding (716/1302 nt)	rstB	sensory histidine kinase of RstAB two‑component system

AACTGCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGA  >  NC_000913/1683546‑1683663                                                    |                                                                   aaCTGCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGaa                                                    <  1:387889/68‑1 (MQ=255)      tcgataacgcaggcgcactaacggtgatcgcagttcgtgagagataACCTGCGCGCCGCCGAATCCCa                                               >  2:715290‑M2/47‑68 (MQ=255)         acgccGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGc                                            >  1:29691/4‑68 (MQ=255)          cacCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCg                                           >  1:346372/1‑68 (MQ=255)                agtcgcactaacggtgttcgcagttcgtgagcgataACCTGAGCGCCGCCGAATCCCAGGCGTTGAAt                                     <  1:715289‑M2/32‑1 (MQ=255)                   tGCGCCTGCGTTATCGACTGGAGATGAGCGATACCCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGt                                  >  2:186637/1‑68 (MQ=255)                   tGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGt                                  >  2:178877/1‑68 (MQ=255)                   tGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGt                                  >  2:49174/1‑68 (MQ=255)                    gcactaacggtgttcgcagttcgtgagcgataACCTGCGCGCCGCCGAATCCCAGGCGTTGAATCg                                   >  2:742928‑M2/33‑66 (MQ=255)                    gcactaacggtgttcgcagttcgtgagcgataACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCg                                   <  1:742927‑M2/34‑1 (MQ=255)                                tCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAAc                     <  2:449134/68‑1 (MQ=255)                                      ggAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAg               <  1:361967/68‑1 (MQ=255)                                             tgcgataTCCTGAGCGCCGCCGAATCCCAGTCGTTGAATCGTGATTTCAGTCAACTTGAAGCTTTTAt        <  2:500702‑M2/61‑1 (MQ=255)                                             agcgataACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCt             <  2:35123‑M2/56‑1 (MQ=255)                                             agcgataACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCt             >  1:35123‑M2/8‑63 (MQ=255)                                             agcgataACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAAt        >  2:404748‑M2/8‑68 (MQ=255)                                                gataACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGa     >  1:496477‑M2/5‑68 (MQ=255)                                                   aACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAaga  <  2:51660‑M2/67‑1 (MQ=255)                                                    |                                                                   AACTGCGAACACCGTTAGTGCGCCTGCGTTATCGACTGGAGATGAGCGATAACCTGAGCGCCGCCGAATCCCAGGCGTTGAATCGTGATATCAGTCAACTTGAAGCTTTAATTGAAGA  >  NC_000913/1683546‑1683663

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