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 42 [0] [0] 42 [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)41 (0.800) 34/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

ACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGC  >  NC_000913/2517903‑2517999
                                               |                                                 
gcaagctgggcaaactaagtatctgaccccgcataaggaatagaacaCAg                                                 >  1:1165764‑M2/48‑50 (MQ=255)
  aagctgggcaaactaagtatctgaccccgcctaaggaatagaacaCAggc                                               >  1:1324663‑M2/46‑50 (MQ=255)
   agctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcg                                              <  1:523325‑M2/6‑1 (MQ=255)
    gctgggcaaactaagtatctgaccccgcataaggaatagaacaCAggcgg                                             >  1:2945353‑M2/44‑50 (MQ=255)
        ggcaaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCa                                           <  1:3640162‑M2/9‑1 (MQ=255)
          caaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTaa                                       >  1:1136174‑M2/38‑50 (MQ=255)
          caaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTaa                                       >  1:2503675‑M2/38‑50 (MQ=255)
           aaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAc                                      >  1:1073027‑M2/37‑50 (MQ=255)
           aaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAc                                      <  1:1535325‑M2/14‑1 (MQ=255)
           aaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAc                                      <  1:1868628‑M2/14‑1 (MQ=255)
           aaactaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAc                                      >  1:466833‑M2/37‑50 (MQ=255)
             actaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAccc                                    <  1:588384‑M2/16‑1 (MQ=255)
             actaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAAccc                                    >  1:530485‑M2/35‑50 (MQ=255)
              ctaagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCg                                   >  1:732334‑M2/34‑50 (MQ=255)
               taagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCg                                   >  1:3510752‑M2/33‑49 (MQ=255)
                aagtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGa                                 >  1:3271289‑M2/32‑50 (MQ=255)
                 agtatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGAc                                <  1:1064386‑M2/20‑1 (MQ=255)
                   tatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACg                               >  1:342008‑M2/29‑49 (MQ=255)
                   tatctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGt                              <  1:302132‑M2/22‑1 (MQ=255)
                     tctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTa                             <  1:1201401‑M2/23‑1 (MQ=255)
                      ctgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGa                           >  1:1849949‑M2/26‑50 (MQ=255)
                       tgaccccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGAc                          >  1:1314843‑M2/25‑50 (MQ=255)
                         accccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGAc                          <  1:138018‑M2/26‑1 (MQ=255)
                         accccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACt                         <  1:2034707‑M2/27‑1 (MQ=255)
                            ccgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTAt                     >  1:3400253‑M2/20‑50 (MQ=255)
                             cgcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTAt                     >  1:3076021‑M2/19‑49 (MQ=255)
                              gcataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTa                   >  1:598214‑M2/18‑50 (MQ=255)
                                ataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTAc                  >  1:3004657‑M2/16‑49 (MQ=255)
                                ataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTAc                  >  1:457951‑M2/16‑49 (MQ=255)
                                ataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTAc                  >  1:502817‑M2/16‑49 (MQ=255)
                                ataaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTAc                  >  1:1172088‑M2/16‑49 (MQ=255)
                                  aaggaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACaa                >  1:1455236‑M2/14‑49 (MQ=255)
                                     gaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGc            <  1:1863262‑M2/40‑1 (MQ=255)
                                     gaatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGc            >  1:771823‑M2/11‑50 (MQ=255)
                                      aatagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCg           <  1:1591938‑M2/41‑1 (MQ=255)
                                       atagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCgg          >  1:626689‑M2/9‑50 (MQ=255)
                                       atagaacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCgg          >  1:847139‑M2/9‑50 (MQ=255)
                                           aacaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaa      <  1:627785‑M2/46‑1 (MQ=255)
                                            acaCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaa      >  1:962675‑M2/4‑49 (MQ=255)
                                             caCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCaaa      >  1:3608509‑M2/3‑48 (MQ=255)
                                              aCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAg   <  1:3876694‑M2/49‑1 (MQ=255)
                                               cAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGc  <  1:1597869/50‑1 (MQ=255)
                                               |                                                 
ACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGC  >  NC_000913/2517903‑2517999

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