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 30 [0] [0] 29 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)27 (0.740) 20/96 0.4 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:2315895‑M1/1‑49 (MQ=255)
 gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTTACTCtc                                               >  1:325826‑M1/1‑48 (MQ=255)
   gTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAACACTCtacc                                             >  1:2639477‑M1/1‑46 (MQ=255)
        ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaat                                        <  1:2530568‑M1/50‑10 (MQ=255)
          gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatg                                       <  1:2405175‑M1/49‑11 (MQ=255)
           aaGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCtacccccatgtg                                     >  1:955705‑M1/1‑38 (MQ=255)
           aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtg                                     <  1:2193639‑M1/50‑13 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCtacctccatgtg                                     >  1:77386‑M1/1‑37 (MQ=255)
            aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgt                                    <  1:1020065‑M1/50‑14 (MQ=255)
              atatGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtat                                  <  1:438012‑M1/50‑16 (MQ=255)
               tatGCGAGCACCTGTAAAGTGGCGGGGGTCACTCtacctccatgtgtat                                  >  1:975325‑M1/1‑34 (MQ=255)
               tatGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatc                                 <  1:130658‑M1/50‑17 (MQ=255)
               tatGCGAGCACCTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatc                                 <  1:2182685‑M1/50‑17 (MQ=255)
               tatGCGAGCACCTGTAAAGTGGCGGGGATCAATCtacctcaatgtgttt                                  >  1:1787643‑M1/1‑34 (MQ=255)
                atGCGAGCACCTGTAAAGTGGCGGGGGTCACTCtaccaccatgtggatga                                >  1:2495576‑M1/1‑33 (MQ=255)
                          cTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccat                      <  1:2411518‑M1/50‑28 (MQ=255)
                          cTGTAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccat                      <  1:410015‑M1/50‑28 (MQ=255)
                             tAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccatatt                   <  1:664547‑M1/50‑31 (MQ=255)
                             tAAAGTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccata                     <  1:1113063‑M1/48‑29 (MQ=255)
                                 gTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccatattcttt               <  1:1769874‑M1/50‑35 (MQ=255)
                                 gTGGCGGGGATCACTCtacctcaatgtgtatcacaatatccatattcttt               <  1:1281895‑M1/50‑35 (MQ=255)
                                      ggggATCACTCtacctcaatgtgtatcacaatatccatattctttgtggg          >  1:1496529‑M1/1‑11 (MQ=255)
                                      ggggATCACTCtacctcaatgtgtatcacaatatccatattctttgtggg          >  1:1012411‑M1/1‑11 (MQ=255)
                                       gggATCACTCtacctcaatgtgtatcacaatatccatattctttgtgggg         <  1:1821160‑M1/50‑41 (MQ=255)
                                         gATCACTCtacctcaatgtgtatcacaatatccatattctttgtggggga       <  1:91238‑M1/50‑43 (MQ=255)
                                             aCTCtgcctcaatgtgtatcacaatatccatattctttgtgggggagtc    <  1:687115‑M1/49‑46 (MQ=255)
                                             aCTCtacctcaatgtgtatcacaatatccatattctttgtgggggagtct   >  1:2474991‑M1/1‑4 (MQ=255)
                                               tctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:124388‑M1/1‑2 (MQ=255)
                                                ctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:2018981‑M1/1‑1 (MQ=255)
                                                ctacctcaatgtgtatcacaatatccatattctttgtgggggagtctg  >  1:1560701‑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 < 16 ≤ ATCG/ATCG < 32 ≤ 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.