New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 458317057 (0.610)8 (0.090) 8/188 NT 11.8% intergenic (‑122/+79) hsdM/hsdR DNA methyltransferase M/endonuclease R Type I restriction enzyme
?NC_000913 = 4583190 66 (0.750)intergenic (‑142/+59) hsdM/hsdR DNA methyltransferase M/endonuclease R Type I restriction enzyme

AGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTC  >  NC_000913/4583100‑4583336
                                                                                          |                                                                                                                                                  
agaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGAccg                                                                                                                                                 >  1:403533‑M2/1‑91 (MQ=255)
agaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGAccg                                                                                                                                                 <  2:403533‑M2/94‑4 (MQ=255)
 gaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACgtg                                                                                          >  2:979562/1‑148 (MQ=255)
  aaaTTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGAccgg                                                                                                                                                <  4:236209‑M2/93‑5 (MQ=255)
  aaaTTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGAccgg                                                                                                                                                >  3:236209‑M2/1‑89 (MQ=255)
            gggTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCt                                                                               >  4:568113/1‑148 (MQ=255)
                aaTATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGAccggg                                                                                                                                               >  4:422235‑M2/1‑75 (MQ=255)
                aaTATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGAccggg                                                                                                                                               <  3:422235‑M2/80‑6 (MQ=255)
                     acacAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCt                                                                               >  4:582504/1‑139 (MQ=255)
                     acacAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCt                                                                               <  3:582504/139‑1 (MQ=255)
                        caAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                             >  2:953555/1‑106 (MQ=255)
                        caAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                             <  1:953555/106‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAAc                                                                                                                                         <  4:322953/76‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAAc                                                                                                                                         >  3:322953/1‑76 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAg                                                                                                                                      >  2:281228/1‑79 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAg                                                                                                                                      <  1:281228/79‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAg                                                                                                                                      <  3:988447/79‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAg                                                                                                                                      >  4:988447/1‑79 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACgtg                                                                                          >  4:459060/1‑123 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACgtg                                                                                          <  3:459060/123‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAgg                                                                                     >  4:958996/1‑128 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAgg                                                                                     <  3:958996/128‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAgg                                                                                     <  1:688959/128‑1 (MQ=255)
                          aaaTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAgg                                                                                     >  2:688959/1‑128 (MQ=255)
                                 tCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGa                                                          >  1:108804/1‑148 (MQ=255)
                                    ggTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGAAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAAtt                                                       <  2:108804/148‑1 (MQ=255)
                                      ttGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaatga                                                                                                         >  4:752259/1‑93 (MQ=255)
                                      ttGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAATGa                                                                                                         <  3:752259/96‑1 (MQ=255)
                                       tGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaa                                                                                                            <  3:655057/92‑1 (MQ=255)
                                       tGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaa                                                                                                            >  4:655057/1‑92 (MQ=255)
                                       tGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                             >  4:1092219/1‑91 (MQ=255)
                                       tGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                             <  3:1092219/91‑1 (MQ=255)
                                         cATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaa                                                                                                            <  1:857247/90‑1 (MQ=255)
                                         cATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGaa                                                                                                            >  2:857247/1‑90 (MQ=255)
                                            cAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTg                                                           >  4:273624/1‑136 (MQ=255)
                                            cAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTg                                                           <  3:273624/136‑1 (MQ=255)
                                            cAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagc                                               >  2:653388/1‑148 (MQ=255)
                                              aGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagcag                                             >  3:561750/1‑148 (MQ=255)
                                               ggcggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACg                                                                                            >  2:329592/1‑100 (MQ=255)
                                               ggcggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACg                                                                                            <  1:329592/100‑1 (MQ=255)
                                                   gcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTAt                                        >  4:208303/1‑148 (MQ=255)
                                                    cAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATc                                       >  4:259125/1‑148 (MQ=255)
                                                        tgagtgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGa                                                                                                      <  2:943370/81‑1 (MQ=255)
                                                        tgagtgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGa                                                                                                      >  1:943370/1‑81 (MQ=255)
                                                         gagtgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGAt                                  >  2:406828/1‑148 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGtata                                                                                              <  4:1001189/85‑1 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGtata                                                                                              >  3:1001189/1‑85 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAAt                                                                >  4:215145/1‑115 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAAt                                                                <  3:215145/115‑1 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTAt                               <  1:979562/148‑1 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTAt                               >  4:899842/1‑148 (MQ=255)
                                                            tgaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTAt                               <  3:568113/148‑1 (MQ=255)
                                                             gaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                              >  3:1124686/1‑148 (MQ=255)
                                                             gaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                              >  4:330713/1‑148 (MQ=255)
                                                             gaATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATc                              <  3:259125/148‑1 (MQ=255)
                                                              aaTCCCCGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttg                                                      <  2:873755‑M2/123‑95 (MQ=255)
                                                              aaTCCCCGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttg                                                      >  1:873755‑M2/1‑29 (MQ=255)
                                                              aaTCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCg                             <  3:330713/148‑1 (MQ=255)
                                                               aTCCCCGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcgg                                                                   >  2:412075‑M2/1‑28 (MQ=255)
                                                               aTCCCCGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcgg                                                                   <  1:412075‑M2/109‑82 (MQ=255)
                                                                 ccccGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttct                                                                             >  4:326648‑M2/1‑26 (MQ=255)
                                                                 ccccGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttct                                                                             <  3:326648‑M2/97‑72 (MQ=255)
                                                                 ccccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTg                                                           >  3:774870/1‑115 (MQ=255)
                                                                 ccccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTg                                                           <  4:774870/115‑1 (MQ=255)
                                                                  cccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCaa                         <  1:406828/148‑1 (MQ=255)
                                                                   ccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGtata                                                                                              >  2:359580/1‑78 (MQ=255)
                                                                   ccGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGtata                                                                                              <  1:359580/78‑1 (MQ=255)
                                                                    cGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttgggccgctgcttaatttacggaactcacaatg                       >  2:870412‑M2/1‑23 (MQ=255)
                                                                    cGGGAGCGTACAGAAGTACGTGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttgggccgctgcttaatttacggaactcacaatg                       <  1:870412‑M2/148‑126 (MQ=255)
                                                                    cGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTc                                                                           <  4:779815/96‑1 (MQ=255)
                                                                    cGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTc                                                                           >  3:779815/1‑96 (MQ=255)
                                                                     gggAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGTAGCCTGa                                                                                                             <  3:800352/61‑1 (MQ=255)
                                                                     gggAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGa                                                                                                             >  4:800352/1‑61 (MQ=255)
                                                                      ggAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAg                                                                       <  4:342767/98‑1 (MQ=255)
                                                                      ggAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAg                                                                       >  3:342767/1‑98 (MQ=255)
                                                                       gAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGtt                    <  3:208303/148‑1 (MQ=255)
                                                                             aCAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGa                                                                      >  1:781133/1‑92 (MQ=255)
                                                                             aCAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGa                                                                      <  2:781133/92‑1 (MQ=255)
                                                                                  aGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCaa                         <  3:677825/132‑1 (MQ=255)
                                                                                  aGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCaa                         >  4:677825/1‑132 (MQ=255)
                                                                                  aGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCAt         <  4:561750/148‑1 (MQ=255)
                                                                                    tACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCg       >  2:312158/1‑148 (MQ=255)
                                                                                       gTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagcag                                             <  2:623677/107‑1 (MQ=255)
                                                                                       gTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagcag                                             >  1:623677/1‑107 (MQ=255)
                                                                                        tGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttgggccgctgcttaatttacggaactcaca                          <  2:455927‑M2/125‑123 (MQ=255)
                                                                                        tGAccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttgggccgctgcttaatttacggaactcaca                          >  1:455927‑M2/1‑3 (MQ=255)
                                                                                         gACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTc  <  4:1124686/148‑1 (MQ=255)
                                                                                          |                                                                                                                                                  
AGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTC  >  NC_000913/4583100‑4583336

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