New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4621981	65 (0.760)	5 (0.060)	5/166	NT	8.2%	coding (213/1332 nt)	yjjJ	putative protein kinase
	?	NC_000913	= 4622003	52 (0.650)					coding (235/1332 nt)	yjjJ	putative protein kinase
Rejected: Frequency below cutoff threshold.

GGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGAC  >  NC_000913/4621906‑4622102                                                                                                  |                                                                                                    ggCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACAtc                                                                                                   <  1:348407‑M2/100‑3 (MQ=255)   cAACGCTATATGCACTGCTGCGTCCTTATTGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACAtccg                                                                                                 >  1:23075‑M2/1‑96 (MQ=255)   cAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACAtccg                                                                                                 <  2:1158343‑M2/100‑5 (MQ=255)    aaCGCGATACGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACAtccgg                                                                                                >  2:502180‑M2/1‑95 (MQ=255)     aCGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACAtccggt                                                                                               >  1:95699‑M2/1‑94 (MQ=255)     aCGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACAtccggt                                                                                               <  1:1542588‑M2/100‑7 (MQ=255)      cgcgATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCCTAAATTCGCCGACATCCGGtt                                                                                              >  1:745788/1‑100 (MQ=255)       gcgATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTtg                                                                                             <  1:1524052/100‑1 (MQ=255)         gATATGCACTGCTGCGTCCTTATCACGGAATTGAGCGTATTCCCGTCTGGCGGGTGGTCCGTTCCGGAAAGGCGCATAAATTCGCCGACCTCCCGTtgtg                                                                                           >  2:1665319/1‑100 (MQ=255)          atatGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGc                                                                                          <  1:1868183/100‑1 (MQ=255)          atatGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGc                                                                                          <  1:1642052/100‑1 (MQ=255)               cACTGCTGCGTCCTTATCGCGGGATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGcc                                                                                     <  1:569166/100‑1 (MQ=255)                aCTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCg                                                                                    <  2:1618093/100‑1 (MQ=255)                 ctgctgCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGc                                                                                   >  2:1805816/1‑100 (MQ=255)                   gctgtGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAg                                                                                 <  1:119917/100‑1 (MQ=255)                   gctgcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAg                                                                                 >  1:1477973/1‑100 (MQ=255)                    ctgcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAgg                                                                                <  1:1869104/100‑1 (MQ=255)                     tgcGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAggg                                                                               >  2:95780/1‑100 (MQ=255)                        gTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAg                                                                            <  1:2238499/100‑1 (MQ=255)                          ccTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGtt                                                                          >  2:1086798/1‑100 (MQ=255)                              aTCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCt                                                                      <  1:1147460/100‑1 (MQ=255)                                cgcgGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTg                                                                     <  2:336181/99‑1 (MQ=255)                                   ggAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTaa                                                                 >  1:1439619/1‑100 (MQ=255)                                    gAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAAc                                                                <  2:833588/100‑1 (MQ=255)                                    gAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAAc                                                                <  2:883603/100‑1 (MQ=255)                                          aGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGgcgc                                                          >  2:830131/1‑100 (MQ=255)                                           gCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCa                                                         <  1:1158349/100‑1 (MQ=255)                                           gCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCa                                                         <  1:1462242/100‑1 (MQ=255)                                             gTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGa                                                       <  2:2286661/100‑1 (MQ=255)                                             gTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGa                                                       >  1:1527956/1‑100 (MQ=255)                                                  cccGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcg                                                  >  1:1179034/1‑100 (MQ=255)                                                   ccGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcga                                                 <  2:1913163/100‑1 (MQ=255)                                                   ccGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcga                                                 >  1:567488/1‑100 (MQ=255)                                                    cGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcgac                                                <  1:1504273/100‑1 (MQ=255)                                                     gtatGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcgacg                                               <  2:1516371/97‑1 (MQ=255)                                                     gTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcgacg                                               <  1:1805831/100‑1 (MQ=255)                                                      tCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGcgacga                                              >  1:581018/1‑100 (MQ=255)                                                       cTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGaa                                             >  1:970336/1‑100 (MQ=255)                                                        tGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAAc                                            >  1:1729429/1‑100 (MQ=255)                                                        tGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAAc                                            <  2:2175969/100‑1 (MQ=255)                                                           cGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGt                                         >  1:431009/1‑100 (MQ=255)                                                            gtgTGGACGATACCGTAAAGGCGCATATAGTCGCCGACATCCGGTTGTGATGTCCGCAGGGGAGTTGTCGGGTAACAGACGCAGATGGcgacg                                               <  2:116633/91‑1 (MQ=255)                                                            gggTGGACGATACCGGAAAGGCGCATTAATTCGCCGACATCCGGTTGTGCTTGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAAc                                            <  2:821548/96‑1 (MQ=255)                                                                ggACGATACCGGAAAGGCGCATAAATTCGCCGACAtccggtatcgtccacccgccagacgggaatacgctcaattccgcgataaggacgcagcagtgcat                                    <  1:1428213‑M2/100‑66 (MQ=255)                                                                 gACGATACCGGAAAGGCGCATAAATTCGCCGACAtccggtatcgtccacccgccagccgggaatacgctcaattccgcgataaggacgcagcagtgcata                                   <  2:2377266‑M2/100‑67 (MQ=255)                                                                 gACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTg                                   >  1:623668/1‑100 (MQ=255)                                                                 gACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTg                                   >  1:419962/1‑100 (MQ=255)                                                                   cGATACCGGAAAGGCGCATAAATTCGCCGACAtccggtatcgtccacccgccagacgggaatacgctcaattccgcgataaggacgccgcagtgcatatc                                 <  2:1822413‑M2/100‑69 (MQ=255)                                                                   cGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGAt                                 >  2:1906035/1‑100 (MQ=255)                                                                    gATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATg                                <  2:373763/100‑1 (MQ=255)                                                                     aTACCGGAAAGGCGCATAAATTCGCCGACAtccggtatcgtccacccgccagacgggaatacgctcaattccgcgataaggacgcagcagtgcatatcgc                               <  1:2178890‑M2/100‑71 (MQ=255)                                                                       aCCGGAAAGGCGCATAAATTCGCCGACAtccggtatcgtccacccgccagacgggaatacgctcaattccgcgataaggacgcagcagtgcatatcgcgt                             <  1:133693‑M2/100‑73 (MQ=255)                                                                       aCCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGtt                             <  2:1092868/100‑1 (MQ=255)                                                                          ggAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGc                          <  2:1729176/100‑1 (MQ=255)                                                                          ggAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGc                          <  1:900519/100‑1 (MQ=255)                                                                               ggCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTgg                     <  1:155713/100‑1 (MQ=255)                                                                                 cgcATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTa                   >  1:393158/1‑100 (MQ=255)                                                                                  gcATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTAt                  >  1:552690/1‑100 (MQ=255)                                                                                   cATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTAtt                 >  1:1911629/1‑100 (MQ=255)                                                                                     taccttCGCCGACATCCGGTTGCGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGTTAttt                <  2:790737/95‑1 (MQ=255)                                                                                     tAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTg               <  1:1394036/100‑1 (MQ=255)                                                                                       aaTTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGAc             >  1:2137752/1‑100 (MQ=255)                                                                                        aTTCGCCGATATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGAc             <  1:82849/99‑1 (MQ=255)                                                                                        aTTCGCCGACATCTGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGAcc            <  2:911020/100‑1 (MQ=255)                                                                                           cGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAt         >  1:827061/1‑100 (MQ=255)                                                                                             ccGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAtct       <  2:1332162/100‑1 (MQ=255)                                                                                             ccGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAtct       >  1:821327/1‑100 (MQ=255)                                                                                             ccGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAtc        <  2:770633/99‑1 (MQ=255)                                                                                              cGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGAtctc      >  1:1370321/1‑100 (MQ=255)                                                                                                aCATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCg    <  2:300005/100‑1 (MQ=255)                                                                                                 cATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGa   >  1:2222420/1‑100 (MQ=255)                                                                                                  aTCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGAc  >  1:978735/1‑100 (MQ=255)                                                                                                  aTCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGAc  >  1:1552167/1‑100 (MQ=255)                                                                                                  |                                                                                                    GGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGAC  >  NC_000913/4621906‑4622102

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