New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979271 =	1 (0.020)	22 (0.490)	20/490	1.0	95.7%	intergenic (‑56/‑482)	insA5/uspC	IS1 family protein InsA/universal stress protein C
	?	NC_000913	= 2102947	NA (NA)					pseudogene (443/446 nt)	wbbL	interrupted rhamnosyltransferase WbbL

GCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAAC  >  NC_000913/1979027‑1979517                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       gCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCaccgac                                                                                                                                                                                                                                                  >  1:95559/1‑246 (MQ=11) ccagcggataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgctccttccttaaCATAAAt                                                                                                                                                                                                                                                  >  2:147755‑M1/245‑251 (MQ=255)                         taccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACttttgttttg                                                                                                                                                                                                                           >  1:333535‑M1/221‑250 (MQ=255)                          accagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTgg                                                                                                                                                                                                                          >  1:399196‑M1/220‑250 (MQ=255)                                            cggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                                                                                                                                                                        <  1:136623‑M1/49‑1 (MQ=255)                                                     gacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                                                                                                              <  1:64558‑M1/59‑1 (MQ=255)                                                                           ggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                                                                                                                        <  2:373708‑M1/81‑1 (MQ=255)                                                                                    tgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                                                                                                                                                                <  1:318007‑M1/89‑1 (MQ=255)                                                                                     gctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAGGAAACAg                                                                                                                                                                 >  2:309788‑M1/161‑248 (MQ=255)                                                                                        ccatgcgggacaagaaaatctcttgtctggtttgacggcgcttactgctgaattcactgtcggcgaaggtaagtttatgactcatgatgaatcctgttctatgtctccagatgacaaacatgatctcatattagggacttgttcgtaccttccttaaCATAAATGTATAATTCATACTGTTTTTTTGGGTGTATTTTAATCTGTTAAAAAGTTTTTGGCTACGCTGGCAAGCAAAaattaaacgtgaata                                                                                                                                                            >  2:307379‑M1/158‑243 (MQ=255)                                                                                                     gaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                                                                                                              >  2:41897‑M1/145‑251 (MQ=255)                                                                                                             cttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaatcctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcac                                                                                                                                      <  1:333146‑M1/115‑1 (MQ=255)                                                                                                                              cgctaaacggtacaaacactgtcggcgacggtcagctgatgactcatgatgaatcctgttatatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc                                                                                                                     <  2:276580‑M1/132‑1 (MQ=255)                                                                                                                                                                                    cctgttctatggctccagatgacaaacatgatcgcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTATTTCGCTACGCTAGCAAGCAAAAAGGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCCGCTATAAATAATCACTTTCATGCAATACCACATAAGCTa                                                                >  2:229554‑M1/66‑250 (MQ=255)                                                                                                                                                                                               gctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGa                                                    >  2:318007‑M1/55‑251 (MQ=255)                                                                                                                                                                                                          caaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACAAGATAAGCTATTTTTAAACAGACACTTACCGc                                          <  2:19296‑M1/207‑1 (MQ=255)                                                                                                                                                                                                             acatgatctcatatcagggacttgttcgcaccttccttaaCCTAACTGTATCAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCACGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCACAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGcaca                                       <  2:429555‑M1/210‑1 (MQ=255)                                                                                                                                                                                                              catgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCacaac                                     <  2:133947‑M1/212‑1 (MQ=255)                                                                                                                                                                                                                           cagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTaaa                         >  2:136623‑M1/27‑250 (MQ=255)                                                                                                                                                                                                                                   tgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaag                  >  1:171411‑M1/19‑249 (MQ=255)                                                                                                                                                                                                                                                cttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa                                                                                                                                                      <  2:259813‑M1/99‑1 (MQ=255)                                                                                                                                                                                                                                                cttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa                                                                                                                                                      >  1:259813‑M1/6‑104 (MQ=255)                                                                                                                                                                                                                                                  taaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAAc  >  2:7119‑M1/4‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       GCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAAC  >  NC_000913/1979027‑1979517

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