New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 369327123 (1.050)3 (0.150) 3/152 1.6 12.3% intergenic (+18/+69) bcsG/ldrD cellulose biosynthesis protein BcsG/type I toxin‑antitoxin system toxic polypeptide LdrD
?NZ_CP009273 = 3693281 22 (1.130)intergenic (+28/+59) bcsG/ldrD cellulose biosynthesis protein BcsG/type I toxin‑antitoxin system toxic polypeptide LdrD
Rejected: Frequency below/above cutoff threshold.

AATACCAGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGTTCAGCCA  >  NZ_CP009273/3693192‑3693360
                                                                               |                                                                                         
aaTACCAGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTccccc                                                                                 >  2:90176/1‑90 (MQ=255)
  tACCAGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCct                                                                               <  1:113432‑M2/90‑3 (MQ=255)
    ccAGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgc                                                                             >  2:375787‑M2/1‑86 (MQ=255)
      aGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTcccc                                                                                  <  2:121994/83‑1 (MQ=255)
      aGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTcccc                                                                                  >  1:121994/1‑83 (MQ=255)
      aGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgcgg                                                                           <  2:560238‑M2/90‑7 (MQ=255)
      aGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgcgg                                                                           <  2:534698‑M2/90‑7 (MQ=255)
      aGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgcgg                                                                           <  1:556196‑M2/90‑7 (MQ=255)
         aTAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGGGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgcggcct                                                                        <  1:450978‑M2/90‑10 (MQ=255)
           aaaCCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTgg                                                                      >  1:380722/1‑90 (MQ=255)
           aaaCCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTgg                                                                      >  2:528004/1‑90 (MQ=255)
           aaaCCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTgg                                                                      >  2:270221/1‑90 (MQ=255)
           aaaCCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTgg                                                                      >  2:155456/1‑90 (MQ=255)
           aaaCCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTgg                                                                      >  2:127946/1‑90 (MQ=255)
            aaCCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGt                                                                     >  2:513382/1‑90 (MQ=255)
              ccGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgcggcctttttc                                                                   >  2:473480‑M2/1‑76 (MQ=255)
                gTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctg                                                                              >  2:597660‑M2/1‑74 (MQ=255)
                gTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctg                                                                              <  1:597660‑M2/77‑4 (MQ=255)
                  aCGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGgc                                                               >  2:134780/1‑90 (MQ=255)
                  aCGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGgc                                                               >  2:102444/1‑90 (MQ=255)
                    gTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGgcgc                                                             <  2:35207/90‑1 (MQ=255)
                      tCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGggggggaactctgcggcctttttcgttctttc                                                           >  2:156034‑M1/1‑58 (MQ=255)
                       cGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCaaa                                                          <  1:270221/90‑1 (MQ=255)
                       cGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCaaa                                                          <  2:512004/90‑1 (MQ=255)
                              aCGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCctgcggcctttttcgttcttactgcgggta                                                   >  1:472904‑M2/1‑60 (MQ=255)
                                   ggcgACTGGGTGCCTTCCCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCCGGTCCGGGCGCAAATTGCCATTACgg                                              <  1:281973/90‑1 (MQ=255)
                                                  tACCCGCAGTAAGAACGAAAAAGGCCGCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgt                               >  1:472904‑M1/1‑30 (MQ=255)
                                                  tACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGc                               <  1:514545/90‑1 (MQ=255)
                                                       gCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGac                          <  1:8223/90‑1 (MQ=255)
                                                       gCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGac                          <  1:79826/90‑1 (MQ=255)
                                                          gaaagaACGAAAAAGGCCGCAGAGTTCCCCCCctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcga                       >  2:156034‑M2/3‑32 (MQ=255)
                                                             aGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTa                    >  2:131783/1‑90 (MQ=255)
                                                                  gAAAAAGGCCGCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtacgg               >  2:473480‑M1/1‑14 (MQ=255)
                                                                  gAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTcc               >  2:541861/1‑90 (MQ=255)
                                                                       aGGCCGCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcccccagtcgccgccgttcaggcgaacgtacggtttat          <  1:450978‑M1/90‑82 (MQ=255)
                                                                        ggCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGt         >  2:579640/1‑90 (MQ=255)
                                                                          ccGCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtacggtttatcct       <  1:556196‑M1/90‑85 (MQ=255)
                                                                          ccGCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtacggtttatcct       <  2:534698‑M1/90‑85 (MQ=255)
                                                                          ccGCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtacggtttatcct       <  2:560238‑M1/90‑85 (MQ=255)
                                                                            gCAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtacggtttatcctgg     >  2:375787‑M1/1‑4 (MQ=255)
                                                                             cAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtac                 <  1:597660‑M1/77‑75 (MQ=255)
                                                                             cAGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtac                 >  2:597660‑M1/1‑3 (MQ=255)
                                                                              agagTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGTTCAGcc   <  2:96329/90‑1 (MQ=255)
                                                                              aGgggggaaactctgcggcctttttcgttcttactgcgggtaaggcacccagtcgccgccgttcaggcgaacgtacggtttatcctggta   <  1:113432‑M1/90‑89 (MQ=255)
                                                                               gagCTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGTTCAGCCa  >  1:194118/1‑90 (MQ=255)
                                                                               |                                                                                         
AATACCAGGATAAACCGTACGTTCGCCTGAACGGCGGCGACTGGGTGCCTTACCCGCAGTAAGAACGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGTTCAGCCA  >  NZ_CP009273/3693192‑3693360

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