New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2482514	39 (1.510)	5 (0.230)	4/112	1.8	14.7%	coding (826/1164 nt)	emrK	multidrug resistance efflux pump membrane fusion protein
	?	NC_000913	= 2482573	25 (1.140)					coding (767/1164 nt)	emrK	multidrug resistance efflux pump membrane fusion protein
Rejected: Frequency below/above cutoff threshold.

AACACAACATTTTCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGCCATTAACG  >  NC_000913/2482507‑2482636                                                                   |                                                                aaCACAACATTTTCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAg                                                                <  1:138852‑M2/68‑2 (MQ=255)  acacAACATTTTCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGCGAgt                                                               >  2:173835‑M2/1‑66 (MQ=255)     caacaTTTTCACCATAAAGATCGCTGATATTATTGTCCGTTTGACCAATCCGTTCATCCGTGAgttgt                                                            >  2:946734‑M2/1‑63 (MQ=255)          ttttCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttc                                                       >  1:916150‑M2/1‑58 (MQ=255)          ttttCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttc                                                       >  2:958155‑M2/1‑58 (MQ=255)          ttttCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttc                                                       >  1:590956‑M2/1‑58 (MQ=255)              cACCATAAAGCTCGCAGAAAATATTGAACGATTGACCAATACGAACATCCGTGAGTTGTGTTTCTTTa                                                   >  2:954184/1‑68 (MQ=255)              cACCATAAAGATCGCTGATAATATTGCCCGATTGCCCAATCCGTACATCCGTGAGTTGTGTTTCTTTa                                                   >  2:601830/1‑68 (MQ=255)              cACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttcatgg                                                   <  1:13475‑M2/68‑15 (MQ=255)              cACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttcatgg                                                   <  1:498674‑M2/68‑15 (MQ=255)              cACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTa                                                   >  2:411446/1‑68 (MQ=255)               aCCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTag                                                  >  1:740934/1‑67 (MQ=255)                      aGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGc                                           <  2:262697/68‑1 (MQ=255)                       gATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttc                                                       >  1:218767‑M2/1‑45 (MQ=255)                       gATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttc                                                       <  2:218767‑M2/55‑11 (MQ=255)                                taataTTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCcac                                 <  1:869180/68‑1 (MQ=255)                                 aataTTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGGTTCTTTAAAGTTGGCATTAACCcaca                                <  2:220009/68‑1 (MQ=255)                                    atTGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttcatggtcgggtgacagggatcaatatg                             >  2:13475‑M2/1‑32 (MQ=255)                                     ttGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttcatggtcgggtgacagggatcaatatgg                            <  1:263406‑M2/68‑38 (MQ=255)                                     ttGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTg                            >  1:151160/1‑68 (MQ=255)                                     ttGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTg                            >  2:733200/1‑68 (MQ=255)                                     ttGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTg                            >  1:809856/1‑68 (MQ=255)                                     ttGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTg                            >  1:685158/1‑68 (MQ=255)                                      tGACCGATTGGCCAGTCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGa                           >  2:655384/1‑68 (MQ=255)                                      tGACCGATTGACCAATCCGTACATCCGTGAgttgtgtttcatggtcgggtgacagggatcaatatggg                           >  2:263406‑M2/1‑30 (MQ=255)                                        aCCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACg                         >  1:858504/1‑68 (MQ=255)                                         ccGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGCCGt                        >  2:222022/1‑68 (MQ=255)                                         ccGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGt                        >  2:516549/1‑68 (MQ=255)                                              agACCAATACGTGCGGCCGTGTATAGTGTTAATATAAAGTTGGCATTAAAACACCTAAGACGTGACgg                   >  2:166296/2‑68 (MQ=255)                                                 ccAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTAc                >  2:690749/1‑68 (MQ=255)                                                 ccAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTAc                >  1:198755/1‑68 (MQ=255)                                                      ccGTACATCCGGGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAg           <  2:212127/68‑1 (MQ=255)                                                       cgaaCATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGc          <  1:819544/65‑1 (MQ=255)                                                          aCATCCGTGAGTTGTGTTTCTTTAAAGTTTGCATTAACCCACATTTGACGTGCCGGTACGACAGCCAt       <  1:485473/68‑1 (MQ=255)                                                          aCATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGCCAt       <  2:587827/68‑1 (MQ=255)                                                          aCATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGCCAt       <  2:925512/68‑1 (MQ=255)                                                           cATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGCCAtt      >  2:866551/1‑68 (MQ=255)                                                               cGTGAGATGTGTATCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGCCATTAACg  >  1:676157/1‑68 (MQ=255)                                                                   |                                                                AACACAACATTTTCACCATAAAGATCGCTGATAATATTGACCGATTGACCAATCCGTACATCCGTGAGTTGTGTTTCTTTAAAGTTGGCATTAACCCACATTTGACGTGCCGGTACGACAGCCATTAACG  >  NC_000913/2482507‑2482636

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