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 128 [0] [0] 128 [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)124 (0.830) 65/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

Only 100 of 128 total aligned reads displayed.

TTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGT  >  NC_000913/2509581‑2509679
                                                 |                                                 
tttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                   >  1:3685046/1‑50 (MQ=255)
 ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                   <  1:73408/49‑1 (MQ=255)
 ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                   <  1:4312115/49‑1 (MQ=255)
 ttGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                   >  2:1475846/1‑49 (MQ=255)
  tGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAc                                                  <  1:5391399‑M1/49‑2 (MQ=255)
   gCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                >  2:2379318‑M1/1‑47 (MQ=255)
   gCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAca                                                 <  1:5303164‑M1/49‑3 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               <  1:1127229‑M1/50‑5 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               >  1:3252880‑M1/1‑46 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcagg                                               <  1:771085‑M1/50‑5 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  2:3423636‑M1/49‑4 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  1:2035846‑M1/49‑4 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  2:4512651‑M1/49‑4 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  2:2135779‑M1/49‑4 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  2:3770942‑M1/49‑4 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  1:1326741‑M1/49‑4 (MQ=255)
    cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  1:1105003‑M1/49‑4 (MQ=255)
     aaGCTGGGCAAACTAAGTATCTGACCCCGTATAAGGAATAGAACAcaggc                                              <  1:908955‑M1/50‑6 (MQ=255)
      aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                             >  2:3418356‑M1/1‑44 (MQ=255)
      aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                             <  1:4442674‑M1/50‑7 (MQ=255)
      aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                             <  2:3334762‑M1/50‑7 (MQ=255)
      aGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcg                                             <  2:4294192‑M1/50‑7 (MQ=255)
       gCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                            >  2:3833828‑M1/1‑43 (MQ=255)
       gCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                            >  2:2212680‑M1/1‑43 (MQ=255)
       gCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                            <  1:5138094‑M1/50‑8 (MQ=255)
        cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggc                                           <  2:3116755‑M1/50‑9 (MQ=255)
        cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                            >  1:4959885‑M1/1‑42 (MQ=255)
        cTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcgg                                            >  2:4149617‑M1/1‑42 (MQ=255)
         tGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggca                                          >  2:3294237‑M1/1‑41 (MQ=255)
         tGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggca                                          >  1:3584593‑M1/1‑41 (MQ=255)
         tGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggca                                          >  2:4937160‑M1/1‑41 (MQ=255)
           ggCAAACTAAGTATCTGACCCCGCATAATGAATAGAACAcaggcggcac                                         >  2:2276113‑M1/1‑39 (MQ=255)
           ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcac                                         >  2:5367571‑M1/1‑39 (MQ=255)
           ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcac                                         <  2:1232897‑M1/49‑11 (MQ=255)
            gCAAACTAAGTATATGCCCCCGCATAAGGAATAGAACAaaggcggcactc                                       <  2:5320412‑M1/50‑13 (MQ=255)
             cAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaa                                      <  1:5697350‑M1/50‑14 (MQ=255)
             cAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaa                                      >  2:4154310‑M1/1‑37 (MQ=255)
             cAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaa                                      <  1:3562004‑M1/50‑14 (MQ=255)
              aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                     >  2:3040529‑M1/1‑36 (MQ=255)
              aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                     >  2:1547665‑M1/1‑36 (MQ=255)
               aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    <  2:1717864‑M1/50‑16 (MQ=255)
               aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    <  1:2427452‑M1/50‑16 (MQ=255)
                aCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaaccc                                   <  2:3258641‑M1/50‑17 (MQ=255)
                 cTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  <  2:1060480‑M1/50‑18 (MQ=255)
                 cTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  <  2:5076629‑M1/50‑18 (MQ=255)
                  tAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  >  1:4178659‑M1/1‑32 (MQ=255)
                   aaGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccgga                                >  1:3939928‑M1/1‑31 (MQ=255)
                   aaGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccgga                                >  1:3921587‑M1/1‑31 (MQ=255)
                    aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               <  1:2814895‑M1/50‑21 (MQ=255)
                    aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               <  1:1943472‑M1/50‑21 (MQ=255)
                    aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               >  1:2315266‑M1/1‑30 (MQ=255)
                    aGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggac                               >  2:342040‑M1/1‑30 (MQ=255)
                      tATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  2:3762074‑M1/1‑28 (MQ=255)
                        tCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  2:1003689‑M1/1‑26 (MQ=255)
                           gACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagact                        >  2:5275675‑M1/1‑23 (MQ=255)
                           gACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagact                        >  2:546686‑M1/1‑23 (MQ=255)
                            aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtaga                          <  1:3248382‑M1/47‑26 (MQ=255)
                            aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtag                           <  2:4591284‑M1/46‑25 (MQ=255)
                            aCCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtag                           <  2:4353215‑M1/46‑25 (MQ=255)
                              cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgta                     <  1:4194819‑M1/50‑31 (MQ=255)
                              cccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagacggta                     <  2:4721789‑M1/50‑31 (MQ=255)
                               ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    >  1:479414‑M1/1‑19 (MQ=255)
                               ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    <  1:2754838‑M1/50‑32 (MQ=255)
                               ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    >  2:5370464‑M1/1‑19 (MQ=255)
                               ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    <  2:2242219‑M1/50‑32 (MQ=255)
                               ccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    >  1:1607585‑M1/1‑19 (MQ=255)
                                cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  2:1125922‑M1/50‑33 (MQ=255)
                                cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  1:730193‑M1/50‑33 (MQ=255)
                                cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  2:1217239‑M1/50‑33 (MQ=255)
                                cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  2:1671624‑M1/50‑33 (MQ=255)
                                cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtat                    >  2:2925217‑M1/1‑18 (MQ=255)
                                 gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  <  1:4685600‑M1/50‑34 (MQ=255)
                                 gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  1:2555110‑M1/1‑17 (MQ=255)
                                 gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  1:2359032‑M1/1‑17 (MQ=255)
                                 gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  <  1:2251293‑M1/50‑34 (MQ=255)
                                  cATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  2:2599823‑M1/1‑16 (MQ=255)
                                   aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca                <  1:2639597‑M1/50‑36 (MQ=255)
                                   aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattac                 >  1:3526930‑M1/1‑15 (MQ=255)
                                    tAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               <  2:1755467‑M1/50‑37 (MQ=255)
                                    tAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca                >  1:3056733‑M1/1‑14 (MQ=255)
                                    tAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca                >  1:3968461‑M1/1‑14 (MQ=255)
                                    tAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaaaa               <  1:5519893‑M1/50‑37 (MQ=255)
                                     aaGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               >  2:458530‑M1/1‑13 (MQ=255)
                                      aGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               >  2:814464‑M1/1‑12 (MQ=255)
                                      aGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               >  2:169426‑M1/1‑12 (MQ=255)
                                       ggAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaag            >  1:3817128‑M1/1‑11 (MQ=255)
                                       ggAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaag            >  1:231025‑M1/1‑11 (MQ=255)
                                        gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           <  1:46041‑M1/50‑41 (MQ=255)
                                        gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           >  2:3840992‑M1/1‑10 (MQ=255)
                                        gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           >  1:5055271‑M1/1‑10 (MQ=255)
                                        gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           <  2:316828‑M1/50‑41 (MQ=255)
                                        gAATAGAACAcaggcggcactaacccggacgtagactgtattacaaaagc           <  1:1228698‑M1/50‑41 (MQ=255)
                                          aTAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcgg         >  1:1815644‑M1/1‑8 (MQ=255)
                                           tAGAACAcaggcggcactaacccggacgtagactgtattacaaaagcggc        <  1:3659053‑M1/50‑44 (MQ=255)
                                             gAACAcaggcggcactaacccggacgtagactgtattacaaaagcggca       <  2:1757482‑M1/49‑45 (MQ=255)
                                              aaCAcaggcggcactaacccggacgtagactgtattacaaaagcggcaaa     <  2:2490537‑M1/50‑47 (MQ=255)
                                               acacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaa    <  1:4976167‑M1/50‑48 (MQ=255)
                                                 acaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaag  >  2:1168839‑M1/1‑1 (MQ=255)
                                                 acaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaag  >  2:5493217‑M1/1‑1 (MQ=255)
                                                 acaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaag  <  2:2752012‑M1/50‑50 (MQ=255)
                                                 |                                                 
TTTGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGT  >  NC_000913/2509581‑2509679

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