New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2532466	51 (1.260)	3 (0.080)	3/272	6.1	5.9%	coding (58/972 nt)	cysK	cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
	?	NC_000913	= 2532494	47 (1.210)					coding (86/972 nt)	cysK	cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTG  >  NC_000913/2532324‑2532603                                                                                                                                               |                                                                                                                                          gATATTCTAAATCCTTACTTCCTCATATTCTCTGCGCGTGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGcc                                                                                                                                           <  2:221204/143‑1 (MQ=255) gATATGCTAGATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGcc                                                                                                                                           <  2:165109/143‑1 (MQ=255)        tAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCGGTAGAATCCCAATTTCATACAGTTAAGGCCAGGCCATGAGTCAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCg                                                                                                                                    >  1:521753/1‑143 (MQ=255)        tAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATACTACCCGTTTCTTCCCAATTTCATACAGTTAAGAACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCg                                                                                                                                    >  1:502955/1‑143 (MQ=255)         aaaTCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGc                                                                                                                                   >  2:159926/1‑143 (MQ=255)          aaTCCGTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTGACACGCCGCTGGTTCGCCTGAATCGCa                                                                                                                                  <  2:563723/143‑1 (MQ=255)             ccTTACGTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGTCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCg                                                                                                                               <  2:231870/143‑1 (MQ=255)               ttACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGCTTCACAATTTCATACAGTTAAGGACAGGCCATGCGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCggc                                                                                                                             >  2:597758/1‑142 (MQ=255)                 aCTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGACTGAATCGCATCGGTaa                                                                                                                           >  2:302557/1‑143 (MQ=255)                      cGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGAc                                                                                                                      <  2:586823/143‑1 (MQ=255)                        cATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTCACGGACGc                                                                                                                    >  1:74342/1‑143 (MQ=255)                              ctctGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCCGTAACGGACGACTTCTg                                                                                                              >  1:109444/1‑143 (MQ=255)                               tctGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTgg                                                                                                             >  1:638986/1‑143 (MQ=255)                                  gAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGa                                                                                                          <  1:341575/143‑1 (MQ=255)                                       ggTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGcccaaggc                                                                                                     >  2:507925/1‑136 (MQ=255)                                       ggTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTg                                                                                                     <  2:203902/143‑1 (MQ=255)                                            gCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAAtc                                                                                                >  2:565557/1‑143 (MQ=255)                                             cTACCTGTTGTATCCCAATTTCATACAGGTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGGTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGGGGAAtct                                                                                               >  1:522623/1‑143 (MQ=255)                                              tACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTCACGGACGCATTCTGGCGAAGGTGGAAtctc                                                                                              >  1:283464/1‑143 (MQ=255)                                                   gttgtATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAAc                                                                                         >  1:1045/1‑143 (MQ=255)                                                   gttgtATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAAc                                                                                         >  1:699023/1‑143 (MQ=255)                                                         ttccAATTTCATGCAGTTAAGGACAGGCCATGTGTAAGATTTTTGAAGATACCTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGc                                                                                   <  2:533012/141‑1 (MQ=255)                                                         tCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCCTTCTGGCGAAGGTGGAATCTCGTAACCCCAGc                                                                                   >  1:661850/1‑143 (MQ=255)                                                         tCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTGCGCCTGAATCGCATCGGTAACGGACGCCTTCTGGAGAAGGTGGAATCGCGTAACCCCCGc                                                                                   >  1:422296/1‑143 (MQ=255)                                                           ccAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGCAGGTGGAATCTAGTAACCCCAGcct                                                                                 >  2:675360/1‑141 (MQ=255)                                                           ccAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCtt                                                                                 >  2:194869/1‑143 (MQ=255)                                                             aaTTTCATACAGTTAGGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCGGCGTCa                                                                               >  2:406448/1‑143 (MQ=255)                                                               ttccATACAGTTAAGGACAGGCCATGAGTAAGCTTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGc                                                                             <  1:319102/140‑1 (MQ=255)                                                                  cATGCAGTTAAGGACAGGCAAGGAGCAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGtt                                                                          <  2:513616/143‑1 (MQ=255)                                                                   acaCATTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATACCTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTa                                                                         <  2:560086/141‑1 (MQ=255)                                                                         ttAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGTACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGGTAAGTGcc                                                                   >  2:128545/1‑143 (MQ=255)                                                                         tgaaGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGcc                                                                   <  2:186156/141‑1 (MQ=255)                                                                          tAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCg                                                                  <  1:159926/143‑1 (MQ=255)                                                                           caGGACGGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGt                                                                 <  2:66373/142‑1 (MQ=255)                                                                               ccAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGCCACACGCCGCTGGTTCGCCTGAATCGCATCGGGAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCg                                                             <  2:200650/142‑1 (MQ=255)                                                                                cAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCgg                                                            <  1:352531/143‑1 (MQ=255)                                                                                   gCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGT‑ACCCCAGCCTCAGCGTTAAGTGCCGTATCGGTGcc                                                        >  2:143231/1‑143 (MQ=255)                                                                                      ctGAGGAAGATGTTTGAGGGTAACTCGCGGGCGGGCGGTCACACGCCGCTGGTTCGCCGGAGTCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCGCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCaa                                                      <  2:12238/142‑1 (MQ=255)                                                                                             aGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGagt                                               >  1:107830/1‑141 (MQ=255)                                                                                               aTTTTGGAGGATAACTCGCTGACTATCTGTCACACGCCGCTGGTTCGGCTGAATCGAGCCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTg                                             <  2:352474/143‑1 (MQ=255)                                                                                                tttttGACGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCCTCGGTAGCGGTTGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTgg                                            <  1:535033/143‑1 (MQ=255)                                                                                                tttttGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTgg                                            <  1:433371/143‑1 (MQ=255)                                                                                                      aaGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGCAGGTGGAATCTCGTAACCCCAGCTTCCGCGTGAAGTGCCGTAACGGTGCCCACACGCATCGGGATGcc                                      >  2:498472/1‑143 (MQ=255)                                                                                                              tCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAgcg                              >  1:285237/1‑143 (MQ=255)                                                                                                                gCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGcgcg                            <  2:283464/143‑1 (MQ=255)                                                                                                                           gTCACGCGTCGTTGGTTCCGCTGAATCGCATAGGGAACGGACGCATCCGGGCGCAGGTGGAATCTCGCAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGaaa                 <  2:152791/143‑1 (MQ=255)                                                                                                                            tCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAGGCGCGACGTGCTGacgc                >  2:386514/1‑140 (MQ=255)                                                                                                                                acGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTTGAAACAgg            >  1:267529/1‑143 (MQ=255)                                                                                                                                acGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAgg            >  2:731486/1‑143 (MQ=255)                                                                                                                                   ccgcTGGGTGGCCTGACTCGCATCGGTAACGGACGCATTCGGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATGGGTGTCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGt         <  1:519025/143‑1 (MQ=255)                                                                                                                                    cgcCGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGACAGTGGAAGCTCGTAACCCCAGCTCCAGCGTTAACTGCCGGATCACTGCCAACATCATTGACGACGCCGAAAATCGCGACGTGCAGAAAGCAcgcgct        >  2:118091/1‑141 (MQ=255)                                                                                                                                         gTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTACGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACt   <  1:597758/143‑1 (MQ=255)                                                                                                                                          ttCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGGTGAACTg  >  2:597842/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                          GATATGCTAAATCCTTACTTCCGCATATTCTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTG  >  NC_000913/2532324‑2532603

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