Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 15,389 Δ5,175 bp insL‑1insA‑1 insL‑1, hokC, mokC, sokC, nhaA, nhaR, insB‑1, insA‑1

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 15389 20563 5175 50 [0] [0] 49 insL‑1–insA‑1 insL‑1,hokC,mokC,sokC,nhaA,nhaR,insB‑1,insA‑1

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 153880 (0.000)45 (0.900) 30/96 0.2 100% noncoding (2/1345 nt) IS186 repeat region
?NC_000913 20564 = 0 (0.000)intergenic (‑56/+251) insA‑1/rpsT IS1 protein InsA/30S ribosomal subunit protein S20

AGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTT  >  NC_000913/15340‑15435
                                                |                                               
aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTC                                                 >  1:2211895‑M1/1‑49 (MQ=255)
aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTC                                                 >  1:1746587‑M1/1‑49 (MQ=255)
aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAACACTC                                                 >  1:780940‑M1/1‑49 (MQ=255)
 gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCt                                                <  1:1953573‑M1/49‑2 (MQ=255)
      cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCtacccc                                           >  1:52573‑M1/1‑43 (MQ=255)
       gTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCcacctcct                                         >  1:1240364‑M1/1‑42 (MQ=255)
       gTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaa                                         <  1:3900528‑M1/50‑9 (MQ=255)
       gTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaa                                         <  1:2379865‑M1/50‑9 (MQ=255)
        ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaat                                        <  1:99072‑M1/50‑10 (MQ=255)
        ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaat                                        <  1:7367‑M1/50‑10 (MQ=255)
         tGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCtacctcaatg                                       >  1:1262740‑M1/1‑40 (MQ=255)
          gAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTTACTCttcctccatgt                                      >  1:1006709‑M1/1‑39 (MQ=255)
          gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatg                                       <  1:2179943‑M1/49‑11 (MQ=255)
           aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtg                                     <  1:3509947‑M1/50‑13 (MQ=255)
           aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtg                                     <  1:2707079‑M1/50‑13 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCcacctccatgtg                                     >  1:3640134‑M1/1‑37 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgt                                    <  1:1390436‑M1/50‑14 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgt                                    <  1:745226‑M1/50‑14 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgt                                    <  1:2739631‑M1/50‑14 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgt                                    <  1:1732126‑M1/50‑14 (MQ=255)
             gATATGCGAGCACCTGTAAAGTGGCGGGGGTCAACCtacctcaaagtgta                                   >  1:1436153‑M1/1‑36 (MQ=255)
              atatGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtat                                  <  1:897340‑M1/50‑16 (MQ=255)
               tatGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatc                                 <  1:685296‑M1/50‑17 (MQ=255)
               tatGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatc                                 <  1:1145161‑M1/50‑17 (MQ=255)
                 tGCGAGCACCTGTAAAGTGGCGGGGATTACACcacctcaatgtgtatta                                >  1:139144‑M1/1‑32 (MQ=255)
                 tGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatca                                <  1:1941870‑M1/49‑18 (MQ=255)
                   cGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaa                             <  1:3057689‑M1/50‑21 (MQ=255)
                   cGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcaca                              <  1:2205860‑M1/49‑20 (MQ=255)
                     aGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaata                           <  1:1429566‑M1/50‑23 (MQ=255)
                       cACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatc                         <  1:1629297‑M1/50‑25 (MQ=255)
                         ccTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatcca                       <  1:3243566‑M1/50‑27 (MQ=255)
                             tAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccata                     <  1:1199590‑M1/48‑29 (MQ=255)
                               aaGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccatattc                  <  1:2996193‑M1/49‑32 (MQ=255)
                               aaGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccatattc                  <  1:1373827‑M1/49‑32 (MQ=255)
                                 gTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccatattcttt               <  1:3145514‑M1/50‑35 (MQ=255)
                                   ggCGGGGATCACTCtacctcaatgtgtatcacaatatccatattctttgt             >  1:1101201‑M1/1‑14 (MQ=255)
                                    gCGGGGATCACTCtacctcaatgtgtatcacaatatccatattctttgtg            >  1:3494708‑M1/1‑13 (MQ=255)
                                      ggggATCACTCtacctcaatgtgtatcacattatccatattctttgtggg          >  1:1354537‑M1/1‑11 (MQ=255)
                                      ggggATCACTCtacctcaatgtgtatcacaatatccatattctttgtggg          >  1:1708503‑M1/1‑11 (MQ=255)
                                      ggggATCACTCtacctcaatgtgtatcacaatatccatattctttgtggg          >  1:2453659‑M1/1‑11 (MQ=255)
                                         gATCACTCtacctcaatgtgtatcacaatatccatattctttgtgggg         >  1:1112647‑M1/1‑8 (MQ=255)
                                           tCACTCtacctcaatgtgtatcacaatatccatattctttgtgggggagt     <  1:2975548‑M1/50‑45 (MQ=255)
                                            cACTCtacctcaatgtgtatcacaatatccatattctttgtgggggagtc    >  1:208463‑M1/1‑5 (MQ=255)
                                             aCTCtacctcaatgtgtatcacaatatccatattctttgtgggggagtc    <  1:745221‑M1/49‑46 (MQ=255)
                                              ctctacctcaatgtgtatcacaatatccatattatttgtgggggagtctg  <  1:3770063‑M1/50‑48 (MQ=255)
                                               tctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:35571‑M1/1‑2 (MQ=255)
                                               tctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:3502003‑M1/1‑2 (MQ=255)
                                               tctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:798557‑M1/1‑2 (MQ=255)
                                                ctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:1034914‑M1/1‑1 (MQ=255)
                                                ctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:1301412‑M1/1‑1 (MQ=255)
                                                |                                               
AGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTT  >  NC_000913/15340‑15435

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