New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979271 =	1 (0.020)	39 (0.680)	39/488	0.4	97.5%	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

TCCAACGCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGG  >  NC_000913/1979021‑1979506                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            cccAACGCCATTCAGGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCAACATACACTAAATCAGTAAGTTGGCAGCATCAccc                                                                                                                                                                                                                                             <  2:53293/250‑1 (MQ=11)           cggataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGtat                                                                                                                                                                                                                                   >  1:447981‑M1/241‑251 (MQ=255)            ggataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGtata                                                                                                                                                                                                                                  >  1:294289‑M1/240‑251 (MQ=255)              ataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggaacagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacctgatctcatatcagggacttgttcgcaccttccttaaCATAAATGtata                                                                                                                                                                                                                                  <  2:180160‑M1/12‑1 (MQ=255)                aaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccctggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgacttcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATCAATGTATAAGTc                                                                                                                                                                                                                              <  2:231410‑M1/16‑1 (MQ=255)                aagggcgccggccattaccagccttggggtaaaacggctcgatgacttcctccatgtttttccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcgtcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagtgacttgttctccccttccttacCATAAATGTATAAGt                                                                                                                                                                                                                               >  1:174059‑M1/236‑250 (MQ=255)                     cgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaaaatgatctcatctcagggtcttgttcgcaccctccttaaCATAAATGTATAAGTAATACt                                                                                                                                                                                                                         <  2:281824‑M1/21‑1 (MQ=255)                         ggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACttttg                                                                                                                                                                                                                     >  1:333813‑M1/227‑251 (MQ=255)                                        tggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgttctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                                                                                                        >  1:506080‑M1/212‑249 (MQ=255)                                           ggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                                                                                                   <  2:294671‑M1/43‑1 (MQ=255)                                              aaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                                                                                                                <  1:183944‑M1/46‑1 (MQ=255)                                                 acggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                                                                                                                                                             >  1:231410‑M1/203‑251 (MQ=255)                                                   ggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattaactgccggcgaaggtcagttgatgactcatgatgaaccctgttctatggatccagatgacaaacatgatctcctatcagggccttgttcgcaccttccttaaCATAACTGTATAAGGCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                                                                                                                                                            <  1:356462‑M1/50‑1 (MQ=255)                                                            acttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                                                                                                                  >  1:284073‑M1/192‑251 (MQ=255)                                                                            gccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcaaatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                                                                                                                                                   >  2:454846‑M1/176‑250 (MQ=255)                                                                             ccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                                                                                                                 >  2:261399‑M1/175‑251 (MQ=255)                                                                                      aatctgctccatgcgggacaagaaaatatcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAAc                                                                                                                                                        <  2:57825‑M1/86‑1 (MQ=255)                                                                                            ctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                                                                                                                  <  1:282683‑M1/92‑1 (MQ=255)                                                                                               cattcgggacaagaaaatcgcttttctggtctgacggcgcttcctgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggccccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTCGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                                                                                                               >  2:418975‑M1/157‑251 (MQ=255)                                                                                                atgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                                                                                                              <  1:5438‑M1/96‑1 (MQ=255)                                                                                                  gcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggcgccagatgacaaacatgatctcatatcagggacttggtcgcacctttcttaaCATAAATGTATAAGACATACTTGTGTGGTGGGGGTAGTGCAAGCTGTGAAAAAGTTTTGCGCTTCGGGAGCAAGCAAAAATGAAACAGGAATAATCg                                                                                                                                             >  2:498159‑M1/154‑250 (MQ=255)                                                                                                     ggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                                                                                                         <  2:423802‑M1/101‑1 (MQ=255)                                                                                                                   cttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcac                                                                                                                           <  1:356276‑M1/115‑1 (MQ=255)                                                                                                                        ctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                                                                                                                       <  2:329549‑M1/119‑1 (MQ=255)                                                                                                                               gacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGGTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                                                                                                                  >  1:325258‑M1/125‑248 (MQ=255)                                                                                                                                 cggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGGTGCGCACAGTCAAAATAACTc                                                                                                             >  1:423802‑M1/123‑251 (MQ=255)                                                                                                                                     gcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                                                                                                                                                  <  2:411300‑M1/76‑1 (MQ=255)                                                                                                                                     gcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                                                                                                                                                  >  1:411300‑M1/119‑194 (MQ=255)                                                                                                                                            tctgaattcactgtcggcgaaggtaaagtgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccctccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCGCCGTAaata                                                                                                   >  1:73075‑M1/112‑250 (MQ=255)                                                                                                                                                aattcactgtcgtcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat                                                                                              >  2:409986‑M1/108‑251 (MQ=255)                                                                                                                                                aattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat                                                                                              <  2:292598‑M1/144‑1 (MQ=255)                                                                                                                                                          cggcgaaggtaagttgaagactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataa                                                                                     >  2:183944‑M1/98‑250 (MQ=255)                                                                                                                                                                    aagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACttt                                                                          >  1:245873‑M1/88‑251 (MQ=255)                                                                                                                                                                                         ccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAAATAAGCTa                                                     >  1:398195‑M1/67‑251 (MQ=255)                                                                                                                                                                                                    ggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATTCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTTAAAATAACTCACCGTAAATAATCAGCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACa                                           >  1:340402‑M1/56‑250 (MQ=255)                                                                                                                                                                                                     gctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGa                                         >  1:294671‑M1/55‑251 (MQ=255)                                                                                                                                                                                                                   acatgatctcatatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGcaca                            <  2:66758‑M1/210‑1 (MQ=255)                                                                                                                                                                                                                            catatcagggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTCCCGCACAACAAACt                     >  2:166638‑M1/32‑248 (MQ=255)                                                                                                                                                                                                                                  agggacttgttcgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATa            >  2:5438‑M1/26‑251 (MQ=255)                                                                                                                                                                                                                                             cgcaccttccttaaCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAgg  >  1:3408‑M1/15‑250 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            TCCAACGCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGG  >  NC_000913/1979021‑1979506

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