New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3653230	NA (NA)	50 (0.800)	36/276	0.6	98.0%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 3766089	1 (0.020)					coding (3907/4134 nt)	rhsA	rhs element protein RhsA

TGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCG  >  NC_000913/3653105‑3653368                                                                                                                              |                                                                                                                                           tGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAt                                                                                                                              >  1:248752/1‑140 (MQ=35)   gTCTGACGGCGCTTACTGATGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttt                                                                                                                            >  1:432321/1‑140 (MQ=37)    tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtttt                                                                                                                           <  1:457447/140‑1 (MQ=38)      tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtttttt                                                                                                                         <  2:221441/140‑1 (MQ=40)       gACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttttttt                                                                                                                        <  2:803842/140‑1 (MQ=255)           gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACAt                                                                                                                    <  1:884815/140‑1 (MQ=255)            cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATa                                                                                                                   <  2:1017934/140‑1 (MQ=255)              cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaaa                                                                                                                 >  1:935770/1‑140 (MQ=255)                 aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGaa                                                                                                              <  2:781098/140‑1 (MQ=255)                 aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGaa                                                                                                              <  2:515683/140‑1 (MQ=255)                 aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGaa                                                                                                              <  2:1044085/140‑1 (MQ=255)                  ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAac                                                                                                             >  2:381560/1‑140 (MQ=255)                     ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAcc                                                                                                          <  2:935770/140‑1 (MQ=255)                      tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAccc                                                                                                         <  2:432321/140‑1 (MQ=255)                          ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcc                                                                                                     >  2:532943/1‑140 (MQ=255)                          ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcc                                                                                                     >  2:980929/1‑140 (MQ=255)                           tCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccc                                                                                                    <  1:532943/140‑1 (MQ=255)                                gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                               <  2:410916/140‑1 (MQ=255)                                gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                               <  1:967948/140‑1 (MQ=255)                                  cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAAtt                                                                                             >  2:56937/1‑140 (MQ=255)                                                tgatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGAc                                                                               <  2:241517/140‑1 (MQ=255)                                                 gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACa                                                                              <  1:980929/140‑1 (MQ=255)                                                          tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGa                                                                     <  1:153880/140‑1 (MQ=255)                                                                 cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGacaaacaa                                                              >  1:425087/1‑140 (MQ=255)                                                                   cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAAt                                                            >  2:320908/1‑140 (MQ=255)                                                                     gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCg                                                          >  2:150943/1‑140 (MQ=255)                                                                      ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGa                                                         >  2:773283/1‑140 (MQ=255)                                                                           tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGtt                                                    >  1:50236/1‑140 (MQ=255)                                                                           tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGtt                                                    >  2:522847/1‑140 (MQ=255)                                                                            ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGttt                                                   <  1:381560/140‑1 (MQ=255)                                                                                ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                               >  2:963901/1‑140 (MQ=255)                                                                                ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                               >  1:36460/1‑140 (MQ=255)                                                                                ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                               >  1:67514/1‑140 (MQ=255)                                                                                ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                               >  2:701593/1‑140 (MQ=255)                                                                                 cAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAtttt                                              <  2:248752/140‑1 (MQ=255)                                                                                   gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCa                                            >  1:642447/1‑140 (MQ=255)                                                                                      gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACa                                         <  2:36460/140‑1 (MQ=255)                                                                                             aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTaa                                  <  2:642447/140‑1 (MQ=255)                                                                                             aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTaa                                  <  2:425087/140‑1 (MQ=255)                                                                                                aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCa                               >  2:13719/1‑140 (MQ=255)                                                                                                aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCa                               <  2:67514/140‑1 (MQ=255)                                                                                                aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCa                               >  1:8979/1‑140 (MQ=255)                                                                                                 tctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                              <  1:963901/140‑1 (MQ=255)                                                                                                 tctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                              <  1:56937/140‑1 (MQ=255)                                                                                                  ctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAtt                             <  2:50236/140‑1 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAtt                             >  1:818358/1‑130 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAtt                             <  2:818358/130‑1 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGAc                       <  2:358068/136‑1 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGAc                       >  1:358068/1‑136 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                   >  2:388442/1‑140 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                   >  1:147856/1‑140 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                   >  1:782383/1‑140 (MQ=255)                                                                                                            ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                   >  1:205351/1‑140 (MQ=255)                                                                                                                ttGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCtttt               >  1:949295/1‑140 (MQ=255)                                                                                                                     cGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAg          <  1:320908/140‑1 (MQ=255)                                                                                                                            tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGc   >  2:477621/1‑140 (MQ=255)                                                                                                                            tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGc   <  2:8979/140‑1 (MQ=255)                                                                                                                            tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGc   >  2:8189/1‑140 (MQ=255)                                                                                                                            tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGc   <  1:701593/140‑1 (MQ=255)                                                                                                                             ccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCg  >  1:842417/1‑140 (MQ=255)                                                                                                                              |                                                                                                                                           TGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCG  >  NC_000913/3653105‑3653368

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