New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3801110 (0.000)11 (0.410)
+CGAGAT
10/84 0.9 40.5% coding (732/774 nt) yaiO outer membrane protein YaiO
?NC_000913 4070444 = 37 (1.220)intergenic (‑43/+71) yihR/yihS putative aldose 1‑epimerase YihR/sulfoquinovose isomerase

AAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTC  >  NC_000913/380063‑380153
                                                |                                          
aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                            >  1:1880628/1‑49 (MQ=255)
 aGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAc                                           >  1:512182‑M1/1‑48 (MQ=255)
    aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgag                                        <  1:2331555‑M1/49‑5 (MQ=255)
         gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgga                                   >  1:202778‑M1/1‑40 (MQ=255)
           aTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacg                                 <  1:1028474‑M1/49‑12 (MQ=255)
                cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtccc                            >  1:1947580‑M1/1‑33 (MQ=255)
                cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtccc                            <  1:437464‑M1/49‑17 (MQ=255)
                    tCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcg                        >  1:2251460‑M1/1‑29 (MQ=255)
                           ggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt                 >  1:305736‑M1/1‑22 (MQ=255)
                           ggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt                 >  1:752089‑M1/1‑22 (MQ=255)
                             cTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgg               >  1:444940‑M1/1‑20 (MQ=255)
                                 aaGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg              >  1:395756‑M1/1‑16 (MQ=255)
                                        aGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtta    <  1:7476‑M1/49‑41 (MQ=255)
                                          aCCGTTAcgagatggacggtcccctcgcccccttggggagagggttagg  >  1:646038‑M1/1‑7 (MQ=255)
                                                |                                          
AAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTC  >  NC_000913/380063‑380153

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