New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1661537 =30 (0.570)3 (0.060) 3/158 3.5 7.5% intergenic (‑61/+64) bioD/mlc dethiobiotin synthase/DNA‑binding transcriptional repressor Mlc
?NZ_CP009273 3358915 = 45 (0.900)intergenic (+3/+146) BW25113_RS16670/BW25113_RS16680 YhcE family protein/IS5‑like element IS5 family transposase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACA  >  NZ_CP009273/1661448‑1661626
                                                                                         |                                                                                         
gTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTc                                                                                           <  2:163967/90‑1 (MQ=255)
gTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTc                                                                                           <  2:48290/90‑1 (MQ=255)
 tACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCt                                                                                          >  2:1161851/1‑90 (MQ=255)
     ggTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGt                                                                                      >  1:1266512/1‑90 (MQ=255)
     ggTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGt                                                                                      >  2:1031388/1‑90 (MQ=255)
          taaAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGc                                                                                 <  1:1067199/90‑1 (MQ=255)
                aaCGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCt                                                                           >  2:311920/1‑90 (MQ=255)
                     ttCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAATCCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGccc                                                                      <  1:1471765/90‑1 (MQ=255)
                     ttCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGccc                                                                      <  2:1283284/90‑1 (MQ=255)
                     ttCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGccc                                                                      <  2:1449570/90‑1 (MQ=255)
                      tCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCa                                                                     >  2:500476/1‑90 (MQ=255)
                                        ccGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAAt                                                   >  2:1002805/1‑90 (MQ=255)
                                            tATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTg                                               >  2:664760/1‑90 (MQ=255)
                                             aTGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTgg                                              >  2:545519/1‑90 (MQ=255)
                                             aTGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTgg                                              >  2:1027970/1‑90 (MQ=255)
                                             aTGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTgg                                              >  2:1043504/1‑90 (MQ=255)
                                               gCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTa                                            >  1:1492009/1‑90 (MQ=255)
                                                  tCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGaa                                         >  1:558691/1‑90 (MQ=255)
                                                   cacaAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAAc                                        <  1:243603/90‑1 (MQ=255)
                                                   cacaAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAAc                                        <  2:1245314/90‑1 (MQ=255)
                                                    acaAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACa                                       >  2:112636/1‑90 (MQ=255)
                                                       aaTATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGtt                                    >  1:1291588/1‑90 (MQ=255)
                                                           taAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTaaaa                                <  2:1492009/90‑1 (MQ=255)
                                                            aaaCCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTaaaaa                               <  1:545519/90‑1 (MQ=255)
                                                                 aGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGtt                          >  2:216514/1‑90 (MQ=255)
                                                                 aGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGtt                          >  2:863437/1‑90 (MQ=255)
                                                                  ggAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCCGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAAc                                        >  1:882574/1‑75 (MQ=255)
                                                                  ggAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCCGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAAc                                        <  2:882574/75‑1 (MQ=255)
                                                                      aataatTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAAccc                     <  2:341869/90‑1 (MQ=255)
                                                                                     ttcccTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATc      >  1:862469/5‑90 (MQ=255)
                                                                                        tCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATcaac   >  1:245305/1‑90 (MQ=255)
                                                                                         cTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATcaaca  >  1:716440/1‑90 (MQ=255)
                                                                                         |                                                                                         
GTACCGGTAATAAAGAAACGCTTCAGCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACA  >  NZ_CP009273/1661448‑1661626

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