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 80 [0] [0] 75 [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)74 (0.830) 51/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

TTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCG  >  NC_000913/2509582‑2509678
                                                |                                                
ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                  <  1:2489799/49‑1 (MQ=255)
ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                  <  1:1553945/49‑1 (MQ=255)
ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                  <  1:3515490/49‑1 (MQ=255)
ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                  <  1:3512322/49‑1 (MQ=255)
ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                  >  1:1662701/1‑49 (MQ=255)
 tGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                  >  1:6693136/1‑48 (MQ=255)
 tGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAc                                                 <  1:4076988‑M1/49‑2 (MQ=255)
     aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                            <  1:4235808‑M1/50‑7 (MQ=255)
     aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                            >  1:3272818‑M1/1‑44 (MQ=255)
       cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                           >  1:1808471‑M1/1‑42 (MQ=255)
           gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcacta                                      <  1:2204554‑M1/50‑13 (MQ=255)
           gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcacta                                      <  1:2143114‑M1/50‑13 (MQ=255)
           gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                       >  1:157225‑M1/1‑38 (MQ=255)
           gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                       >  1:3616802‑M1/1‑38 (MQ=255)
           gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                       >  1:6291621‑M1/1‑38 (MQ=255)
           gCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                       >  1:6405932‑M1/1‑38 (MQ=255)
            cAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaa                                     >  1:6667472‑M1/1‑37 (MQ=255)
            cAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaa                                     <  1:4086283‑M1/50‑14 (MQ=255)
            cAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaa                                     >  1:6469691‑M1/1‑37 (MQ=255)
             aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                    <  1:320939‑M1/50‑15 (MQ=255)
             aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                    >  1:2752180‑M1/1‑36 (MQ=255)
              aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                   >  1:4018175‑M1/1‑35 (MQ=255)
               aCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccc                                  <  1:6734236‑M1/50‑17 (MQ=255)
               aCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccc                                  <  1:4547497‑M1/50‑17 (MQ=255)
               aCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccc                                  >  1:2685215‑M1/1‑34 (MQ=255)
               aCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccc                                  >  1:830450‑M1/1‑34 (MQ=255)
                cTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccct                                 >  1:3576074‑M1/1‑33 (MQ=255)
                cTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                 <  1:1357173‑M1/50‑18 (MQ=255)
                 tAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                 >  1:575479‑M1/1‑32 (MQ=255)
                 tAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                 >  1:2732988‑M1/1‑32 (MQ=255)
                   aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                              <  1:4095547‑M1/50‑21 (MQ=255)
                   aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                              <  1:1965174‑M1/50‑21 (MQ=255)
                   aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                              <  1:355973‑M1/50‑21 (MQ=255)
                   aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccgga                               <  1:6598098‑M1/49‑20 (MQ=255)
                   aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccgg                                <  1:3823106‑M1/48‑19 (MQ=255)
                    gTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                             >  1:6560395‑M1/1‑29 (MQ=255)
                     tATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                             >  1:5651959‑M1/1‑28 (MQ=255)
                     tATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                             >  1:3622358‑M1/1‑28 (MQ=255)
                      aTCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                             >  1:4183629‑M1/1‑27 (MQ=255)
                      aTCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                             >  1:428075‑M1/1‑27 (MQ=255)
                       tCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                             >  1:1643231‑M1/1‑26 (MQ=255)
                       tCTGACCCCGCATAAGGAATAGAACAcaggcggcacgaacccggacgta                           <  1:5308958‑M1/49‑24 (MQ=255)
                         tGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagac                        >  1:1622897‑M1/1‑24 (MQ=255)
                           aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactg                      >  1:4598731‑M1/1‑22 (MQ=255)
                           aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactg                      >  1:1334420‑M1/1‑22 (MQ=255)
                           aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactg                      <  1:1974920‑M1/50‑29 (MQ=255)
                           aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagac                        <  1:5305913‑M1/48‑27 (MQ=255)
                           aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtag                          <  1:2192495‑M1/46‑25 (MQ=255)
                            ccccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgt                     <  1:5979675‑M1/50‑30 (MQ=255)
                            ccccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactg                      >  1:5720305‑M1/1‑21 (MQ=255)
                             cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                    >  1:5957093‑M1/1‑20 (MQ=255)
                             cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                    >  1:2845880‑M1/1‑20 (MQ=255)
                             cccGAATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                    >  1:3006419‑M1/1‑20 (MQ=255)
                             cccACATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                    <  1:6083478‑M1/50‑31 (MQ=255)
                              ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   >  1:6411348‑M1/1‑19 (MQ=255)
                              ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   <  1:324115‑M1/50‑32 (MQ=255)
                              ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   <  1:2037325‑M1/50‑32 (MQ=255)
                              ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   <  1:689067‑M1/50‑32 (MQ=255)
                              ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   >  1:1258334‑M1/1‑19 (MQ=255)
                               cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                  <  1:914906‑M1/50‑33 (MQ=255)
                               cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   >  1:5136435‑M1/1‑18 (MQ=255)
                               cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                   >  1:4767975‑M1/1‑18 (MQ=255)
                                gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                 <  1:1568528‑M1/50‑34 (MQ=255)
                                 cATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattac                <  1:1436825‑M1/50‑35 (MQ=255)
                                 cATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattac                <  1:6209306‑M1/50‑35 (MQ=255)
                                  aTAAGGAATAGAACAcaggcggcgctaacccggacgtagactgtattaca               <  1:303490‑M1/50‑36 (MQ=255)
                                  aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca               <  1:2345500‑M1/50‑36 (MQ=255)
                                   tAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca               >  1:1234047‑M1/1‑14 (MQ=255)
                                    aaGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa              >  1:1825962‑M1/1‑13 (MQ=255)
                                     aGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa              >  1:5661421‑M1/1‑12 (MQ=255)
                                     aGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa              >  1:877946‑M1/1‑12 (MQ=255)
                                      ggAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaag           >  1:5843051‑M1/1‑11 (MQ=255)
                                      ggAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaag           <  1:1084841‑M1/50‑40 (MQ=255)
                                        aaTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcg         <  1:94627‑M1/50‑42 (MQ=255)
                                         aTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcgg        >  1:232667‑M1/1‑8 (MQ=255)
                                            gAACAcaggcggcactaacccggacgtagactgtattacaaaagcggcaa     <  1:2667453‑M1/50‑46 (MQ=255)
                                             aaCAcaggcggcactaacccggacgtagactgtattacaaaagcggcaaa    <  1:487746‑M1/50‑47 (MQ=255)
                                             aaCAcaggcggcactaacccggacgtagactgtattacaaaagcggcaa     >  1:2006075‑M1/1‑4 (MQ=255)
                                              acacaggcggcactaacccggacgtagactgtattacaaaagcggcaaa    >  1:6662226‑M1/1‑3 (MQ=255)
                                               cacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaa  <  1:2968477‑M1/50‑49 (MQ=255)
                                                |                                                
TTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCG  >  NC_000913/2509582‑2509678

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: