Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,453,928 Δ15,926 bp gspA[chiA] 16 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3453928 3469853 15926 42 [0] [0] 41 gspA–[chiA] gspA,gspC,gspD,gspE,gspF,gspG,gspH,gspI,gspJ,gspK,gspL,gspM,gspO,bfr,bfd,[chiA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 34539270 (0.000)38 (1.330) 32/92 0.0 100% coding (1/420 nt) gspB putative general secretion pathway protein B
?NC_000913 3469854 = 0 (0.000)intergenic (‑1/+291) chiA/tufA endochitinase/translation elongation factor Tu 1

GGGTGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATGTCATTTCCCCTTCACTGTTGTTTCCGGGGAAATATTTGCTTCATCCC  >  NC_000913/3453878‑3453975
                                                 |                                                
gggTGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACAT                                                  <  1:2561042‑M1/50‑1 (MQ=255)
gggTGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACAT                                                  >  1:1856750‑M1/1‑50 (MQ=255)
  gTGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATta                                                >  1:348419‑M1/1‑48 (MQ=255)
   tGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtat                                               <  1:2037696‑M1/50‑4 (MQ=255)
    gACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtat                                               >  1:359023‑M1/1‑46 (MQ=255)
     aCCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtataa                                             <  1:2470732‑M1/50‑6 (MQ=255)
     aCCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtata                                              <  1:1035780‑M1/49‑5 (MQ=255)
      ccTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtataat                                            <  1:224651‑M1/50‑7 (MQ=255)
          tttctttcTGCTCACGCGCGGCGATATAGAATTCAAACATtataatccct                                        <  1:1424369‑M1/50‑11 (MQ=255)
          tttctttcTGCTCACGCGCGGCGATATAGAATTCAAACATtataatcc                                          <  1:1008142‑M1/48‑9 (MQ=255)
            tctttctGCTCACGCGCGGCGATATAGAATTCAAACATtataatccctt                                       <  1:192830‑M1/49‑12 (MQ=255)
             ctttctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgt                                     >  1:1136261‑M1/1‑37 (MQ=255)
              tttctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgt                                     >  1:1810072‑M1/1‑36 (MQ=255)
                tctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgac                                  >  1:1169340‑M1/1‑34 (MQ=255)
                tctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgac                                  >  1:1522825‑M1/1‑34 (MQ=255)
                tctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtg                                    <  1:1455139‑M1/48‑15 (MQ=255)
                 ctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacg                                 >  1:2283186‑M1/1‑33 (MQ=255)
                     tCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacgtaa                              >  1:1621373‑M1/1‑29 (MQ=255)
                      cACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacgtaaa                             >  1:425721‑M1/1‑28 (MQ=255)
                      cACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacgtaaa                             >  1:91798‑M1/1‑28 (MQ=255)
                          cgcgGCGATATAGAATTCAAACATtataatcccttgtgacgtaaaaactg                        >  1:1454945‑M1/1‑24 (MQ=255)
                           gcggcgATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgc                       <  1:2682069‑M1/50‑28 (MQ=255)
                            cggcgATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgc                       >  1:1513457‑M1/1‑22 (MQ=255)
                             ggcgATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaa                     <  1:274385‑M1/50‑30 (MQ=255)
                                gATATCGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaac                  <  1:1143559‑M1/50‑33 (MQ=255)
                                gATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaa                   >  1:921993‑M1/1‑18 (MQ=255)
                                 atatAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaaca                 <  1:462520‑M1/50‑34 (MQ=255)
                                   ataGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaa                >  1:215447‑M1/1‑15 (MQ=255)
                                   ataGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaa                >  1:2613327‑M1/1‑15 (MQ=255)
                                    taGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaa              <  1:2166755‑M1/50‑37 (MQ=255)
                                    taGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaa               <  1:93488‑M1/49‑36 (MQ=255)
                                     aGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaag             >  1:391930‑M1/1‑13 (MQ=255)
                                     aGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaag             >  1:1570638‑M1/1‑13 (MQ=255)
                                     aGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaag             >  1:1049277‑M1/1‑13 (MQ=255)
                                       aaTTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtc           >  1:1686225‑M1/1‑11 (MQ=255)
                                         ttCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcct         <  1:2080550‑M1/50‑42 (MQ=255)
                                         ttCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcct         <  1:111006‑M1/50‑42 (MQ=255)
                                          tCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtccta        >  1:243012‑M1/1‑8 (MQ=255)
                                           cAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctac       >  1:1320774‑M1/1‑7 (MQ=255)
                                              aCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgca    <  1:2442914‑M1/50‑47 (MQ=255)
                                                aTtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgcag   >  1:1353347‑M1/1‑2 (MQ=255)
                                                 ttataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgcagg  >  1:860493‑M1/1‑1 (MQ=255)
                                                 |                                                
GGGTGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATGTCATTTCCCCTTCACTGTTGTTTCCGGGGAAATATTTGCTTCATCCC  >  NC_000913/3453878‑3453975

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