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] 42 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.080) 27/92 0.1 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

TGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATGTCATTTCCCCTTCACTGTTGTTTCCGGGGAAATATTTGCTTCATCCC  >  NC_000913/3453881‑3453975
                                              |                                                
tGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtat                                               <  1:1178484‑M1/50‑4 (MQ=255)
tGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtat                                               <  1:2743068‑M1/50‑4 (MQ=255)
 gACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtat                                               >  1:1651134‑M1/1‑46 (MQ=255)
  aCCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtataa                                             <  1:1503593‑M1/50‑6 (MQ=255)
  aCCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATtata                                              <  1:2370851‑M1/49‑5 (MQ=255)
         tctttctGCTCACGCGCGGCGATATAGAATTCAAACATtataatccctt                                       <  1:658520‑M1/49‑12 (MQ=255)
         tctttctGCTCACGCGCGGCGATATAGAATTCAAACATtataatccctt                                       <  1:2757264‑M1/49‑12 (MQ=255)
             tctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtg                                    <  1:1956723‑M1/48‑15 (MQ=255)
              ctGCTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacg                                 >  1:1384914‑M1/1‑33 (MQ=255)
                 cTCACGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacgta                               >  1:1714476‑M1/1‑30 (MQ=255)
                    aCGCGCGGCGATATAGAATTCAAACATtataatcccttgtgacgtaaaaa                           >  1:1290987‑M1/1‑27 (MQ=255)
                      gcgcgGCGATATAGAATTCAAACATtataatcccttgtgacgtaaaaact                         <  1:287541‑M1/50‑26 (MQ=255)
                         cggcgATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgca                      <  1:2020252‑M1/50‑29 (MQ=255)
                           gcgATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaa                    <  1:2721239‑M1/50‑31 (MQ=255)
                            cgATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaa                   <  1:762551‑M1/50‑32 (MQ=255)
                             gATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaac                  <  1:2729104‑M1/50‑33 (MQ=255)
                             gATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaac                  <  1:1043106‑M1/50‑33 (MQ=255)
                             gATATAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaa                   >  1:1933247‑M1/1‑18 (MQ=255)
                              atatAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaac                  >  1:247213‑M1/1‑17 (MQ=255)
                              atatAGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaac                  >  1:1845159‑M1/1‑17 (MQ=255)
                                ataGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaa                >  1:2415455‑M1/1‑15 (MQ=255)
                                ataGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaa                >  1:2089625‑M1/1‑15 (MQ=255)
                                ataGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaa                >  1:1863523‑M1/1‑15 (MQ=255)
                                ataGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaa                >  1:1484600‑M1/1‑15 (MQ=255)
                                 taGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaa               <  1:2758679‑M1/49‑36 (MQ=255)
                                 taGAATTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaa               <  1:1381649‑M1/49‑36 (MQ=255)
                                     aTTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcc          >  1:25208‑M1/1‑10 (MQ=255)
                                     aTTCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcc          <  1:2749510‑M1/50‑41 (MQ=255)
                                      ttCAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcct         <  1:1850889‑M1/50‑42 (MQ=255)
                                        cAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctac       >  1:2187982‑M1/1‑7 (MQ=255)
                                        cAAACATtataatcccttgtgacgtaaaaactgcaaaacaaaagtccta        <  1:751795‑M1/49‑43 (MQ=255)
                                         aaaCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacg      >  1:1605397‑M1/1‑6 (MQ=255)
                                         aaaCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacg      >  1:1271616‑M1/1‑6 (MQ=255)
                                         aaaCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacg      >  1:1122453‑M1/1‑6 (MQ=255)
                                         aaaCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctac       <  1:2756237‑M1/49‑44 (MQ=255)
                                         aaaCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctac       <  1:2210167‑M1/49‑44 (MQ=255)
                                          aaCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgc     <  1:485863‑M1/50‑46 (MQ=255)
                                           aCATtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacg      >  1:284761‑M1/1‑4 (MQ=255)
                                             aTtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgcagg  <  1:1101413‑M1/50‑49 (MQ=255)
                                             aTtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgcagg  <  1:99235‑M1/50‑49 (MQ=255)
                                             aTtataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgcagg  <  1:1238337‑M1/50‑49 (MQ=255)
                                              ttataatcccttgtgacgtaaaaactgcaaaacaaaagtcctacgcagg  >  1:2143131‑M1/1‑1 (MQ=255)
                                              |                                                
TGACCTGTTTCTTTCTGCTCACGCGCGGCGATATAGAATTCAAACATGTCATTTCCCCTTCACTGTTGTTTCCGGGGAAATATTTGCTTCATCCC  >  NC_000913/3453881‑3453975

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: