New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)22 (0.840)
+CGAGAT
18/86 0.2 54.5% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 4070444 = 42 (1.400)intergenic (‑43/+71) yihR/yihS putative aldose 1‑epimerase YihR/sulfoquinovose isomerase

CGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACC  >  NC_000913/4070396‑4070493
                                                |                                                 
cagaatttccaggtcagacgggctgcaagttgcagaccgttacgagatg                                                   >  1:632363‑M2/49‑49 (MQ=255)
 gACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGa                                                 >  1:2039731/1‑50 (MQ=255)
 gACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGa                                                 <  1:2132775/50‑1 (MQ=255)
 gACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGa                                                 <  1:336539/50‑1 (MQ=255)
 agaatctccaggtcagacgggctgcaagttgcagaccgttacgagatgg                                                  <  1:1208150‑M2/2‑1 (MQ=255)
  gaatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGAc                                                >  1:2233346‑M2/47‑50 (MQ=255)
  gaatttccaggtcagacgggctgcaagttgcagaccgttacgagatGGAc                                                <  1:1261839‑M2/4‑1 (MQ=255)
  aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAc                                                >  1:2170482/1‑50 (MQ=255)
    atttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACgg                                              <  1:2227528‑M2/6‑1 (MQ=255)
    atttccaggtcagacgggctgcaagttgcagaccgttacgagatGGACgg                                              >  1:913277‑M2/45‑50 (MQ=255)
    aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                              >  1:649612/1‑50 (MQ=255)
    aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                              >  1:115830/1‑50 (MQ=255)
    aTACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                              >  1:2195535/1‑50 (MQ=255)
       tccaggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTc                                            <  1:1261561‑M2/8‑1 (MQ=255)
        aGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTccc                                          <  1:785320/50‑1 (MQ=255)
         caggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTcccc                                         >  1:1664135‑M2/40‑50 (MQ=255)
         caggtcagacgggctgcaagttgcagaccgttacgagatGGACGGTcccc                                         >  1:171193‑M2/40‑50 (MQ=255)
           cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTc                                       <  1:2026163/50‑1 (MQ=255)
           cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTc                                       <  1:1780708/50‑1 (MQ=255)
           cTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTc                                       >  1:45990/1‑50 (MQ=255)
             cccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                     >  1:540621/1‑50 (MQ=255)
              ccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                    >  1:1671507/1‑50 (MQ=255)
              ccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                    >  1:2056401/1‑50 (MQ=255)
              cagacgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGcc                                    >  1:940143‑M2/35‑50 (MQ=255)
                gTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcccc                                  >  1:1289718/1‑50 (MQ=255)
                  gAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                >  1:1542962/1‑50 (MQ=255)
                  cgggctgcaagttgcagaccgttacgagatGGACGGTCCCCTcgccccct                                <  1:348486‑M2/20‑3 (MQ=255)
                   gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                               <  1:81795‑M2/21‑1 (MQ=255)
                   gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt                               <  1:1352855‑M2/21‑1 (MQ=255)
                   aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt                               >  1:124924/1‑50 (MQ=255)
                   aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                <  1:16434/49‑1 (MQ=255)
                    ggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTg                              >  1:1398580‑M2/29‑50 (MQ=255)
                     gctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgg                             <  1:1025323‑M2/23‑1 (MQ=255)
                       tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:1804110/1‑49 (MQ=255)
                        gATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:1959212/1‑48 (MQ=255)
                        gATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:306879/1‑48 (MQ=255)
                         caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgggg                           <  1:145944‑M2/25‑1 (MQ=255)
                         aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:1111183/1‑47 (MQ=255)
                         aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                            >  1:951111/1‑47 (MQ=255)
                          aagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGga                          >  1:1528390‑M2/23‑48 (MQ=255)
                           agttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGagag                       <  1:1636427‑M2/29‑1 (MQ=255)
                            gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg                      <  1:462883‑M2/30‑1 (MQ=255)
                              aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                    >  1:1373806/1‑50 (MQ=255)
                               gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:1115953‑M2/18‑49 (MQ=255)
                               aaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGtt                   <  1:318898/50‑1 (MQ=255)
                               aaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                    >  1:158113/1‑49 (MQ=255)
                               aaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                    >  1:1683566/1‑49 (MQ=255)
                                cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGtt                   <  1:636809‑M2/33‑1 (MQ=255)
                                cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt                    >  1:2130177‑M2/17‑48 (MQ=255)
                                  ccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                <  1:92094/50‑1 (MQ=255)
                                   cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg               <  1:530260/50‑1 (MQ=255)
                                   cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg               <  1:604543/50‑1 (MQ=255)
                                   cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg               <  1:1721008/50‑1 (MQ=255)
                                   cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                >  1:1368552/1‑49 (MQ=255)
                                   cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                >  1:413833/1‑49 (MQ=255)
                                    aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                >  1:950213/1‑48 (MQ=255)
                                    aaaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                >  1:1028823/1‑48 (MQ=255)
                                      aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga            <  1:1181569/50‑1 (MQ=255)
                                        aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg          >  1:1163452/1‑50 (MQ=255)
                                            cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa      <  1:1855386/50‑1 (MQ=255)
                                            cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa       >  1:1392417/1‑49 (MQ=255)
                                              ccGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa    <  1:819394/50‑1 (MQ=255)
                                               cGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc   >  1:2179333/1‑50 (MQ=255)
                                               cGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc   <  1:96078/50‑1 (MQ=255)
                                                ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc  >  1:1556195/1‑50 (MQ=255)
                                                |                                                 
CGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACC  >  NC_000913/4070396‑4070493

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