Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 389,898 Δ9,908 bp yaiT[yaiZ] yaiT, insF‑2, insE‑2, yaiT, iprA, ampH, sbmA, yaiW, yaiY, [yaiZ]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 389898 399805 9908 35 [0] [0] 33 yaiT–[yaiZ] yaiT,insF‑2,insE‑2,yaiT,iprA,ampH,sbmA,yaiW,yaiY,[yaiZ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3898970 (0.000)30 (0.990) 24/96 0.1 100% intergenic (‑170/‑354) hemB/yaiT porphobilinogen synthase/putative autotransporter YaiT
?NC_000913 399806 = 0 (0.000)intergenic (+1/+23) yaiZ/ddlA DUF2754 domain‑containing protein YaiZ/D‑alanine‑‑D‑alanine ligase A

GTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  NC_000913/399760‑399854
                                              |                                                
tagaaacagcaataattttatattcactgaaaatatttttaatcttGaa                                                >  1:2211955‑M2/47‑49 (MQ=255)
  gaaacagcaataattttatattcactgaaaatatttttaatcttGAata                                              >  1:686942‑M2/45‑49 (MQ=255)
  gaaacagcaataattttatattcactgaaaatatttttaatcttGAat                                               <  1:320449‑M2/4‑1 (MQ=255)
     acagcaataattttatattcactgaaaatatttttaatcttGAAtatta                                           <  1:1872481‑M2/8‑1 (MQ=255)
     acagcaataattttatattcactgaaaatatttttaatcttGAAtatta                                           <  1:1048375‑M2/8‑1 (MQ=255)
     acagcaataattttatattcactgaaaatatttttaatcttGAAtatta                                           <  1:2188292‑M2/8‑1 (MQ=255)
       agcaataattttatattcactgaaaatatttttaatcttGAATATTATc                                         >  1:1370158‑M2/40‑49 (MQ=255)
        gcaataattttatattcactgaaaatatttttaatcttGAATATTATCg                                        >  1:994724‑M2/39‑49 (MQ=255)
        gcaataattttatattcactgaaaatatttttaatcttGAATATTATCg                                        <  1:1703568‑M2/11‑1 (MQ=255)
        gcaataattttatattcactgaaaatatttttaatcttGAATATTATCg                                        <  1:2038070‑M2/11‑1 (MQ=255)
               ttttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTa                                 <  1:1233893‑M2/18‑1 (MQ=255)
                tttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTa                                 >  1:2232624‑M2/31‑48 (MQ=255)
                tttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTa                                 >  1:886144‑M2/31‑48 (MQ=255)
                 ttatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTaa                                <  1:950334‑M2/19‑1 (MQ=255)
                  tatattcactgaaaatatttttaatcttGAATATTATCGGGCGTTaaaa                              >  1:420382‑M2/29‑49 (MQ=255)
                    tattcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAata                            <  1:485005‑M2/23‑1 (MQ=255)
                    tattcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAata                            >  1:434343‑M2/27‑49 (MQ=255)
                    tattcactgaaaatatttttaatcttGAATATTATCGGGCGTTAAAat                             <  1:695642‑M2/22‑1 (MQ=255)
                           tgaaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTACAtt                     >  1:1891572‑M2/20‑49 (MQ=255)
                           tgaaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTACAtt                     >  1:1693258‑M2/20‑49 (MQ=255)
                             aaaatatttttaatcttGAATATTATCGGGCGTTAAAATATTACATtg                    <  1:2306003‑M2/31‑1 (MQ=255)
                               aatatttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTgg                 <  1:1054895‑M2/34‑1 (MQ=255)
                               aatatttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTgg                 <  1:389487‑M2/34‑1 (MQ=255)
                                 tatttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGtt               <  1:2269460‑M2/36‑1 (MQ=255)
                                  atttttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGttt              >  1:262503‑M2/13‑49 (MQ=255)
                                      ttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCaa          >  1:2188742‑M2/9‑49 (MQ=255)
                                      ttaatcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCa           <  1:2224553‑M2/40‑1 (MQ=255)
                                         atcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGc       <  1:446488‑M2/44‑1 (MQ=255)
                                         atcttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGc       >  1:642694‑M2/6‑49 (MQ=255)
                                           cttGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGt     >  1:1958028‑M2/4‑49 (MQ=255)
                                              gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTa   <  1:1211425/48‑1 (MQ=255)
                                              gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAt  >  1:1749042/1‑49 (MQ=255)
                                              gAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAt  >  1:812554/1‑49 (MQ=255)
                                              |                                                
GTTGATTATTCGTCGTCGCGATGAAGAAACTGAAAACGCCCAATAAGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  NC_000913/399760‑399854

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: 
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.