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.

ACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  NZ_CP009273/3248569‑3248739
                                                                                    |                                                                                      
atcgtgccggatgcgatgctggcgcatcttgtccggcctacgggttgtttaattcgttaatgttgttggattttgtaggccggaCAAGAt                                                                                   <  2:531031‑M2/6‑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)
                         atcttgtccggcctacgggttgtttaattcgttaatgttgttggattttgtaggccggaCAAGATGCGCCAGCATCGCATCCGGCACGAt                                                          <  1:716748‑M2/31‑1 (MQ=255)
                         atcttgTCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGAt                                                          <  2:531031‑M1/84‑1 (MQ=255)
                                   gcctacgggttgtttaattcgttaatgttgttggattttgtaggccggaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAAcc                                                <  2:645363‑M2/41‑1 (MQ=255)
                                   GCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAAcc                                                <  2:240502‑M1/90‑1 (MQ=255)
                                                 acaacaTTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTgg                                  >  2:518110/1‑90 (MQ=255)
                                                  aattcgttaatgttgttggattttgtaggccggaCAAGATGCGCCAGCATCGCATCCGGCACTATCCCCAAAACCTGGCGTTTATCTGGc                                 <  2:581561‑M2/56‑1 (MQ=255)
                                                      aTTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCttc                             >  1:474729/1‑89 (MQ=255)
                                                           cGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGc                        >  1:194816/1‑90 (MQ=255)
                                                              aTTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCAt                     >  1:190794/1‑90 (MQ=255)
                                                              aTTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCAt                     >  2:96401/1‑90 (MQ=255)
                                                              aTTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCAt                     >  1:521517/1‑90 (MQ=255)
                                                                      acCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGa             >  2:338367/1‑90 (MQ=255)
                                                                      acCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGa             >  1:273937/1‑90 (MQ=255)
                                                                           taggccagaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATc        >  2:443437‑M2/10‑90 (MQ=255)
                                                                                 ggaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGa  >  2:253702‑M2/4‑90 (MQ=255)
                                                                                 ggaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGa  >  1:86569‑M2/4‑90 (MQ=255)
                                                                                 ggaCAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGa  >  1:536738‑M2/4‑90 (MQ=255)
                                                                                    |                                                                                      
ACAGGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCCAACAACATTAACGAATTAAACAACCCGTAGGCCGGACAAGATGCGCCAGCATCGCATCCGGCACGATCCCCAAAACCTGGCGTTTATCTGGCCTTGCTCGCCATAATCTCGATAATCTGCCGA  >  NZ_CP009273/3248569‑3248739

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