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 33 [0] [0] 32 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)31 (1.210) 23/90 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

CAGCCCCATACCAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAG  >  NC_000913/3469809‑3469900
                                             |                                              
acctgtttctttctgctcacgcgcggcgatatagaattcaaacattata                                             <  1:673024‑M2/4‑1 (MQ=255)
 cctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATaa                                            >  1:1015039‑M2/45‑49 (MQ=255)
  ctgtttctttctgctcacgcgcggcgatatagaattcaaacatTATAAt                                           <  1:445285‑M2/6‑1 (MQ=255)
     tttctttctgctcacgcgcggcgatatagaattcaaacatTATAATccc                                        <  1:810620‑M2/9‑1 (MQ=255)
           tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTtg                                     <  1:888001‑M2/12‑1 (MQ=255)
           tctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGa                                  >  1:2090975‑M2/35‑49 (MQ=255)
            ctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGa                                  <  1:1181671‑M2/15‑1 (MQ=255)
            ctgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGa                                  <  1:407854‑M2/15‑1 (MQ=255)
             tgctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACg                                >  1:1858793‑M2/33‑49 (MQ=255)
              gctcacgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACg                                >  1:459927‑M2/32‑48 (MQ=255)
                  acgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaaa                           <  1:1445025‑M2/22‑1 (MQ=255)
                   cgcgcggcgatatagaattcaaacatTATAATCCCTTGTGACGTaaaaa                          <  1:1186749‑M2/23‑1 (MQ=255)
                        ggcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGca                     <  1:1393587‑M2/28‑1 (MQ=255)
                         gcgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGca                     >  1:1090209‑M2/21‑48 (MQ=255)
                          cgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaa                   <  1:768359‑M2/30‑1 (MQ=255)
                          cgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaa                    >  1:1549159‑M2/20‑48 (MQ=255)
                          cgatatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaa                    >  1:473719‑M2/20‑48 (MQ=255)
                             tatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaaca                <  1:1904507‑M2/33‑1 (MQ=255)
                             tatagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaaca                <  1:1288589‑M2/33‑1 (MQ=255)
                              atagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaaca                >  1:864064‑M2/16‑48 (MQ=255)
                               tagaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa               >  1:1913549‑M2/15‑48 (MQ=255)
                                agaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa               >  1:448003‑M2/14‑47 (MQ=255)
                                agaattcaaacatTATAATCCCTTGTGACGTAAAAACTGcaaaacaa               >  1:546098‑M2/14‑47 (MQ=255)
                                    ttcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc         >  1:1362704‑M2/10‑49 (MQ=255)
                                     tcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc         >  1:284168‑M2/9‑48 (MQ=255)
                                     tcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc         >  1:984255‑M2/9‑48 (MQ=255)
                                     tcaaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTcc         >  1:1907915‑M2/9‑48 (MQ=255)
                                       aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc      <  1:727979‑M2/43‑1 (MQ=255)
                                       aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc      <  1:421292‑M2/43‑1 (MQ=255)
                                       aaacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTAc      <  1:1072493‑M2/43‑1 (MQ=255)
                                        aacatTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACg     <  1:2009409‑M2/44‑1 (MQ=255)
                                           atTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAg  <  1:89404‑M2/47‑1 (MQ=255)
                                             |                                              
CAGCCCCATACCAATCATAGATTTAGTAAATATATTTAATTTCATTATAATCCCTTGTGACGTAAAAACTGCAAAACAAAAGTCCTACGCAG  >  NC_000913/3469809‑3469900

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