Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,509,631 Δ8,319 bp [yfeO]pdeA [yfeO], ypeC, mntH, nupC, insL‑3, pdeA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2509631 2517949 8319 60 [0] [0] 59 [yfeO]–pdeA [yfeO],ypeC,mntH,nupC,insL‑3,pdeA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 25096300 (0.000)58 (1.000) 48/98 0.0 100% coding (1/1257 nt) yfeO putative transport protein YfeO
?NC_000913 2517950 = 0 (0.000)intergenic (‑118/+91) pdeA/alaX putative c‑di‑GMP phosphodiesterase PdeA/tRNA‑Ala

AACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAG  >  NC_000913/2517902‑2517998
                                                |                                                
tgcaagctgggcaaactaagtatctgaccccgcataaggaatagaacac                                                  <  1:326032‑M2/1‑1 (MQ=255)
  caagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAgg                                               <  1:2853993‑M2/4‑1 (MQ=255)
  caagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAgg                                               >  1:281793‑M2/47‑50 (MQ=255)
  caagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAg                                                <  1:3046339‑M2/3‑1 (MQ=255)
   aagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggc                                              >  1:303098‑M2/46‑50 (MQ=255)
   aagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggc                                              >  1:4390861‑M2/46‑50 (MQ=255)
   aagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggc                                              <  1:2495275‑M2/5‑1 (MQ=255)
    agctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcg                                             >  1:290723‑M2/45‑50 (MQ=255)
    agctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcg                                             <  1:1217605‑M2/6‑1 (MQ=255)
    agctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcg                                             >  1:2008165‑M2/45‑50 (MQ=255)
     gctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcgg                                            >  1:3931272‑M2/44‑50 (MQ=255)
      ctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcgg                                            >  1:2700900‑M2/43‑49 (MQ=255)
      ctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcgg                                            >  1:2445841‑M2/43‑49 (MQ=255)
       tgggcaaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCa                                          >  1:3591913‑M2/42‑50 (MQ=255)
         ggcaaactaagtatctgaccccgcataaggaatagcacaCAGGCGGCACt                                        <  1:3900734‑M2/11‑1 (MQ=255)
         ggcaaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCAc                                         >  1:3092484‑M2/40‑49 (MQ=255)
         ggcaaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACt                                        <  1:3608705‑M2/11‑1 (MQ=255)
          gcaaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTa                                       <  1:1861509‑M2/12‑1 (MQ=255)
           caaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTaa                                      <  1:3807261‑M2/13‑1 (MQ=255)
           caaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTaa                                      >  1:3611155‑M2/38‑50 (MQ=255)
           caaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTaa                                      >  1:4139778‑M2/38‑50 (MQ=255)
            aaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAc                                     >  1:116401‑M2/37‑50 (MQ=255)
            aaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAc                                     <  1:2521432‑M2/14‑1 (MQ=255)
             aactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAcc                                    <  1:209471‑M2/15‑1 (MQ=255)
             aactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAcc                                    >  1:1327296‑M2/36‑50 (MQ=255)
                  agtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGAc                               <  1:1171095‑M2/20‑1 (MQ=255)
                  agtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGAc                               >  1:4241408‑M2/31‑50 (MQ=255)
                   gtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCggatg                              <  1:1400251‑M2/21‑3 (MQ=255)
                      tctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTa                            <  1:1014819‑M2/23‑1 (MQ=255)
                      tctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTa                            <  1:2462571‑M2/23‑1 (MQ=255)
                        tgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGAc                         >  1:1900850‑M2/25‑50 (MQ=255)
                          accccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAg                           <  1:3978302‑M2/24‑1 (MQ=255)
                          accccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACt                        <  1:4206310‑M2/27‑1 (MQ=255)
                          accccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACt                        <  1:4112041‑M2/27‑1 (MQ=255)
                          accccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTg                       >  1:1081325‑M2/23‑50 (MQ=255)
                           ccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGt                      <  1:375615‑M2/29‑1 (MQ=255)
                            cccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTa                     <  1:2253012‑M2/30‑1 (MQ=255)
                            cccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTa                     >  1:3635247‑M2/21‑50 (MQ=255)
                              cgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTAt                    >  1:3510246‑M2/19‑49 (MQ=255)
                               gcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTa                  >  1:2112416‑M2/18‑50 (MQ=255)
                                cataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTa                  >  1:209282‑M2/17‑49 (MQ=255)
                                 ataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTAc                 >  1:2581757‑M2/16‑49 (MQ=255)
                                 ataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACa                <  1:2641946‑M2/35‑1 (MQ=255)
                                  taaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACa                >  1:29231‑M2/15‑49 (MQ=255)
                                   aaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACaa               >  1:1498942‑M2/14‑49 (MQ=255)
                                    aggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACaaaa             <  1:477170‑M2/38‑1 (MQ=255)
                                    aggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACaaa              <  1:896944‑M2/37‑1 (MQ=255)
                                    aggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACaa               >  1:2214033‑M2/13‑48 (MQ=255)
                                     ggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAg            <  1:1853185‑M2/39‑1 (MQ=255)
                                      gaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGc           >  1:571268‑M2/11‑50 (MQ=255)
                                        atagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCgg         >  1:2009036‑M2/9‑50 (MQ=255)
                                          agaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGc        >  1:1483786‑M2/7‑49 (MQ=255)
                                             acaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaa     >  1:1953879‑M2/4‑49 (MQ=255)
                                              caCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaaaa   <  1:475053‑M2/48‑1 (MQ=255)
                                              caCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaaaa   <  1:542660‑M2/48‑1 (MQ=255)
                                              caCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaa     >  1:4347501‑M2/3‑48 (MQ=255)
                                               aCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAg  >  1:1200463‑M2/2‑50 (MQ=255)
                                               aCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAg  >  1:685661‑M2/2‑50 (MQ=255)
                                                cAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAg  >  1:1951930/1‑49 (MQ=255)
                                                |                                                
AACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAG  >  NC_000913/2517902‑2517998

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