New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 86651837 (1.190)4 (0.140) 4/158 2.7 9.9% coding (96/921 nt) gsiC glutathione ABC transporter permease GsiC
?NZ_CP009273 = 866536 39 (1.370)coding (114/921 nt) gsiC glutathione ABC transporter permease GsiC
Rejected: Frequency below/above cutoff threshold.

CTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACT  >  NZ_CP009273/866446‑866624
                                                                                          |                                                                                        
cTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCgattgatt                                                                                          >  1:525106/1‑91 (MQ=255)
  tACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgc                                                                                        <  2:365398‑M2/91‑3 (MQ=255)
  tACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgc                                                                                        <  1:18999‑M2/91‑3 (MQ=255)
   aCTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgcc                                                                                       >  1:661526‑M2/1‑88 (MQ=255)
        ggTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgg                                                                                     >  1:183545‑M2/1‑83 (MQ=255)
        ggTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgg                                                                                     <  2:183545‑M2/88‑6 (MQ=255)
            gATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGaa                                                                              >  1:409549/1‑91 (MQ=255)
              ttCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggca                                                                                  >  2:600946‑M2/1‑77 (MQ=255)
              ttCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggca                                                                                  <  1:600946‑M2/85‑9 (MQ=255)
                ccGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTg                                                                          >  1:571958/1‑91 (MQ=255)
                      cTGTTTATCGTCTGGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCGGGCGACCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  2:298428/91‑1 (MQ=255)
                      cTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  2:493921/91‑1 (MQ=255)
                      cTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  2:645705/91‑1 (MQ=255)
                      cTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  1:461059/91‑1 (MQ=255)
                      cTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  1:108148/91‑1 (MQ=255)
                      cTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  1:236985/91‑1 (MQ=255)
                      cTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  1:244056/91‑1 (MQ=255)
                        gTTTATCGTCTCGGTGCTGGTTTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAg                                                                  <  2:160390/91‑1 (MQ=255)
                        gTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAg                                                                  <  2:511070/91‑1 (MQ=255)
                        gTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAg                                                                  <  1:196315/91‑1 (MQ=255)
                          ttATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggcagcatatggacaaataaaa                                                                >  1:453045‑M2/1‑65 (MQ=255)
                            aTCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggcagcatatggacaaataaaaac                                                              >  1:457900‑M2/1‑63 (MQ=255)
                                  tcGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAAc                                                        >  1:290469/1‑91 (MQ=255)
                                    ggTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTg                                                      <  2:525106/91‑1 (MQ=255)
                                               ttttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAg                                                                  >  2:764829/1‑68 (MQ=255)
                                               ttttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAg                                                                  <  1:764829/68‑1 (MQ=255)
                                               ttttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTcagca                                           >  2:10186/1‑91 (MQ=255)
                                                  tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCt                                        <  2:114501/91‑1 (MQ=255)
                                                             atGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATc                             <  1:863047/91‑1 (MQ=255)
                                                                ctgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGc                          >  1:450070/1‑91 (MQ=255)
                                                                 tgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCt                                        <  1:805506/76‑1 (MQ=255)
                                                                 tgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCt                                        >  2:805506/1‑76 (MQ=255)
                                                                 tgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGcc                         >  2:538488/1‑91 (MQ=255)
                                                                 tgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGcc                         >  1:535192/1‑91 (MQ=255)
                                                                 tgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGcc                         >  2:801434/1‑91 (MQ=255)
                                                                  gcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCg                        >  2:659007/1‑91 (MQ=255)
                                                                  gcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCg                        <  1:182859/91‑1 (MQ=255)
                                                                  gcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCg                        >  2:220238/1‑91 (MQ=255)
                                                                  gcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCg                        >  1:491176/1‑91 (MQ=255)
                                                                   cccGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    >  1:754722/1‑46 (MQ=255)
                                                                   cccGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATgcgc                                                                    <  2:754722/46‑1 (MQ=255)
                                                                            ccGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCAcc              >  2:817764/1‑91 (MQ=255)
                                                                             cGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCa             <  2:217910/91‑1 (MQ=255)
                                                                                   gattgattGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTg       >  2:144043/1‑91 (MQ=255)
                                                                                   gattgattGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTg       <  1:45508/91‑1 (MQ=255)
                                                                                    attgattgCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTgg      >  2:284460/1‑91 (MQ=255)
                                                                                     ttgattgCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGc     <  2:767279/91‑1 (MQ=255)
                                                                                     ttgattgCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGc     <  1:365398/91‑1 (MQ=255)
                                                                                        attgCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACt  <  2:317292/91‑1 (MQ=255)
                                                                                          |                                                                                        
CTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGGGTTGGATCAGCCGCTGTATCACCAGTTCTGGCACT  >  NZ_CP009273/866446‑866624

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.