New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 253246453 (1.260)3 (0.070) 3/272 6.3 5.7% coding (58/972 nt) cysK cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
?NC_000913 = 2532492 48 (1.190)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/2532351‑2532631
                                                                                                                                             |                                                                                                                                           
ctctGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCCGTAACGGACGACTTCTg                                                                                                                                            >  2:109444/1‑143 (MQ=255)
 tctGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTgg                                                                                                                                           >  2:638986/1‑143 (MQ=255)
    gAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGa                                                                                                                                        <  2:341575/143‑1 (MQ=255)
         ggTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGcccaaggc                                                                                                                                   >  1:507925/1‑136 (MQ=255)
         ggTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTg                                                                                                                                   <  1:203902/143‑1 (MQ=255)
              gCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAAtc                                                                                                                              >  1:565557/1‑143 (MQ=255)
               cTACCTGTTGTATCCCAATTTCATACAGGTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGGTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGGGGAAtct                                                                                                                             >  2:522623/1‑143 (MQ=255)
                tACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTCACGGACGCATTCTGGCGAAGGTGGAAtctc                                                                                                                            >  2:283464/1‑143 (MQ=255)
                     gttgtATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAAc                                                                                                                       >  2:699023/1‑143 (MQ=255)
                     gttgtATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAAc                                                                                                                       >  2:1045/1‑143 (MQ=255)
                           ttccAATTTCATGCAGTTAAGGACAGGCCATGTGTAAGATTTTTGAAGATACCTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGc                                                                                                                 <  1:533012/141‑1 (MQ=255)
                           tCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCCTTCTGGCGAAGGTGGAATCTCGTAACCCCAGc                                                                                                                 >  2:661850/1‑143 (MQ=255)
                           tCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTGCGCCTGAATCGCATCGGTAACGGACGCCTTCTGGAGAAGGTGGAATCGCGTAACCCCCGc                                                                                                                 >  2:422296/1‑143 (MQ=255)
                             ccAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGCAGGTGGAATCTAGTAACCCCAGcct                                                                                                               >  1:675360/1‑141 (MQ=255)
                             ccAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCtt                                                                                                               >  1:194869/1‑143 (MQ=255)
                               aaTTTCATACAGTTAGGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCGGCGTCa                                                                                                             >  1:406448/1‑143 (MQ=255)
                                 ttccATACAGTTAAGGACAGGCCATGAGTAAGCTTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGc                                                                                                           <  2:319102/140‑1 (MQ=255)
                                    cATGCAGTTAAGGACAGGCAAGGAGCAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGtt                                                                                                        <  1:513616/143‑1 (MQ=255)
                                     acaCATTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATACCTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTa                                                                                                       <  1:560086/141‑1 (MQ=255)
                                           ttAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGTACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGGTAAGTGcc                                                                                                 >  1:128545/1‑143 (MQ=255)
                                           tgaaGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGcc                                                                                                 <  1:186156/141‑1 (MQ=255)
                                            tAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCg                                                                                                <  2:159926/143‑1 (MQ=255)
                                             caGGACGGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGt                                                                                               <  1:66373/142‑1 (MQ=255)
                                                gACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGCTACCGGACGCATTCTGGCGAAAGTGGAATCTCGTACCCCCAGCCCCAGCgctgcgtgcagcccg                                                                                            >  1:673122/1‑129 (MQ=255)
                                                 ccAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGCCACACGCCGCTGGTTCGCCTGAATCGCATCGGGAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCg                                                                                           <  1:200650/142‑1 (MQ=255)
                                                  cAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCgg                                                                                          <  2:352531/143‑1 (MQ=255)
                                                     gCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTACCCCAGCCTCAGCGTTAAGTGCCGTATCGGTGcc                                                                                      >  1:143231/1‑143 (MQ=255)
                                                        ctGAGGAAGATGTTTGAGGGTAACTCGCGGGCGGGCGGTCACACGCCGCTGGTTCGCCGGAGTCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCGCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCaa                                                                                    <  1:12238/142‑1 (MQ=255)
                                                               aGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGagt                                                                             >  2:107830/1‑141 (MQ=255)
                                                                 aTTTTGGAGGATAACTCGCTGACTATCTGTCACACGCCGCTGGTTCGGCTGAATCGAGCCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTg                                                                           <  1:352474/143‑1 (MQ=255)
                                                                  tttttGACGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCCTCGGTAGCGGTTGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTgg                                                                          <  2:535033/143‑1 (MQ=255)
                                                                  tttttGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTgg                                                                          <  2:433371/143‑1 (MQ=255)
                                                                        aaGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGCAGGTGGAATCTCGTAACCCCAGCTTCCGCGTGAAGTGCCGTAACGGTGCCCACACGCATCGGGATGcc                                                                    >  1:498472/1‑143 (MQ=255)
                                                                                tCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAgcg                                                            >  2:285237/1‑143 (MQ=255)
                                                                                  gCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGcgcg                                                          <  1:283464/143‑1 (MQ=255)
                                                                                             gTCACGCGTCGTTGGTTCCGCTGAATCGCATAGGGAACGGACGCATCCGGGCGCAGGTGGAATCTCGCAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGaaa                                               <  1:152791/143‑1 (MQ=255)
                                                                                              tCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAGGCGCGACGTGCTGacgc                                              >  1:386514/1‑140 (MQ=255)
                                                                                                  acGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTTGAAACAgg                                          >  2:267529/1‑143 (MQ=255)
                                                                                                  acGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAgg                                          >  1:731486/1‑143 (MQ=255)
                                                                                                     ccgcTGGGTGGCCTGACTCGCATCGGTAACGGACGCATTCGGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATGGGTGTCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGt                                       <  2:519025/143‑1 (MQ=255)
                                                                                                      cgcCGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGACAGTGGAAGCTCGTAACCCCAGCTCCAGCGTTAACTGCCGGATCACTGCCAACATCATTGACGACGCCGAAAATCGCGACGTGCAGAAAGCAcgcgct                                      >  1:118091/1‑141 (MQ=255)
                                                                                                           gTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTACGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACt                                 <  2:597758/143‑1 (MQ=255)
                                                                                                            ttCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGGTGAACTg                                >  1:597842/1‑143 (MQ=255)
                                                                                                                   gatccgcgTCGGTACCGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAAc                         <  2:425068/139‑1 (MQ=255)
                                                                                                                          ctCGGTAGCGGACGCATTCTGCCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAACGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAg                  <  2:31357/142‑1 (MQ=255)
                                                                                                                              gTATCGGACTCACTCTGGCGATGGTGGAATCTCGTTACCCCAGCTTTAGCGTTAAGTGCCGTATCGGTGCTAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGt              <  1:731750/143‑1 (MQ=255)
                                                                                                                               tAACGCACGCATTCTGGCGAAGGTGGAATCGCGTAACCCCAGCGTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTa             <  1:1045/143‑1 (MQ=255)
                                                                                                                                  cGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAGCCGACCAGCGGTGATa          >  2:442536/1‑143 (MQ=255)
                                                                                                                                   ggACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGACAAGCGCGGCGTGCTGAAACCAGGCGCTGAACCATATGAACCCACCAGCGGGACTAc         >  1:475224/1‑143 (MQ=255)
                                                                                                                                      cGCATTCTGGCGAAGGTGGACTCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCgg      <  2:299177/143‑1 (MQ=255)
                                                                                                                                       gCATTCCAGAGAACGTGCAATCTCCTAACCGCATCTTCAGCGTTAAGTGCCTTATCGGTGCCAACATTATTTGAGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCggg     <  2:99501/143‑1 (MQ=255)
                                                                                                                                          ttCTGGCGAAGGTGGAAGCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGAtt  <  1:422296/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
CTCTGAGCGGGTATGCTACCTGTTGTATCCCAATTTCATACAGTTAAGGACAGGCCATGAGTAAGATTTTTGAAGATAACTCGCTGACTATCGGTCACACGCCGCTGGTTCGCCTGAATCGCATCGGTAACGGACGCATTCTGGCGAAGGTGGAATCTCGTAACCCCAGCTTCAGCGTTAAGTGCCGTATCGGTGCCAACATGATTTGGGATGCCGAAAAGCGCGGCGTGCTGAAACCAGGCGTTGAACTGGTTGAACCGACCAGCGGTAATACCGGGATT  >  NC_000913/2532351‑2532631

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.