New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3122373 =38 (1.400)4 (0.170) 3/148 2.1 10.9% coding (743/765 nt) glcC transcriptional regulator GlcC
?NZ_CP009273 3122512 = 32 (1.350)coding (1027/1137 nt) yghO hypothetical protein
Rejected: Frequency below/above cutoff threshold.

CGGTAACGCTTGTACGGCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACTGATAC  >  NZ_CP009273/3122424‑3122597
                                                                                        |                                                                                     
tggaacgcattggggcggagccgtacaagcgttcccgattgtatgaaaagcagatttaataccagctaactcaggttcatctccagcgcc                                                                                      <  1:22176‑M2/2‑2 (MQ=255)
  gTAACGCTTGTACGGCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATc                                                                                    <  2:241802/90‑1 (MQ=255)
  gTAACGCTTGTACGGCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATc                                                                                    <  2:443163/90‑1 (MQ=255)
               ggCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCcggcgg                                                                       <  1:253812/90‑1 (MQ=255)
               cggagccgtacaagcgttaccgattgtatgaaaagcagatttaataccagctaactcaggttcatctccagcgCATCGATCTTCcggcgg                                                                       <  1:558134‑M2/17‑1 (MQ=255)
                gCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCggg                                                                      >  2:299484/1‑90 (MQ=255)
                gCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCggg                                                                      >  1:278709/1‑90 (MQ=255)
                 cTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGc                                                                     >  1:600051/1‑90 (MQ=255)
                     gCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAAc                                                                 >  1:22303/1‑90 (MQ=255)
                         cAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGAt                                                             <  2:22303/90‑1 (MQ=255)
                         cAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGAt                                                             <  1:412481/90‑1 (MQ=255)
                             gCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCcggc                                                                         >  1:393736/1‑74 (MQ=255)
                             gCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCcggc                                                                         <  2:393736/74‑1 (MQ=255)
                                     gCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCtt                                                 >  1:94091/1‑90 (MQ=255)
                                      ttgtatgaaaagcagatttaataccagctaactcaggttcatctccagcgCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtc                                                <  2:245990‑M2/40‑1 (MQ=255)
                                      cATGTTATTCCTGCCTGTATTGGTTTCGAGGATCCATGCCATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtc                                                <  2:523657/90‑1 (MQ=255)
                                      cATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtc                                                <  1:478596/90‑1 (MQ=255)
                                      cATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtc                                                <  2:278709/90‑1 (MQ=255)
                                      cATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtc                                                <  2:94091/90‑1 (MQ=255)
                                       aTGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtca                                               >  2:382484/1‑90 (MQ=255)
                                       aTGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtca                                               >  1:211280/1‑90 (MQ=255)
                                        tGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtcaa                                              <  2:669701/90‑1 (MQ=255)
                                        tGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtcaa                                              <  2:215677/90‑1 (MQ=255)
                                        gtatgaaaagcagatttaataccagctaactcaggttcatctccagcgCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTtcaa                                              <  2:636043‑M2/42‑1 (MQ=255)
                                                 tGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACa                                     <  2:266428/90‑1 (MQ=255)
                                                   ccTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAgc                                   <  2:293854/90‑1 (MQ=255)
                                                      gTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCa                                >  1:116572/1‑90 (MQ=255)
                                                      gTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCa                                >  2:734355/1‑90 (MQ=255)
                                                      gTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCa                                >  1:390943/1‑90 (MQ=255)
                                                       tATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTCCAATCAATAACAGCGCaa                               <  1:624658/90‑1 (MQ=255)
                                                           ggTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCAc                           >  1:468349/1‑90 (MQ=255)
                                                             tttCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCg                         >  1:17959/1‑90 (MQ=255)
                                                             tttCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCg                         >  1:425620/1‑90 (MQ=255)
                                                                         aTGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGc             <  1:382484/90‑1 (MQ=255)
                                                                            aCATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGa          <  2:425620/90‑1 (MQ=255)
                                                                            aCATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGa          <  2:17959/90‑1 (MQ=255)
                                                                               tctccagcgCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACt       >  2:762247‑M2/10‑90 (MQ=255)
                                                                                ctccagcgCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACTg      >  2:496137‑M2/9‑90 (MQ=255)
                                                                                  ccagcgCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACTGAt    >  2:152894‑M2/7‑90 (MQ=255)
                                                                                    agcgCATCGATCTTCCGGCGGGCAAATGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACTGATAc  <  2:22176‑M2/86‑1 (MQ=255)
                                                                                    agcgCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACTGATAc  >  2:306862‑M2/5‑90 (MQ=255)
                                                                                        |                                                                                     
CGGTAACGCTTGTACGGCTCCGCCCCAATGCGTTCCAGCATGTTATTCATGCCTGTATTGGTTTCGAGGATCCATGACATCTCCAGCGCATCGATCTTCCGGCGGGCAAACGGATCGCGTAAGGCTTCAATCAATAACAGCGCAATCACCGGGCCGATGCGGCTGAACTGATAC  >  NZ_CP009273/3122424‑3122597

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