New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 440028827 (0.470)34 (0.610) 26/94 0.5 56.3% coding (1/309 nt) hfq RNA‑binding protein Hfq
?NC_000913 4400296 = 27 (0.490)coding (9/309 nt) hfq RNA‑binding protein Hfq

AAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTG  >  NC_000913/4400247‑4400345
                                                 |                                                 
aaaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAg                                                   <  1:2489722/50‑1 (MQ=255)
 gaatcgaaaggttcaaagtacaaataagcatataaggaaaagagagaagg                                                  >  1:1856445‑M2/49‑50 (MQ=255)
  aGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAggg                                                 <  1:3069665/50‑1 (MQ=255)
    gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAgggg                                                >  1:2841033/1‑49 (MQ=255)
    gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGc                                               <  1:2560771/50‑1 (MQ=255)
     ttCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGc                                               >  1:729324/1‑49 (MQ=255)
     cgaaaggttcaaagtacaaataagcatataaggaaaagagagaaGGGGc                                               >  1:1095954‑M2/45‑49 (MQ=255)
        aaggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATc                                           >  1:4136221‑M2/42‑50 (MQ=255)
         aggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATCt                                          <  1:1011373‑M2/10‑1 (MQ=255)
         aggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATCt                                          <  1:3230738‑M2/10‑1 (MQ=255)
         aggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATCt                                          >  1:3092567‑M2/41‑50 (MQ=255)
         aggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATCt                                          >  1:113958‑M2/41‑50 (MQ=255)
         aggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATCt                                          >  1:1823739‑M2/41‑50 (MQ=255)
         aggttcaaagtacaaataagcatataaggaaaagagagaaGGGGCAATCt                                          >  1:1922270‑M2/41‑50 (MQ=255)
            tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTa                                       >  1:1055156/1‑50 (MQ=255)
              caaagtacaaataagcatataaggaaaagagagaaGGGGCAATCTTTAc                                      >  1:741979‑M2/36‑49 (MQ=255)
               aaagtacaaataagcatataaggaaaagagagaaGGGGCAATCTTTACa                                     <  1:640590‑M2/15‑1 (MQ=255)
                aagtacaaataagcatataaggaaaagagagaaGGGGCAATCTTTACAAg                                   <  1:2706528‑M2/17‑1 (MQ=255)
                aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTAc                                      >  1:3321423/1‑47 (MQ=255)
                aaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTAc                                      >  1:2797066/1‑47 (MQ=255)
                 agtacaaataagcatataaggaaaagagagaaGGGGCAATCTTTAc                                      >  1:1092003‑M2/33‑46 (MQ=255)
                  gtacaaataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGAt                                 <  1:421126‑M2/19‑1 (MQ=255)
                   tacaaataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATc                                <  1:281823‑M2/20‑1 (MQ=255)
                    acaaataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATcc                               >  1:1145575‑M2/30‑50 (MQ=255)
                    acaaataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATcc                               >  1:4002362‑M2/30‑50 (MQ=255)
                     gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCg                              >  1:2967756/1‑50 (MQ=255)
                       atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGtt                            >  1:1074065/1‑50 (MQ=255)
                       atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGtt                            >  1:952172/1‑50 (MQ=255)
                       aataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGtt                            <  1:2764118‑M2/24‑1 (MQ=255)
                       aataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGtt                            <  1:3297098‑M2/24‑1 (MQ=255)
                       aataagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGtt                            >  1:3567401‑M2/27‑50 (MQ=255)
                        tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTc                           >  1:2805927/1‑50 (MQ=255)
                         taagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTcc                          >  1:3520705‑M2/25‑50 (MQ=255)
                          taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCt                         <  1:1785898/50‑1 (MQ=255)
                          taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCt                         >  1:4120210/1‑50 (MQ=255)
                          aagcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCt                         >  1:3901616‑M2/24‑50 (MQ=255)
                           agcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTg                        <  1:309927‑M2/28‑1 (MQ=255)
                           aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCt                         >  1:2155894/1‑49 (MQ=255)
                            gcatataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGa                       <  1:3685817‑M2/29‑1 (MQ=255)
                             ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGaa                      >  1:28539/1‑50 (MQ=255)
                             catataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGaa                      >  1:1253868‑M2/21‑50 (MQ=255)
                             catataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGaa                      >  1:2558981‑M2/21‑50 (MQ=255)
                             catataaggaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGaa                      >  1:3496865‑M2/21‑50 (MQ=255)
                              gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAAc                     >  1:443112/1‑50 (MQ=255)
                                aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                    >  1:2264178/1‑49 (MQ=255)
                                aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                    >  1:3509884/1‑49 (MQ=255)
                                aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACg                    >  1:2934001/1‑49 (MQ=255)
                                     gaaaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGc              >  1:1092520‑M2/13‑50 (MQ=255)
                                      agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCg             >  1:2530431/1‑50 (MQ=255)
                                       gaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGt            <  1:2300924/50‑1 (MQ=255)
                                       aaagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGt            >  1:2464158‑M2/11‑50 (MQ=255)
                                        aagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTc           >  1:2888018‑M2/10‑50 (MQ=255)
                                        aagagagaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTc           <  1:2554831‑M2/41‑1 (MQ=255)
                                        aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTc           <  1:4119589/50‑1 (MQ=255)
                                         aTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCg          >  1:1082183/1‑50 (MQ=255)
                                          tGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCgg         >  1:2697021/1‑50 (MQ=255)
                                           agagaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCgg         >  1:959899‑M2/7‑49 (MQ=255)
                                             cTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGaa      >  1:2471256/1‑50 (MQ=255)
                                             cTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGaa      >  1:2384529/1‑50 (MQ=255)
                                             agaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGaa      >  1:615601‑M2/5‑50 (MQ=255)
                                              gaaGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAAc     <  1:3440068‑M2/47‑1 (MQ=255)
                                                aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAAc     <  1:2031445‑M2/47‑1 (MQ=255)
                                                aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGt   <  1:414222‑M2/49‑1 (MQ=255)
                                                aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGt   >  1:3234463‑M2/2‑50 (MQ=255)
                                                 ggggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACgtg  >  1:3624298/1‑50 (MQ=255)
                                                 |                                                 
AAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTG  >  NC_000913/4400247‑4400345

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