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 58 [0] [0] 53 [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)52 (0.560)
+GA
40/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

CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACG  >  NC_000913/3136333‑3136428
                                                 |                                              
ggcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaa                                                >  1:7170707‑M2/50‑50 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                               <  1:4169679‑M2/2‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                               <  1:2274679‑M2/2‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                               <  1:6350781‑M2/2‑1 (MQ=255)
 gcttctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaAc                                               <  1:6567206‑M2/2‑1 (MQ=255)
    tctcagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGc                                             >  1:5052981‑M2/46‑49 (MQ=255)
       cagcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtt                                         >  1:177479‑M2/43‑50 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAttt                                        >  1:5511504‑M2/42‑50 (MQ=255)
        agcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAt                                          <  1:4847144‑M2/7‑1 (MQ=255)
         gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                       >  1:6173640‑M2/41‑50 (MQ=255)
         gcagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                       <  1:4580590‑M2/10‑1 (MQ=255)
          cagcactgacatcaatttgtgtcacccgcagcgcatagaACGCCAtttt                                       >  1:2109856‑M2/40‑49 (MQ=255)
           agcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTg                                      >  1:4180869‑M2/39‑49 (MQ=255)
            gcactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGaa                                    <  1:6362382‑M2/13‑1 (MQ=255)
             cactgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAt                                   >  1:5820561‑M2/37‑50 (MQ=255)
              actgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                  >  1:4928561‑M2/36‑50 (MQ=255)
              actgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                  >  1:4227211‑M2/36‑50 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                  >  1:798434‑M2/35‑49 (MQ=255)
               ctgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAAtt                                  >  1:2200317‑M2/35‑49 (MQ=255)
                tgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                >  1:4540480‑M2/34‑50 (MQ=255)
                tgacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTa                                >  1:6785073‑M2/34‑50 (MQ=255)
                 gacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                               >  1:1981548‑M2/33‑50 (MQ=255)
                 gacatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTAc                               >  1:1664367‑M2/33‑50 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                              <  1:1525403‑M2/19‑1 (MQ=255)
                  acatcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACg                              <  1:3101157‑M2/19‑1 (MQ=255)
                     tcaatttgtgtctcccgcagcgcatagaACGCCATTTTGAATTTACGAAc                           >  1:6444694‑M2/29‑50 (MQ=255)
                     tcaatttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAAc                           <  1:2968449‑M2/22‑1 (MQ=255)
                         tttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGAc                       <  1:2707646‑M2/26‑1 (MQ=255)
                          ttgtgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACg                      >  1:7410459‑M2/24‑50 (MQ=255)
                             tgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                   <  1:7150999‑M2/30‑1 (MQ=255)
                             tgtcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAt                   <  1:2145221‑M2/30‑1 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCAtt                  <  1:6259342‑M2/31‑1 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTa                 <  1:5134136‑M2/32‑1 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTa                 <  1:6869232‑M2/32‑1 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTa                 >  1:6525512‑M2/19‑50 (MQ=255)
                               tcacccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTa                 >  1:7308394‑M2/19‑50 (MQ=255)
                                 acccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTAt                <  1:5899367‑M2/33‑1 (MQ=255)
                                   ccgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAg             >  1:205763‑M2/15‑50 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt            <  1:1256468‑M2/37‑1 (MQ=255)
                                    cgcagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGt            >  1:3240581‑M2/14‑50 (MQ=255)
                                      cagcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGa          <  1:6239733‑M2/39‑1 (MQ=255)
                                       agcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGAc         >  1:3228434‑M2/11‑50 (MQ=255)
                                        gcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACt        <  1:6840485‑M2/41‑1 (MQ=255)
                                        gcgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACt        >  1:3321882‑M2/10‑50 (MQ=255)
                                         cgcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt       <  1:6313132‑M2/42‑1 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt      <  1:46166‑M2/43‑1 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt       >  1:5247135‑M2/8‑49 (MQ=255)
                                          gcatagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACtt       >  1:7223396‑M2/8‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACttt      >  1:2487175‑M2/7‑49 (MQ=255)
                                           catagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa     <  1:7437890‑M2/44‑1 (MQ=255)
                                            atagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTaa    >  1:4245033‑M2/6‑50 (MQ=255)
                                             tagaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAAc   <  1:6816933‑M2/46‑1 (MQ=255)
                                               gaACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACg  <  1:3657732‑M2/47‑1 (MQ=255)
                                                 |                                              
CCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACG  >  NC_000913/3136333‑3136428

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