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 96 [0] [0] 95 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)85 (1.290) 36/62 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

CAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAA  >  NC_000913/3469820‑3469887
                                  |                                 
tctgctcacgcgcggcgatatagaattcaaacatt                                   >  1:3954604‑M2/35‑35 (MQ=255)
tctgctcacgcgcggcgatatagaattcaaacatt                                   <  1:5840709‑M2/1‑1 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  >  1:5649795‑M2/34‑35 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  >  1:7513953‑M2/34‑35 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  >  1:7143917‑M2/34‑35 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  >  1:4705079‑M2/34‑35 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  >  1:3998105‑M2/34‑35 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  <  1:4012385‑M2/2‑1 (MQ=255)
 ctgctcacgcgcggcgatatagaattcaaacatta                                  >  1:6551740‑M2/34‑35 (MQ=255)
  tgctcacgcgcggcgatatagaattcaaacattat                                 <  1:2197601‑M2/3‑1 (MQ=255)
  tgctcacgcgcggcgatatagaattcaaacattat                                 <  1:5328176‑M2/3‑1 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:6170636‑M2/31‑35 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:2868845‑M2/31‑35 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:278363‑M2/31‑35 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:1188964‑M2/31‑35 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:1363526‑M2/31‑35 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:4757167‑M2/31‑35 (MQ=255)
    ctcacgcgcggcgatatagaattcaaacatTATaa                               >  1:2683767‑M2/31‑35 (MQ=255)
     tcacgcgcggcgatatagaattcaaacatTATAAt                              <  1:6785501‑M2/6‑1 (MQ=255)
     tcacgcgcggcgatatagaattcaaacatTATAAt                              <  1:2551760‑M2/6‑1 (MQ=255)
      cacgcgcggcgatatagaattcaaacatTATAATc                             <  1:3114348‑M2/7‑1 (MQ=255)
       acgcgcggcgatatagaattcaaacatTATAATcc                            >  1:5488597‑M2/28‑35 (MQ=255)
       acgcgcggcgatatagaattcaaacatTATAATcc                            <  1:1234320‑M2/8‑1 (MQ=255)
        cgcgcggcgatatagaattcaaacatTATAATccc                           >  1:4508930‑M2/27‑35 (MQ=255)
        cgcgcggcgatatagaattcaaacatTATAATccc                           >  1:3009794‑M2/27‑35 (MQ=255)
           gcggcgatatagaattcaaacatTATAATCCCTtg                        <  1:6001417‑M2/12‑1 (MQ=255)
           gcggcgatatagaactcaaacatTATAATCCCTtg                        <  1:6960874‑M2/12‑1 (MQ=255)
            cggcgatatagaattcaaacatTATAATCCCTtgt                       <  1:3839433‑M2/13‑1 (MQ=255)
            cggcgatatagaattcaaacatTATAATCCCTtgt                       <  1:5653803‑M2/13‑1 (MQ=255)
            cggcgatatagaattcaaacatTATAATCCCTtgt                       <  1:591822‑M2/13‑1 (MQ=255)
            cggcgatatagaattcaaacatTATAATCCCTtgt                       <  1:6966301‑M2/13‑1 (MQ=255)
            cggcgatatagaattcaaacatTATAATCCCTtgt                       <  1:1163993‑M2/13‑1 (MQ=255)
              gcgatatagaattcaaacatTATAATCCCTTGTGa                     <  1:4790686‑M2/15‑1 (MQ=255)
              gcgatatagaattcaaacatTATAATCCCTTGTGa                     <  1:2330935‑M2/15‑1 (MQ=255)
               cgatatagaattcaaacatTATAATCCCTTGTGAc                    <  1:1570244‑M2/16‑1 (MQ=255)
                 atatagaattcaaacatTATAATCCCTTGTGACGt                  >  1:5471168‑M2/18‑35 (MQ=255)
                 atatagaattcaaacatTATAATCCCTTGTGACGt                  >  1:4064958‑M2/18‑35 (MQ=255)
                 atatagaattcaaacatTATAATCCCTTGTGACGt                  >  1:3714587‑M2/18‑35 (MQ=255)
                 atagagaattcaaacatTATAATCCCTTGTGACGt                  >  1:807197‑M2/18‑35 (MQ=255)
                  tatagaattcaaacatTATAATCCCTTGTGACGTa                 >  1:755357‑M2/17‑35 (MQ=255)
                  tatagaattcaaacatTATAATCCCTTGTGACGTa                 <  1:3631321‑M2/19‑1 (MQ=255)
                  tatagaattcaaacatTATAATCCCTTGTGACGTa                 <  1:6978479‑M2/19‑1 (MQ=255)
                  tatagaattcaaacatTATAATCCCTTGTGACGTa                 >  1:2012872‑M2/17‑35 (MQ=255)
                  tatagaattcaaacatTATAATCCCTTGTGACGTa                 <  1:3378968‑M2/19‑1 (MQ=255)
                  tatagaattcaaacatTACAATCCCTTGTGACGTa                 >  1:2035151‑M2/17‑35 (MQ=255)
                   atagaattcaaacatTATAATCCCTTGTGACGTaa                >  1:604922‑M2/16‑35 (MQ=255)
                   atagaattcaaacatTATAATCCCTTGTGACGTaa                <  1:3787534‑M2/20‑1 (MQ=255)
                   atagaattcaaacatTATAATCCCTTGTGACGTaa                >  1:6352111‑M2/16‑35 (MQ=255)
                   atagaattcaaacatTATAATCCCTTGTGACGTaa                >  1:3248187‑M2/16‑35 (MQ=255)
                   atagaattcaaacatTATAATCCCTTGTGACGTaa                >  1:4684303‑M2/16‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               <  1:3239699‑M2/21‑1 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               <  1:3801406‑M2/21‑1 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               >  1:6849675‑M2/15‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               <  1:6973579‑M2/21‑1 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               >  1:5739891‑M2/15‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               >  1:6932389‑M2/15‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               >  1:6095796‑M2/15‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               <  1:3239729‑M2/21‑1 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               >  1:7409109‑M2/15‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               >  1:7530627‑M2/15‑35 (MQ=255)
                    tagaattcaaacatTATAATCCCTTGTGACGTaaa               <  1:3683124‑M2/21‑1 (MQ=255)
                     agaattcaaacatTATAATCCCTTGTGACGTaaaa              >  1:431388‑M2/14‑35 (MQ=255)
                      gaattcaaacatTATAATCCCTTGTGACGTaaaaa             >  1:2610937‑M2/13‑35 (MQ=255)
                      gaattcaaacatTATAATCCCTTGTGACGTaaaaa             >  1:2501744‑M2/13‑35 (MQ=255)
                      gaattcaaacatTATAATCCCTTGTGACGTaaaaa             >  1:4353352‑M2/13‑35 (MQ=255)
                       aattcaaacatTATAATCCCTTTTGACGTAAAAAc            >  1:7103496‑M2/12‑35 (MQ=255)
                       aattcaaacatTATAATCCCTTGTGACGTAAAAAc            >  1:1926657‑M2/12‑35 (MQ=255)
                       aattcaaacatTATAATCCCTTGTGACGTAAAAAc            <  1:5445490‑M2/24‑1 (MQ=255)
                        attcaaacatTATAATCCCTTGTGACGTAAAAACt           <  1:6005034‑M2/25‑1 (MQ=255)
                        attcaaacatTATAATCCCTTGTGACGTAAAAACt           >  1:6121451‑M2/11‑35 (MQ=255)
                        attcaaacatTATAATCCCTTGTGACGTAAAAACt           >  1:2260282‑M2/11‑35 (MQ=255)
                        attcaaacatTATAATCCCTTGTGACGTAAAAACt           >  1:373055‑M2/11‑35 (MQ=255)
                         ttcaaacatTATAATCCCTTGTGACGTAAAAACTg          >  1:6399593‑M2/10‑35 (MQ=255)
                          tcaaacatTATAATCCCTTGTGACGTAAAAACTGc         <  1:217600‑M2/27‑1 (MQ=255)
                           caaacatTATAATCCCTTGTGACGTAAAAACTGCa        <  1:1061959‑M2/28‑1 (MQ=255)
                           caaacatTATAATCCCTTGTGACGTAAAAACTGCa        >  1:1290597‑M2/8‑35 (MQ=255)
                           caaacatTATAATCCCTTGTGACGTAAAAACTGCa        >  1:7344429‑M2/8‑35 (MQ=255)
                           caaacatTATAATCCCTTGTGACGTAAAAACTGCa        >  1:6339742‑M2/8‑35 (MQ=255)
                           caaacatTATAATCCCTTGTGACGTAAAAACTGCa        >  1:2244684‑M2/8‑35 (MQ=255)
                           caaacatTATAATCCCTTGTGACGTAAAAACTGCa        <  1:4931231‑M2/28‑1 (MQ=255)
                            aaacatTATAATCCCTTGTGACGTAAAAACTGCaa       <  1:3759121‑M2/29‑1 (MQ=255)
                            aaacatTATAATCCCTTGTGACGTAAAAACTGCaa       >  1:7049721‑M2/7‑35 (MQ=255)
                             aacatTATAATCCCTTGTGACGTAAAAACTGCaaa      >  1:624351‑M2/6‑35 (MQ=255)
                             aacatTATAATCCCTTGTGACGTAAAAACTGCaaa      >  1:1752853‑M2/6‑35 (MQ=255)
                              acatTATAATCCCTTGTGACGTAAAAACTGCaaaa     >  1:1622844‑M2/5‑35 (MQ=255)
                               catTATAATCCCTTGTGACGTAAAAACTGCAAAAc    >  1:2755588‑M2/4‑35 (MQ=255)
                               catTATAATCCCTTGTGACGTAAAAACTGCAAAAc    >  1:1258812‑M2/4‑35 (MQ=255)
                                atTATAATCCCTTGTGACGTAAAAACTGCAAAACa   <  1:4371157‑M2/33‑1 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  <  1:3837952‑M2/34‑1 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  >  1:4135228‑M2/2‑35 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  <  1:3435250‑M2/34‑1 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  >  1:7175969‑M2/2‑35 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  >  1:4719355‑M2/2‑35 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  >  1:2393522‑M2/2‑35 (MQ=255)
                                 tTATAATCCCTTGTGACGTAAAAACTGCAAAACaa  <  1:6301963‑M2/34‑1 (MQ=255)
                                  |                                 
CAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAA  >  NC_000913/3469820‑3469887

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