New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1629411	37 (0.200)	4 (0.020)	4/226	NT	10.1%	coding (197/204 nt)	ydfZ	selenoprotein, function unknown
	?	NC_000913	= 1629426	36 (0.200)					intergenic (+8/+27)	ydfZ/ydfI	selenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase

TTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCAT  >  NC_000913/1629264‑1629559                                                                                                                                                    |                                                                                                                                                     ttATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATg                                                                                                                                                     >  1:787839/1‑149 (MQ=255)   aTGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGaa                                                                                                                                                   <  2:943134/149‑1 (MQ=255)       ttAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTaa                                                                                                                                               <  2:608657/149‑1 (MQ=255)         aGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGc                                                                                                                                             <  3:862309/149‑1 (MQ=255)           cGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgt                                                                                                                                           <  1:1644506/149‑1 (MQ=255)            ggCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtg                                                                                                                                          <  1:481535/149‑1 (MQ=255)              cACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGtgtg                                                                                                                                        <  3:1739692/149‑1 (MQ=255)                ccGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAa                                                                                                                                      >  4:178685‑M2/1‑148 (MQ=255)                ccGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAa                                                                                                                                      >  4:1053224‑M2/1‑148 (MQ=255)                 cGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAat                                                                                                                                     >  4:673677‑M2/1‑147 (MQ=255)                      acacTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATgccgc                                                                                                                                >  2:1038722/1‑149 (MQ=255)                         cTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGa                                                                                                                             <  2:431030/149‑1 (MQ=255)                            gCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATgg                                                                                                                          <  1:1820891/149‑1 (MQ=255)                                 aTCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGa                                                                                                                             >  3:358978/1‑141 (MQ=255)                                 aTCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGa                                                                                                                             <  4:358978/141‑1 (MQ=255)                                    cTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTg                                                                                                                  <  2:787839/149‑1 (MQ=255)                                         gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtttt                                                                                                             <  1:1038722/149‑1 (MQ=255)                                           ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtgt                                                                                                                                         >  3:285964/1‑119 (MQ=255)                                           ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtgt                                                                                                                                         <  4:285964/119‑1 (MQ=255)                                                aCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                        >  2:1469344/1‑131 (MQ=255)                                                aCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                        <  1:1469344/131‑1 (MQ=255)                                                  tGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGaa                                                                                                                                                   <  4:671422/102‑1 (MQ=255)                                                  tGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGaa                                                                                                                                                   >  3:671422/1‑102 (MQ=255)                                                   gAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACg                                                                                                   >  4:398252/1‑149 (MQ=255)                                                           gACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                   <  2:1378729/125‑1 (MQ=255)                                                           gACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                   >  1:1378729/1‑125 (MQ=255)                                                               gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAa                                                                                                                                      <  2:1735337‑M2/102‑2 (MQ=255)                                                               gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAa                                                                                                                                      >  1:1735337‑M2/1‑101 (MQ=255)                                                                   aGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                         >  2:1101991/1‑111 (MQ=255)                                                                   aGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                         <  1:1101991/111‑1 (MQ=255)                                                                       aaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATtcacacgcttagttcatgccgag                                                                                                                               >  4:794603‑M1/1‑78 (MQ=255)                                                                       aaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATtcacacgcttagttcatgccgag                                                                                                                               <  3:794603‑M1/101‑24 (MQ=255)                                                                        aTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                         >  1:55792/1‑106 (MQ=255)                                                                        aTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                         <  2:55792/106‑1 (MQ=255)                                                                             ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAat                                                                                                                                     >  4:522913‑M2/1‑87 (MQ=255)                                                                             ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAat                                                                                                                                     <  3:522913‑M2/89‑3 (MQ=255)                                                                                             gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaa                                                         <  3:398252/149‑1 (MQ=255)                                                                                                 ttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATtcacacgcttagttcatgccga                                                                                                                                <  2:1098489‑M1/74‑23 (MQ=255)                                                                                                 ttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATtcacacgcttagttcatgccga                                                                                                                                >  1:1098489‑M1/1‑52 (MQ=255)                                                                                                              agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAt                                                                                                                    >  1:1557467/1‑73 (MQ=255)                                                                                                              agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAt                                                                                                                    <  2:1557467/73‑1 (MQ=255)                                                                                                               gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTTACGGAATATTTTGCCa                                                                                      >  1:170529/1‑102 (MQ=255)                                                                                                               gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAtttt                                                                                          >  4:64275/1‑98 (MQ=255)                                                                                                               gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAtttt                                                                                          <  3:64275/98‑1 (MQ=255)                                                                                                               gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGccc                                                                                      <  2:170529/102‑2 (MQ=255)                                                                                                                      tGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                            <  3:498629/73‑1 (MQ=255)                                                                                                                      tGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                            >  4:498629/1‑73 (MQ=255)                                                                                                                                 gACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                     >  3:283968/1‑149 (MQ=255)                                                                                                                                    cTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                  >  4:687188/1‑149 (MQ=255)                                                                                                                                     tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAt                 <  3:687188/149‑1 (MQ=255)                                                                                                                                          cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGt                                                                                  >  3:271329/1‑79 (MQ=255)                                                                                                                                          cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGt                                                                                  <  4:271329/79‑1 (MQ=255)                                                                                                                                             ctcGGCATGAACTAAGCGTGTGAatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatt                                                         <  3:794603‑M2/101‑79 (MQ=255)                                                                                                                                             ctcGGCATGAACTAAGCGTGTGAatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatt                                                         >  4:794603‑M2/1‑23 (MQ=255)                                                                                                                                              tcGGCATGAACTAAGCGTGTGAatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaa                                                                                   <  2:1098489‑M2/74‑53 (MQ=255)                                                                                                                                              tcGGCATGAACTAAGCGTGTGAatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaa                                                                                   >  1:1098489‑M2/1‑22 (MQ=255)                                                                                                                                                   aTtcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatttgctcagcggtcagaccttcagtatcaatcgacaggatcttgccagtgtgaccg   >  4:673677‑M1/1‑2 (MQ=255)                                                                                                                                                   aTtcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcgg                                                               >  4:522913‑M1/1‑2 (MQ=255)                                                                                                                                                   aTtcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcgg                                                               <  3:522913‑M1/89‑88 (MQ=255)                                                                                                                                                    ttcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatttgctcagcggtcagaccttcagtatcaatcgacaggatcttgccagtgtgaccgg  >  4:178685‑M1/1‑1 (MQ=255)                                                                                                                                                    ttcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatttgctcagcggtcagaccttcagtatcaatcgacaggatcttgccagtgtgaccgg  >  4:1053224‑M1/1‑1 (MQ=255)                                                                                                                                                    ttcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatttgctcagc                                                 <  2:1735337‑M1/102‑102 (MQ=255)                                                                                                                                                    ttcacacgcttagttcatgccgagacgaatcaggtccagtggtgccagtttctcttcacaaccttcaacaactacggtttttccgcggcggatttgctcagc                                                 >  1:1735337‑M1/1‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCAT  >  NC_000913/1629264‑1629559

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