Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,110,680 25702 bp→GA [yghD]pitB 20 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3110680 3136381 25702 52 [0] [0] 55 [yghD]–pitB [yghD],yghE,yghF,yghG,pppA,yghJ,glcA,glcB,glcG,glcF,glcE,glcD,glcC,yghO,insH‑9,yghQ,yghR,yghS,yghT,pitB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 31106790 (0.000)49 (0.580)
+GA
42/94 0.2 100% coding (448/537 nt) yghD putative type II secretion system M‑type protein
?NC_000913 3136382 = 0 (0.000)intergenic (‑11/+281) pitB/gss metal phosphate:H(+) symporter PitB/fused glutathionylspermidine amidase/glutathionylspermidine synthetase

CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCA  >  NC_000913/3136333‑3136431
                                                 |                                                 
ggcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaa                                                   <  1:5241138‑M2/1‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                  <  1:3239510‑M2/2‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                  >  1:3776501‑M2/49‑50 (MQ=255)
   ttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                                >  1:2112275‑M2/47‑50 (MQ=255)
    tctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                                >  1:5253591‑M2/46‑49 (MQ=255)
     ctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGcc                                               <  1:5281808‑M2/5‑1 (MQ=255)
     ctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGcc                                               <  1:3184393‑M2/5‑1 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCa                                              <  1:3329004‑M2/6‑1 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAttt                                           <  1:1825065‑M2/9‑1 (MQ=255)
         gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                          >  1:3744278‑M2/41‑50 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                          >  1:3178385‑M2/40‑49 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                          >  1:4481313‑M2/40‑49 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                         <  1:6255694‑M2/11‑1 (MQ=255)
           agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGa                                        <  1:1474529‑M2/12‑1 (MQ=255)
            gcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                       >  1:4702539‑M2/38‑50 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                      <  1:6560574‑M2/14‑1 (MQ=255)
              actgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                     <  1:4203854‑M2/15‑1 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAttt                                    <  1:2252139‑M2/16‑1 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAttt                                    <  1:4969437‑M2/16‑1 (MQ=255)
                tgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                   >  1:4127072‑M2/34‑50 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                                  >  1:2679847‑M2/32‑49 (MQ=255)
                     tcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAAc                              <  1:1384657‑M2/22‑1 (MQ=255)
                      caatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                             >  1:449472‑M2/28‑50 (MQ=255)
                       aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                             <  1:1081037‑M2/23‑1 (MQ=255)
                       aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                             <  1:1814944‑M2/23‑1 (MQ=255)
                        atttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGa                           <  1:2191397‑M2/25‑1 (MQ=255)
                        atttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGa                           <  1:4485479‑M2/25‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          >  1:3172111‑M2/25‑50 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:279643‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:5467919‑M2/26‑1 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:6114814‑M2/23‑49 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         <  1:2344601‑M2/27‑1 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGc                        <  1:3434073‑M2/28‑1 (MQ=255)
                             tgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                      >  1:885888‑M2/21‑50 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAtt                     <  1:1007751‑M2/31‑1 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                   >  1:184677‑M2/17‑49 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                   >  1:4286606‑M2/17‑49 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATc                  <  1:3255205‑M2/34‑1 (MQ=255)
                                  cccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                   >  1:3365208‑M2/16‑48 (MQ=255)
                                   ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAg                >  1:2274832‑M2/15‑50 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt               <  1:521312‑M2/37‑1 (MQ=255)
                                     gcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt               >  1:955586‑M2/13‑49 (MQ=255)
                                       agcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGAc            <  1:1260173‑M2/40‑1 (MQ=255)
                                       agcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGAc            >  1:418902‑M2/11‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt          >  1:537897‑M2/9‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt          <  1:930040‑M2/42‑1 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt         >  1:4887275‑M2/7‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa        <  1:429507‑M2/44‑1 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc      <  1:260194‑M2/46‑1 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg    >  1:2099685‑M2/3‑50 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg     <  1:2132872‑M2/47‑1 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg     <  1:5446874‑M2/47‑1 (MQ=255)
                                                aACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg    >  1:249349‑M2/2‑49 (MQ=255)
                                                aACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGc   <  1:5451621‑M2/49‑1 (MQ=255)
                                                 aCGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCa  >  1:5513053/1‑50 (MQ=255)
                                                 |                                                 
CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCA  >  NC_000913/3136333‑3136431

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