New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3248600 =22 (0.730)4 (0.150) 4/148 1.7 19.3% intergenic (+34/+100) BW25113_RS16135/cyuA hypothetical protein/cysteine desulfidase
?NZ_CP009273 3248653 = 14 (0.530)intergenic (+87/+47) BW25113_RS16135/cyuA hypothetical protein/cysteine desulfidase
Rejected: Frequency below/above cutoff threshold.

AGCTATCGTCGCGGATGCTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGAT  >  NZ_CP009273/3248512‑3248683
                                                                                        |                                                                                   
aGCTATCGTCGCGGATGCTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATc                                                                                    <  2:749792/90‑1 (MQ=255)
 tcggcagattatcgagattatggcgagcaaggccagataaacgccaggttttggggatcgtgccggatgcgatgctggcgcatcttgTcc                                                                                   >  1:536738‑M1/88‑90 (MQ=255)
 tcggcagattatcgagattatggcgagcaaggccagataaacgccaggttttggggatcgtgccggatgcgatgctggcgcatcttgTcc                                                                                   >  1:86569‑M1/88‑90 (MQ=255)
 tcggcagattatcgagattatggcgagcaaggccagataaacgccaggttttggggatcgtgccggatgcgatgctggcgcatcttgTcc                                                                                   >  2:253702‑M1/88‑90 (MQ=255)
   tATCGTCGCGGATGCTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCgg                                                                                 >  2:311664/1‑90 (MQ=255)
       gattatcgagattatggcgagcaaggccagataaacgccaggttttggggatcgtgccggatgcgatgctggcgcatcttgTCTGGCCTa                                                                             >  2:443437‑M1/82‑90 (MQ=255)
         cgcgGATGCTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACa                                                                           <  2:134634/90‑1 (MQ=255)
          gcgGATGCTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACaa                                                                          <  1:130396/90‑1 (MQ=255)
                 cTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCcaa                                                                   <  2:404366/90‑1 (MQ=255)
                  tGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCcaac                                                                  <  1:594638/90‑1 (MQ=255)
                          gTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTaa                                                          <  2:131623/90‑1 (MQ=255)
                                gccagataaacgccaggttttggggatagtgccggatgcgatgctggcgcatcttgTCCGGCCTACAAAATCCAACAACATTAACGAAtt                                                    <  2:581561‑M1/34‑1 (MQ=255)
                                  gAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTaa                                                  <  2:558615/90‑1 (MQ=255)
                                     gATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAaaca                                               <  2:849549/90‑1 (MQ=255)
                                     gATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAaaca                                               <  2:769281/90‑1 (MQ=255)
                                      aTGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAaacaa                                              <  1:633021/90‑1 (MQ=255)
                                        gAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAAcc                                            >  2:288209/1‑90 (MQ=255)
                                               ggttttggggatcgtgccggatgcgatgctggcgcatcttgTCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGc                                     <  2:645363‑M1/49‑1 (MQ=255)
                                                 gcgTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGcc                                                                               >  1:846149/1‑46 (MQ=255)
                                                 gcgTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGcc                                                                               <  2:846149/46‑1 (MQ=255)
                                                         atcgtgccggatgcgatgctggcgcatcttgTCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGAt                           <  1:716748‑M1/59‑1 (MQ=255)
                                                           aGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATgc                         >  2:729710/1‑90 (MQ=255)
                                                             gTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATgcgc                       >  1:320687/1‑90 (MQ=255)
                                                              ttGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGcc                      <  2:704560/90‑1 (MQ=255)
                                                               tGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCa                     <  2:77061/90‑1 (MQ=255)
                                                                gCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAg                    <  2:847056/90‑1 (MQ=255)
                                                                 ccGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAgc                   >  2:423816/1‑90 (MQ=255)
                                                                   ggATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAaaca                                               >  2:179921/1‑60 (MQ=255)
                                                                   ggATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAaaca                                               <  1:179921/60‑1 (MQ=255)
                                                                                  atcttgTCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGAt  <  2:531031‑M1/84‑1 (MQ=255)
                                                                                        |                                                                                   
AGCTATCGTCGCGGATGCTGAAACGCGTTTAGCCGAGGATGAAGCTACCGCGTAAAGACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGAT  >  NZ_CP009273/3248512‑3248683

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