New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	2258188 =	33 (0.440)	16 (0.230)	15/254	4.8	24.8%	intergenic (+200/‑192)	PP_1990/asd	transposase/aspartate‑semialdehyde dehydrogenase
	?	NC_002947	= 4022407	66 (0.930)					intergenic (‑58/+610)	PP_3547/emrB	NAD(P)‑binding oxidoreductase/multidrug transporter permease
Rejected: Coverage evenness skew score above cutoff.

GCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGA  >  NC_002947/2258116‑2258336                                                                         |                                                                                                                                                     gCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaaca                                                                         <  1:1034933/150‑1 (MQ=31)   aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagc                                                                       >  2:1150652/1‑150 (MQ=32)    ggACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagca                                                                      <  1:657613/150‑1 (MQ=33)     gcccGGGCCGAAGGCCCCATCCTGCGGCCTCAGCGCTTGACGCCAACCCCTGCCCCTCACACGCCTCCCTAACGCCTGCGTTCCGCCTCGTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaa                                                                     <  2:211669/148‑1 (MQ=11)     gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaa                                                                     >  1:916342/1‑150 (MQ=34)       ccGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaaca                                                                   <  1:331177/150‑1 (MQ=35)       ccGCGCCGAAGGCATCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATTAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaaca                                                                   <  1:1785245/150‑1 (MQ=255)         gcgcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaacagc                                                                 <  1:893514/150‑1 (MQ=38)            ccGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaacagcaac                                                              <  2:1276152/150‑1 (MQ=39)             cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGGCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCTAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaacagcaaca                                                             <  2:312375/150‑1 (MQ=255)                     ccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGc                                                     >  2:348825/1‑150 (MQ=255)                     ccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGc                                                     >  1:1321497/1‑150 (MQ=255)                     ccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATAAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGc                                                     <  2:284576/150‑1 (MQ=255)                          ttGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGAGCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAg                                                <  1:1619524/149‑1 (MQ=255)                            ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCa                                              >  2:15807/1‑150 (MQ=255)                                    aGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCaa                                      <  2:383120/150‑1 (MQ=255)                                       gcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAcagcaacagcaacagcaa                                                         >  2:324004/1‑128 (MQ=255)                                       gcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAcagcaacagcaacagcaa                                                         <  1:324004/128‑1 (MQ=255)                                          tGACGGGCATGGATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTgg                                >  1:1581247/1‑150 (MQ=255)                                          tGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTgg                                <  1:15807/150‑1 (MQ=255)                                          tGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCAGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTgg                                <  2:1156973/150‑1 (MQ=255)                                           gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTggg                               <  2:1321497/150‑1 (MQ=255)                                            aCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTgggg                              <  2:1257591/150‑1 (MQ=255)                                                     cATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATa                                                       >  1:1090736/1‑116 (MQ=255)                                                     cATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATa                                                       <  2:1090736/116‑1 (MQ=255)                                                      aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGCAAGTAAAGGACAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGa                    <  1:691216/150‑1 (MQ=21)                                                           cGTCACCTGCCTCCGCAAGGCCTGCGTTCGGACTCCTGAAGTAGAAAAAGACAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACaaaa               >  1:793200/1‑150 (MQ=25)                                                                 cTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTAc         >  1:283838/1‑150 (MQ=255)                                                                    ccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAgg                          >  1:1765218/1‑130 (MQ=255)                                                                    ccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAgg                          <  2:1765218/130‑1 (MQ=255)                                                                        cGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGCTCAAGATCTTCATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGa  <  1:1725506/150‑1 (MQ=255)                                                                        cGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaa                                                                     <  2:1116279/83‑1 (MQ=38)                                                                        cGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATcaacagcaa                                                                     >  1:1116279/1‑83 (MQ=38)                                                                         |                                                                                                                                                     GCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAACAGCAACAGCAACAGCAATAGCTACAGCAATAGGCAACGCTGGGGGAGGCACAGACAAAACCCTACAATCCGA  >  NC_002947/2258116‑2258336

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.