New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2904081 =	41 (0.630)	5 (0.080)	5/246	NT	11.4%	intergenic (+707/+666)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis
	?	NC_000913	2904099 =	39 (0.640)					intergenic (+725/+648)	ygcE/queE	putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis

ATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAAC  >  NC_000913/2903951‑2904236                                                                                                                                                      |                                                                                                                                          gccagtgttccccgcgccagcggggataaaccgagcagatacacggctttgtattccgtgcgcccgtgttccccgcgccagcggggataaaccgaatagcaatagtccatagatttgcgaaaacaggtgtTCCCCGCTGGC‑GCGGGGaa                                                                                                                                           >  1:433962‑M1/131‑149 (MQ=255) atGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGaa                                                                                                                                           <  3:45683/149‑1 (MQ=255)   gTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAaca                                                                                                                                         >  4:348001/1‑149 (MQ=255)       cTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTg                                                                                                                                     >  4:68346/1‑149 (MQ=255)           ttaCAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGtttt                                                                                                                                 >  4:33949/1‑149 (MQ=255)               aaGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCCACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCa                                                                                                                             >  2:79340/1‑149 (MQ=255)                     tGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTa                                                                                                                       >  4:4055/1‑149 (MQ=255)                       cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCCACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATg                                                                                                                     <  1:79340/149‑1 (MQ=255)                              aaaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTAt                                                                                                                      >  3:124029/1‑141 (MQ=255)                              aaaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTAt                                                                                                                      <  4:124029/141‑1 (MQ=255)                               aaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTg                                                                                                             >  3:16513/1‑149 (MQ=255)                                       agaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCgg                                                                                                     <  3:348001/149‑1 (MQ=255)                                        gaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGt                                                                                                    >  3:284079/1‑149 (MQ=255)                                        gaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGt                                                                                                    >  3:237716/1‑149 (MQ=255)                                          aaGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGttt                                                                                                  <  1:593190/149‑1 (MQ=255)                                            gagaTAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTAt                                                                                                <  2:222155/149‑1 (MQ=255)                                                   cGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCt                                                                                         <  4:43202/149‑1 (MQ=255)                                                       ttATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCaaa                                                                                                                           >  4:281054/1‑111 (MQ=255)                                                       ttATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCaaa                                                                                                                           <  3:281054/111‑1 (MQ=255)                                                           atGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCggg                                                                                 <  4:237716/149‑1 (MQ=255)                                                                         ccGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCAc                                                                   >  1:382817/1‑149 (MQ=255)                                                                         ccGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCAc                                                                   >  3:178071/1‑149 (MQ=255)                                                                          ccgcgccagcggggataaaccgaatagcaatagtccatagatttgcgaaaacaggtgttccccgctggcgcggggaACACCt                                                                                                                                      >  1:189266‑M2/77‑82 (MQ=255)                                                                          ccgcgccagcggggataaaccgaatagcaatagtccatagatttgcgaaaacaggtgttccccgctggcgcggggaACACCt                                                                                                                                      <  2:189266‑M2/6‑1 (MQ=255)                                                                              ggCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAAt                                                              <  4:325870/149‑1 (MQ=255)                                                                                   ggggAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACaaa                                                         <  3:101031/149‑1 (MQ=255)                                                                                   ggggAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACaaa                                                         <  2:336111/149‑1 (MQ=255)                                                                                      gAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGcc                                                      <  4:178071/149‑1 (MQ=255)                                                                                        aCTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCgt                                                    <  3:114929/149‑1 (MQ=255)                                                                                        aCTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACGATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCgc                                                    <  2:9158/149‑2 (MQ=255)                                                                                            gACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTAt                                                <  3:33949/149‑1 (MQ=255)                                                                                             aCAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGcc                                                      <  4:534677/142‑1 (MQ=255)                                                                                             aCAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGcc                                                      >  3:534677/1‑142 (MQ=255)                                                                                                    ggCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCgg                                        >  3:195876/1‑149 (MQ=255)                                                                                                      ccTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGtt                                      >  3:570079/1‑149 (MQ=255)                                                                                                      ccTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGtt                                      <  4:284079/149‑1 (MQ=255)                                                                                                           gtagtCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGttt                                     >  3:349721/1‑145 (MQ=255)                                                                                                           gtagtCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGttt                                     <  4:349721/145‑1 (MQ=255)                                                                                                            tagtCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATccc                                >  1:227724/1‑149 (MQ=255)                                                                                                            tagtCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATccc                                >  2:23750/1‑149 (MQ=255)                                                                                                               tctcGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGc                             <  3:4055/149‑1 (MQ=255)                                                                                                                   gTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGgc                         >  2:4027/1‑149 (MQ=255)                                                                                                                      aGGCTCCTGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGcgcg                      <  3:274439/149‑1 (MQ=255)                                                                                                                          tCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGa                  >  1:398892/1‑149 (MQ=255)                                                                                                                             aggtgtTCCCCGC‑GCCAGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCgg                                                                                  <  2:189266‑M1/76‑1 (MQ=255)                                                                                                                             aggtgtTCCCCGC‑GCCAGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCgg                                                                                  >  1:189266‑M1/7‑82 (MQ=255)                                                                                                                                   ttccccgcgccagcggggaACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGc          >  1:433962‑M2/20‑149 (MQ=255)                                                                                                                                     ccccgcgccagcggggaACACCTGTTTTCGCAAATCTATGGa                                                                                                                   <  4:274111‑M2/25‑1 (MQ=255)                                                                                                                                     ccccgcgccagcggggaACACCTGTTTTCGCAAATCTATGGa                                                                                                                   >  3:274111‑M2/18‑42 (MQ=255)                                                                                                                                     CCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTAt                                                                                                                      <  1:123533‑M1/38‑1 (MQ=255)                                                                                                                                     CCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTAt                                                                                                                      >  2:123533‑M1/1‑38 (MQ=255)                                                                                                                                          tGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGTGCGGGGAACACTGGCTCTGCAAc  <  4:195876/149‑1 (MQ=255)                                                                                                                                                    gaACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAg                                                        >  3:394545‑M2/3‑86 (MQ=255)                                                                                                                                                    gaACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAg                                                        <  4:394545‑M2/84‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                          ATGTTGCTTATTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGC‑GCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGGGAACACTGGCTCTGCAAC  >  NC_000913/2903951‑2904236

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