Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,388,889 Δ9,734 bp ycjG[uspE] ycjG, mpaA, ymjC, ymjD, ycjY, pgrR, mppA, ynaI, insH‑4, ynaJ, [uspE]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1388889 1398622 9734 24 [0] [0] 24 ycjG–[uspE] ycjG,mpaA,ymjC,ymjD,ycjY,pgrR,mppA,ynaI,insH‑4,ynaJ,[uspE]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 13888880 (0.000)23 (0.860) 20/96 0.1 100% intergenic (‑77/‑42) tpx/ycjG lipid hydroperoxide peroxidase/L‑Ala‑D/L‑Glu epimerase
?NC_000913 1398623 = 0 (0.000)intergenic (‑1/+151) uspE/fnr universal stress protein with a role cellular motility/DNA‑binding transcriptional dual regulator FNR

TCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATT  >  NC_000913/1398576‑1398671
                                               |                                                
gtcacacctagtttatgctaactgtcaataacacagcaaacgctatta                                                  >  1:2325116‑M2/48‑48 (MQ=255)
 tcacacctagtttatgctaactgtcaataacacagcaaacgctattACAg                                               <  1:2104027‑M2/4‑1 (MQ=255)
   acacctagtttatgctaactgtcaataacacagcaaacgctattACAggg                                             >  1:1263894‑M2/45‑50 (MQ=255)
      cctagtttatgctaactgtcaataacacagcaaacgctattACAGGGtct                                          >  1:2164304‑M2/42‑50 (MQ=255)
          gtttatgctaactgtcaataacacagcaaacgctattACAGGGTCTcc                                        <  1:918453‑M2/11‑1 (MQ=255)
          gtttatgctaactgtcaataacacagcaaacgctattACAGGGTCTcc                                        <  1:233766‑M2/11‑1 (MQ=255)
          gtttatgctaactgtcaataacacagcaaacgctattACAGGGTCTCCt                                       >  1:1616696‑M2/38‑49 (MQ=255)
             tatgctaactgtcaataacacagcaaacgctattACAGGGTCTCCTTac                                    <  1:1751599‑M2/15‑1 (MQ=255)
                    actgtcaataacacagcaaacgctattACAGGGTCTCCTTACAACAACt                             <  1:1480629‑M2/22‑1 (MQ=255)
                         caataacacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAAc                       >  1:269885‑M2/23‑50 (MQ=255)
                          aataacacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACg                      <  1:2059291‑M2/29‑1 (MQ=255)
                           ataacacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACGc                     <  1:168591‑M2/30‑1 (MQ=255)
                            taacacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACGCa                    <  1:1522789‑M2/31‑1 (MQ=255)
                             aacacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACGCAg                   >  1:1525526‑M2/19‑50 (MQ=255)
                               cacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACGCAGtt                 >  1:734489‑M2/17‑50 (MQ=255)
                               cacagcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACGCAGtt                 <  1:76007‑M2/34‑1 (MQ=255)
                                  agcaaacgctattACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTg               >  1:2262319‑M2/14‑49 (MQ=255)
                                    caaacgctattACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTaa            <  1:2344842‑M2/39‑1 (MQ=255)
                                            attACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaa        >  1:470518‑M2/4‑46 (MQ=255)
                                             ttACAGGGTCTCCTTGCAACAACTGTCAACGCAGTTTGTAATTAAAAGa    <  1:487769‑M2/47‑1 (MQ=255)
                                             ttACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAg     <  1:289608‑M2/46‑1 (MQ=255)
                                             ttACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAg     <  1:714460‑M2/46‑1 (MQ=255)
                                             ttACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAg     <  1:2052753‑M2/46‑1 (MQ=255)
                                               aCAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGAtt  >  1:520013/1‑49 (MQ=255)
                                               |                                                
TCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATT  >  NC_000913/1398576‑1398671

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