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 49 [0] [0] 52 [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)47 (0.850)
+GA
39/94 0.0 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:644628‑M2/50‑50 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                                  <  1:1996824‑M2/2‑1 (MQ=255)
   ttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                                >  1:334426‑M2/47‑50 (MQ=255)
      tcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAt                                             >  1:690755‑M2/44‑50 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCa                                              <  1:1168051‑M2/6‑1 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtt                                            <  1:3614348‑M2/8‑1 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAt                                             <  1:3578971‑M2/7‑1 (MQ=255)
         gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                          <  1:788630‑M2/10‑1 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                         <  1:3113161‑M2/11‑1 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                         <  1:2698015‑M2/11‑1 (MQ=255)
           agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGa                                        <  1:1117214‑M2/12‑1 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                      <  1:1369427‑M2/14‑1 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                      >  1:3266070‑M2/37‑50 (MQ=255)
              actgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                     <  1:2393988‑M2/15‑1 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAttt                                    <  1:4073728‑M2/16‑1 (MQ=255)
                tgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                   >  1:495455‑M2/34‑50 (MQ=255)
                 gacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                                  >  1:887372‑M2/33‑50 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                                 <  1:2121550‑M2/19‑1 (MQ=255)
                   catcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                                 >  1:3114500‑M2/31‑49 (MQ=255)
                       aatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACgg                            <  1:1735542‑M2/24‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:2481287‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:541517‑M2/26‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                          <  1:295902‑M2/26‑1 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:2230135‑M2/23‑49 (MQ=255)
                           tgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                         >  1:4077135‑M2/23‑49 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGc                        <  1:3800215‑M2/28‑1 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCa                       >  1:2305121‑M2/22‑50 (MQ=255)
                            gtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCa                       >  1:3770932‑M2/22‑50 (MQ=255)
                             tgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                      >  1:589152‑M2/21‑50 (MQ=255)
                             tgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                      <  1:693985‑M2/30‑1 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTa                    <  1:1566721‑M2/32‑1 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                   <  1:1094931‑M2/33‑1 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATc                  <  1:4313433‑M2/34‑1 (MQ=255)
                                  cccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCa                 <  1:904484‑M2/35‑1 (MQ=255)
                                  cccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCa                 <  1:4145141‑M2/35‑1 (MQ=255)
                                   ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAg                <  1:1453240‑M2/36‑1 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt               <  1:3129610‑M2/37‑1 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt               >  1:2881401‑M2/14‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt          <  1:939376‑M2/42‑1 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt          >  1:1609205‑M2/8‑49 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt          >  1:566402‑M2/8‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt         >  1:647885‑M2/7‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa        <  1:2301576‑M2/44‑1 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa        <  1:4394357‑M2/44‑1 (MQ=255)
                                            atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa       >  1:1291053‑M2/6‑50 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc      <  1:2998215‑M2/46‑1 (MQ=255)
                                              agaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg     >  1:1111635‑M2/4‑50 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg    >  1:3561340‑M2/3‑50 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg    >  1:800113‑M2/3‑50 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACgg    >  1:1043705‑M2/3‑50 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg     <  1:3446491‑M2/47‑1 (MQ=255)
                                                 aCGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCa  <  1:13421/50‑1 (MQ=255)
                                                 |                                                 
CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCA  >  NC_000913/3136333‑3136431

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