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 35 [0] [0] 35 [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)34 (0.910) 29/98 0.1 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

TGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGT  >  NC_000913/2509583‑2509679
                                               |                                                 
tGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACa                                                   >  1:2624865/1‑48 (MQ=255)
tGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAc                                                  <  1:2783601‑M1/49‑2 (MQ=255)
  cAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcag                                                <  1:571426‑M1/49‑4 (MQ=255)
   aaGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggc                                              <  1:319372‑M1/50‑6 (MQ=255)
       tGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggca                                          >  1:239831‑M1/1‑41 (MQ=255)
         ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcact                                        <  1:484268‑M1/50‑12 (MQ=255)
         ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcac                                         >  1:1185651‑M1/1‑39 (MQ=255)
         ggCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcac                                         <  1:2757845‑M1/49‑11 (MQ=255)
            aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                     <  1:1029688‑M1/50‑15 (MQ=255)
            aaaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaac                                     >  1:212598‑M1/1‑36 (MQ=255)
             aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    >  1:1111146‑M1/1‑35 (MQ=255)
             aaCTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacc                                    <  1:1429040‑M1/50‑16 (MQ=255)
               cTAAGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccg                                  <  1:250048‑M1/50‑18 (MQ=255)
                 aaGTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccgga                                >  1:2785682‑M1/1‑31 (MQ=255)
                   gTATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  1:121201‑M1/1‑29 (MQ=255)
                    tATCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgt                             <  1:1867415‑M1/50‑23 (MQ=255)
                      tCTGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacg                              >  1:720430‑M1/1‑26 (MQ=255)
                        tGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagac                         >  1:622747‑M1/1‑24 (MQ=255)
                        tGACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagac                         >  1:2711910‑M1/1‑24 (MQ=255)
                         gACCCCGCATAAGGAATAGAACAcaggcggcactaacccggacgtagact                        >  1:430376‑M1/1‑23 (MQ=255)
                           ccccGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgt                      <  1:2657249‑M1/50‑30 (MQ=255)
                              cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  1:375914‑M1/50‑33 (MQ=255)
                              cGCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatt                   <  1:2530999‑M1/50‑33 (MQ=255)
                               gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  <  1:1102194‑M1/50‑34 (MQ=255)
                               gCATAAGGAATAGAACAcaggcggcactaacccggacgtagactgtatta                  >  1:1749393‑M1/1‑17 (MQ=255)
                                 aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca                <  1:426528‑M1/50‑36 (MQ=255)
                                 aTAAGGAATAGAACAcaggcggcactaacccggacgtagactgtattaca                <  1:705294‑M1/50‑36 (MQ=255)
                                   aaGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               >  1:245001‑M1/1‑13 (MQ=255)
                                   aaGGAATAGAACAcaggcggcactaacccggacgtagactgtattacaa               >  1:801942‑M1/1‑13 (MQ=255)
                                            aaCAcaggcggcactaacccggacgtagactgtattacaaaagcggcaa      >  1:311715‑M1/1‑4 (MQ=255)
                                            aaCAcaggcggcactaacccggacgtagactgtattacaaaagcggcaa      >  1:146847‑M1/1‑4 (MQ=255)
                                             acacaggcggcactaacccggacgtagactgtattacaaaagcggcaaaa    <  1:2345909‑M1/50‑48 (MQ=255)
                                              cacaggcggcactaacccggacgtagactgtattacaaaagcggcaaa     >  1:2135483‑M1/1‑2 (MQ=255)
                                               acaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaag  >  1:217701‑M1/1‑1 (MQ=255)
                                               acaggcggcactaacccggacgtagactgtattacaaaagcggcaaaaag  <  1:1143055‑M1/50‑50 (MQ=255)
                                               |                                                 
TGCAAGCTGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGT  >  NC_000913/2509583‑2509679

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: