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 41 [0] [0] 41 [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)38 (0.680)
+GA
29/94 0.3 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:2565086‑M2/1‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                  <  1:2822800‑M2/2‑1 (MQ=255)
         gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                          <  1:1893380‑M2/10‑1 (MQ=255)
         gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                          <  1:2748513‑M2/10‑1 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                      <  1:1647158‑M2/14‑1 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                      >  1:2686636‑M2/37‑50 (MQ=255)
              actgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                     >  1:3827363‑M2/36‑50 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                   <  1:1124081‑M2/17‑1 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                   <  1:3472856‑M2/17‑1 (MQ=255)
                   catcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGa                                <  1:4029199‑M2/20‑1 (MQ=255)
                    atcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGaa                               <  1:4175392‑M2/21‑1 (MQ=255)
                    atcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGaa                               >  1:3886938‑M2/30‑50 (MQ=255)
                     tcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAAc                              <  1:3378566‑M2/22‑1 (MQ=255)
                      caatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                             >  1:2490308‑M2/28‑50 (MQ=255)
                       aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACg                             <  1:234455‑M2/23‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:487854‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:3353456‑M2/26‑1 (MQ=255)
                          ttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:715701‑M2/24‑50 (MQ=255)
                          ttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:3984056‑M2/24‑50 (MQ=255)
                          ttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:3266710‑M2/24‑50 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:624514‑M2/23‑49 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGc                        <  1:584120‑M2/28‑1 (MQ=255)
                              gtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                      >  1:184192‑M2/20‑49 (MQ=255)
                              gtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                      >  1:2769856‑M2/20‑49 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTa                    >  1:4327002‑M2/19‑50 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                   >  1:2657384‑M2/17‑49 (MQ=255)
                                  cccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                   >  1:2162596‑M2/16‑48 (MQ=255)
                                   ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAg                <  1:3045923‑M2/36‑1 (MQ=255)
                                   ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAg                <  1:1123644‑M2/36‑1 (MQ=255)
                                     gcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTg              <  1:2545754‑M2/38‑1 (MQ=255)
                                     gcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTg              <  1:428811‑M2/38‑1 (MQ=255)
                                       agcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGAc            >  1:2835615‑M2/11‑50 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt          >  1:3024799‑M2/8‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa        <  1:2077218‑M2/44‑1 (MQ=255)
                                            atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa       >  1:3321635‑M2/6‑50 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc      >  1:115987‑M2/5‑50 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc      <  1:49608‑M2/46‑1 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc      <  1:3126706‑M2/46‑1 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg     <  1:1827199‑M2/47‑1 (MQ=255)
                                                aACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg    >  1:4268855‑M2/2‑49 (MQ=255)
                                                 aCGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCa  >  1:3005365/1‑50 (MQ=255)
                                                 |                                                 
CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCA  >  NC_000913/3136333‑3136431

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