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 97 [0] [0] 99 [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)97 (1.030) 59/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

GCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGT  >  NC_000913/2509584‑2509679
                                              |                                                 
gCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAca                                                 <  1:7070578‑M1/49‑3 (MQ=255)
gCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAca                                                 <  1:1158509‑M1/49‑3 (MQ=255)
gCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAca                                                 <  1:4219046‑M1/49‑3 (MQ=255)
gCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAca                                                 <  1:2248738‑M1/49‑3 (MQ=255)
 cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               <  1:3033970‑M1/50‑5 (MQ=255)
 cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               <  1:4978065‑M1/50‑5 (MQ=255)
 cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               >  1:4929270‑M1/1‑46 (MQ=255)
 cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               >  1:5626969‑M1/1‑46 (MQ=255)
  aaGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggc                                              >  1:6244783‑M1/1‑45 (MQ=255)
  aaGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggc                                              <  1:3271446‑M1/50‑6 (MQ=255)
  aaGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggc                                              <  1:5321652‑M1/50‑6 (MQ=255)
  aaGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggc                                              <  1:5390641‑M1/50‑6 (MQ=255)
   aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                             <  1:4513702‑M1/50‑7 (MQ=255)
   aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                             <  1:4358109‑M1/50‑7 (MQ=255)
     cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggc                                           <  1:3088746‑M1/50‑9 (MQ=255)
     cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggc                                           <  1:1026580‑M1/50‑9 (MQ=255)
     cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                            >  1:1705928‑M1/1‑42 (MQ=255)
      tGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggca                                          >  1:401928‑M1/1‑41 (MQ=255)
        ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                        <  1:334728‑M1/50‑12 (MQ=255)
        ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                        <  1:6745475‑M1/50‑12 (MQ=255)
        ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcac                                         >  1:6508498‑M1/1‑39 (MQ=255)
        ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggca                                          <  1:4656247‑M1/48‑10 (MQ=255)
         gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcacta                                       <  1:6582050‑M1/50‑13 (MQ=255)
         gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                        >  1:3519691‑M1/1‑38 (MQ=255)
           aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcgctaac                                     <  1:2322725‑M1/50‑15 (MQ=255)
           aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                     >  1:528289‑M1/1‑36 (MQ=255)
           aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                     >  1:2812003‑M1/1‑36 (MQ=255)
            aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    >  1:5372180‑M1/1‑35 (MQ=255)
            aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    >  1:4000763‑M1/1‑35 (MQ=255)
            aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    >  1:398688‑M1/1‑35 (MQ=255)
             aCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccc                                   <  1:4200427‑M1/50‑17 (MQ=255)
              cTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  <  1:4539607‑M1/50‑18 (MQ=255)
               tAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  >  1:67429‑M1/1‑32 (MQ=255)
               tAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  >  1:6375277‑M1/1‑32 (MQ=255)
                 aGTATCTGACCCCGCATAAGGGATAGAACAcaggcggcactaacccggac                               <  1:3196678‑M1/50‑21 (MQ=255)
                 aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               <  1:3904482‑M1/50‑21 (MQ=255)
                 aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               <  1:6332519‑M1/50‑21 (MQ=255)
                 aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               <  1:4476307‑M1/50‑21 (MQ=255)
                 aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               >  1:2247222‑M1/1‑30 (MQ=255)
                  gTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              <  1:6697332‑M1/50‑22 (MQ=255)
                  gTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              <  1:5669036‑M1/50‑22 (MQ=255)
                   tATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  1:5345522‑M1/1‑28 (MQ=255)
                   tATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  1:2854232‑M1/1‑28 (MQ=255)
                     tCTGGCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgta                            <  1:5394452‑M1/49‑24 (MQ=255)
                     tCTGGCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgta                            <  1:4528986‑M1/49‑24 (MQ=255)
                     tCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgta                            <  1:4900200‑M1/49‑24 (MQ=255)
                     tCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgta                            <  1:4438024‑M1/49‑24 (MQ=255)
                     tCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  1:4801453‑M1/1‑26 (MQ=255)
                       tGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagac                         >  1:5314751‑M1/1‑24 (MQ=255)
                        gACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagact                        >  1:6672000‑M1/1‑23 (MQ=255)
                         aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtaga                          <  1:1088899‑M1/47‑26 (MQ=255)
                         aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtaga                          <  1:2372402‑M1/47‑26 (MQ=255)
                         aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtag                           <  1:3107276‑M1/46‑25 (MQ=255)
                         aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtag                           <  1:569‑M1/46‑25 (MQ=255)
                           cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                     >  1:6341984‑M1/1‑20 (MQ=255)
                           cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                     >  1:3263930‑M1/1‑20 (MQ=255)
                           cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                     <  1:6075590‑M1/50‑31 (MQ=255)
                            ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    >  1:181167‑M1/1‑19 (MQ=255)
                            ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    <  1:5151781‑M1/50‑32 (MQ=255)
                             cGCATAAGGAATAGAACAcaggcggcactagcccggacgtagactgtatt                   <  1:1028912‑M1/50‑33 (MQ=255)
                             cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  1:3222977‑M1/50‑33 (MQ=255)
                             cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    >  1:5169136‑M1/1‑18 (MQ=255)
                              gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  <  1:6895668‑M1/50‑34 (MQ=255)
                              gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  1:1750775‑M1/1‑17 (MQ=255)
                              gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  <  1:5687818‑M1/50‑34 (MQ=255)
                               cATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattac                 <  1:4963455‑M1/50‑35 (MQ=255)
                               cATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  1:964960‑M1/1‑16 (MQ=255)
                               cATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  1:110985‑M1/1‑16 (MQ=255)
                                aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattac                 >  1:2355108‑M1/1‑15 (MQ=255)
                                aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattac                 >  1:7077151‑M1/1‑15 (MQ=255)
                                 tAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca                >  1:4622950‑M1/1‑14 (MQ=255)
                                  aaGGAATAGAACAcaggcggcactaacccggacgcagactgtattacaa               >  1:1740098‑M1/1‑13 (MQ=255)
                                   aGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaa             <  1:4571708‑M1/50‑39 (MQ=255)
                                   aGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               >  1:5105972‑M1/1‑12 (MQ=255)
                                   aGGAATAGAACAcaggcggcactaacccggacgaagactgtattacaa               >  1:5997362‑M1/1‑12 (MQ=255)
                                     gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           >  1:784340‑M1/1‑10 (MQ=255)
                                     gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           >  1:6469378‑M1/1‑10 (MQ=255)
                                     gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           <  1:690559‑M1/50‑41 (MQ=255)
                                     gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           >  1:7111236‑M1/1‑10 (MQ=255)
                                     gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           >  1:5591438‑M1/1‑10 (MQ=255)
                                      aaTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcg          >  1:6921435‑M1/1‑9 (MQ=255)
                                      aaTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcg          <  1:2101382‑M1/50‑42 (MQ=255)
                                       aTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcgg         >  1:2489387‑M1/1‑8 (MQ=255)
                                       aTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcgg         >  1:4489537‑M1/1‑8 (MQ=255)
                                        tAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcggc        <  1:1539763‑M1/50‑44 (MQ=255)
                                          gAACAcaggcggcactaacccggacgtagactgtattacaaaagcggca       >  1:4104891‑M1/1‑5 (MQ=255)
                                           aaCAcaggcggcactaacccggacgtagactgtattacaaaagcggcaaa     <  1:3645623‑M1/50‑47 (MQ=255)
                                            acacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaa    <  1:5058506‑M1/50‑48 (MQ=255)
                                            acacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaa    <  1:5726499‑M1/50‑48 (MQ=255)
                                            acacaggcggcactaacccggacgtagactgtattacaaaagcggcaaa     >  1:3986156‑M1/1‑3 (MQ=255)
                                             cacaggcggcactaacccggacgtagactttattacaaaagcggcaaaaa   <  1:4619802‑M1/50‑49 (MQ=255)
                                             cacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaa   <  1:2369690‑M1/50‑49 (MQ=255)
                                             cacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaa   <  1:2944271‑M1/50‑49 (MQ=255)
                                             cacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaa   <  1:5559784‑M1/50‑49 (MQ=255)
                                             cacaggcggcactaacccggacgtagactgtattacaaaagcggcaaa     >  1:4263609‑M1/1‑2 (MQ=255)
                                              acaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaag  >  1:4178408‑M1/1‑1 (MQ=255)
                                              acaggcggcactaacccggacgtagactgtattacaaaagcggaaaaagc  <  1:167623‑M1/50‑50 (MQ=255)
                                              |                                                 
GCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGT  >  NC_000913/2509584‑2509679

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