New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 253266228 (0.880)3 (0.100) 3/124 2.8 9.2% coding (254/972 nt) cysK cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
?NC_000913 = 2532707 33 (1.100)coding (299/972 nt) cysK cysteine synthase A, O‑acetylserine sulfhydrolase A subunit
Rejected: Frequency below/above cutoff threshold.

CCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGC  >  NC_000913/2532641‑2532774
                                                                  |                                                                   
ccTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGc                                                                    >  1:440338‑M2/1‑67 (MQ=255)
 cTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGcc                                                                   >  1:731520‑M2/1‑66 (MQ=255)
     gTAGCTGCCGCTCGCGGTTGCAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAGCGcggcgg                                                               <  2:522863‑M2/68‑7 (MQ=255)
     gTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCaa                                                               >  1:310552/1‑68 (MQ=255)
     gTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCaa                                                               >  2:748047/1‑68 (MQ=255)
       aGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAGACCATGAGTATTGAACGCCGCAAgc                                                             >  2:171639/1‑68 (MQ=255)
          tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAgctgc                                                          <  1:722479/68‑1 (MQ=255)
          tGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATAGAACGCCGCAAgctgc                                                          >  1:497468/1‑68 (MQ=255)
                  gcgGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCg                                                  <  1:213525/68‑1 (MQ=255)
                   cgGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGt                                                 <  2:1001363/68‑1 (MQ=255)
                          aaaCTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGc                                          >  1:952615/1‑68 (MQ=255)
                            aCTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCaa                                        >  2:829297/1‑68 (MQ=255)
                             cTCACCCTGACCATGCCAGAAACCATGAGTATTGATCGCCGCATGCTGCTGAAAGCTTTAGGTGCaaa                                       <  2:908177/68‑1 (MQ=38)
                              tCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAc                                      <  1:977400/68‑1 (MQ=255)
                               cACCCTGACCATGCCAGAAACCATGAGTATTGAACTCCGCAAGCTGCTGAAAGCGTTAGGTGCAAAcc                                     <  2:1166887/68‑1 (MQ=255)
                                  ccTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTgg                                  >  1:484810/1‑68 (MQ=255)
                                      aCCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCt                              >  1:451088/1‑68 (MQ=255)
                                        cATGCCAGAAACCATGAGTATTGAACGccgcgagcggcagctacataggccagtgcaatcccggtatt                            <  1:287900‑M2/68‑42 (MQ=255)
                                           gCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACgg                         <  1:889068/68‑1 (MQ=255)
                                            ccAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGa                        <  1:900494/68‑1 (MQ=255)
                                              aGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAg                      >  2:84535/1‑68 (MQ=255)
                                                aaaCCATGAGTATTGAACGCCGCAATCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGt                    <  1:1056409/68‑1 (MQ=255)
                                                aaaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGt                    <  2:1000350/68‑1 (MQ=255)
                                                aaaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGt                    <  1:681502/68‑1 (MQ=255)
                                                aaaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGt                    <  1:51429/68‑1 (MQ=255)
                                                 aaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                           >  2:1027659/1‑44 (MQ=255)
                                                 aaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTg                                           <  1:1027660/44‑1 (MQ=255)
                                                 aaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTg                   >  2:113932/1‑68 (MQ=255)
                                                 aaCCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTg                   <  1:419412/68‑1 (MQ=255)
                                                 aaCCATGAGTAGTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTg                   <  2:441396/68‑1 (MQ=255)
                                                  aCCATGAGTAATGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGc                  >  2:722813/1‑68 (MQ=255)
                                                   ccATGAGTATTGAACGccgcgagcggcagctacataggccagtgcaatcccggtattaccgctggtcg                 >  2:287900‑M2/1‑16 (MQ=255)
                                                       gAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAg             <  2:223483/68‑1 (MQ=255)
                                                        aGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAgg            <  1:829297/68‑1 (MQ=255)
                                                        aGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAgg            <  2:981213/68‑1 (MQ=255)
                                                                aCGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAg    >  1:446638/1‑68 (MQ=255)
                                                                aCGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAg    >  1:107330/1‑68 (MQ=255)
                                                                  gccgcAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGgc  >  1:797721/1‑68 (MQ=255)
                                                                  |                                                                   
CCTATGTAGCTGCCGCTCGCGGTTACAAACTCACCCTGACCATGCCAGAAACCATGAGTATTGAACGCCGCAAGCTGCTGAAAGCGTTAGGTGCAAACCTGGTGCTGACGGAAGGTGCTAAAGGCATGAAAGGC  >  NC_000913/2532641‑2532774

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