New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 73442823 (0.710)3 (0.100) 3/158 2.5 12.7% intergenic (+122/‑29) ybfD/ybgA ISAs1 family transposase/YbgA family protein
?NZ_CP009273 734433 = 19 (0.610)intergenic (+127/‑24) ybfD/ybgA ISAs1 family transposase/YbgA family protein
Rejected: Frequency below/above cutoff threshold.

GTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCACCCTAATCCGTATCTTTGCCCTTCATGAGTTACA  >  NZ_CP009273/734345‑734518
                                                                                        |                                                                                     
gTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTcc                                                                                      <  2:313003/90‑1 (MQ=255)
     tGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTaa                                                                                 <  1:923600/90‑1 (MQ=255)
                  cTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                                                                    <  1:205073/90‑1 (MQ=255)
                     tCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAAc                                                                               >  2:889994/1‑76 (MQ=255)
                     tCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAAc                                                                               <  1:889994/76‑1 (MQ=255)
                        ttcgtaaatttgtcctgctacaacaggattaacttcacaaatatcatttctcaacgtctacactCCTGTAAACCGCTCAGGAGCAGTAAt                                                              <  2:401516‑M2/26‑1 (MQ=255)
                         tcgtaaatttgtcctgctacaacaggattaacttcacaaatatcatttctcaacgtctacactCCTGTAAACCGCTCAGGAGCAGTa                                                                <  1:714634‑M2/24‑1 (MQ=255)
                         tcgtaaatttgtcctgctacaacaggattaacttcacaaatatcatttctcaacgtctacactCCTGTAAACCGCTCAGGAGCAGTa                                                                >  2:714634‑M2/64‑87 (MQ=255)
                         aaaTTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATg                                                             <  1:216970/90‑1 (MQ=255)
                                    gCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACg                                                  >  2:42619/1‑90 (MQ=255)
                                              gATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGAc                                        <  2:547658/90‑1 (MQ=255)
                                                ttAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAg                                      <  2:259622/90‑1 (MQ=255)
                                                 tAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGc                                     <  1:446889/90‑1 (MQ=255)
                                                     ttCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCAccc                                 <  2:73912/90‑1 (MQ=255)
                                                              atCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGc                                     <  1:498882/77‑1 (MQ=255)
                                                              atCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGc                                     >  2:498882/1‑77 (MQ=255)
                                                                 aTTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCACCCTAATCCGTATCt                     >  2:86288/1‑90 (MQ=255)
                                                                         aCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGc                                     >  1:548381/1‑66 (MQ=255)
                                                                         aCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGc                                     <  2:548381/66‑1 (MQ=255)
                                                                           gTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCACCCTAATCCGTATCTTTGCCCTTCa           <  2:492783/90‑1 (MQ=255)
                                                                                    tACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCACCCTAATCCGTATCTTTGCCCTTCATGAGTTACa  >  2:342883/1‑90 (MQ=255)
                                                                                    tACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCACCCTAATCCGTATCTTTGCCCTTCATGAGTTACa  >  2:884215/1‑90 (MQ=255)
                                                                                        |                                                                                     
GTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAACGATTTGATGACAGCACCCTAATCCGTATCTTTGCCCTTCATGAGTTACA  >  NZ_CP009273/734345‑734518

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