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.

CTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATT  >  NC_000913/2532353‑2532633                                                                                                                                              |                                                                                                                                            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)                                                                                                                    gatccgcgTCGGTACCGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAAc                         <  1:425068/139‑1 (MQ=255)                                                                                                                           ctCGGTAGCGGACGCATTCTGCCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAACGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAg                  <  1:31357/142‑1 (MQ=255)                                                                                                                               gTATCGGACTCACTCTGGCGATGGTGGAATCTCGTTACCCCAGCTTTAGCGTTAAGTGCCGTATCGGTGCTAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGt              <  2:731750/143‑1 (MQ=255)                                                                                                                                tAACGCACGCATTCTGGCGAAGGTGGAATCGCGTAACCCCAGCGTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTa             <  2:1045/143‑1 (MQ=255)                                                                                                                                   cGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAGCCGACCAGCGGTGATa          >  1:442536/1‑143 (MQ=255)                                                                                                                                    ggACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGACAAGCGCGGCGTGCTGAAACCAGGCGCTGAACCATATGAACCCACCAGCGGGACTAc         >  2:475224/1‑143 (MQ=255)                                                                                                                                       cGCATTCTGGCGAAGGTGGACTCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCgg      <  1:299177/143‑1 (MQ=255)                                                                                                                                        gCATTCCAGAGAACGTGCAATCTCCTAACCGCATCTTCAGCGTTAAGTGCCTTATCGGTGCCAACATTATTTGAGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCggg     <  1:99501/143‑1 (MQ=255)                                                                                                                                           ttCTGGCGAAGGTGGAAGCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGAtt  <  2:422296/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                            CTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATT  >  NC_000913/2532353‑2532633

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.