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 48 [0] [0] 51 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)48 (1.580) 42/98 0.0 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

AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGC  >  NC_000913/389850‑389943
                                               |                                              
aaTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTg                                               <  1:499936‑M1/49‑2 (MQ=255)
 aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaa                                             >  1:2106533‑M1/1‑47 (MQ=255)
 aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTga                                              <  1:2063767‑M1/49‑3 (MQ=255)
 aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTga                                              <  1:2255480‑M1/49‑3 (MQ=255)
 aTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTga                                              <  1:2607889‑M1/49‑3 (MQ=255)
  tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaa                                             >  1:305744‑M1/1‑46 (MQ=255)
   aGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaata                                           >  1:1184467‑M1/1‑45 (MQ=255)
    gAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaata                                           <  1:2644294‑M1/49‑6 (MQ=255)
     aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatatt                                         <  1:850569‑M1/50‑8 (MQ=255)
     aaaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatatt                                         <  1:746279‑M1/50‑8 (MQ=255)
      aaCAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatatta                                        >  1:427265‑M1/1‑42 (MQ=255)
         aGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcg                                     <  1:2439819‑M1/50‑12 (MQ=255)
         aGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcg                                     <  1:366025‑M1/50‑12 (MQ=255)
         aGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcg                                     >  1:2211946‑M1/1‑39 (MQ=255)
          gCAATAATTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcgg                                    >  1:2036241‑M1/1‑38 (MQ=255)
            aataatTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcgggc                                  <  1:538114‑M1/50‑15 (MQ=255)
             ataatTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcg                                 >  1:270801‑M1/1‑35 (MQ=255)
             ataatTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcg                                 >  1:543475‑M1/1‑35 (MQ=255)
             ataatTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcg                                 <  1:942825‑M1/50‑16 (MQ=255)
                atTTTATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcgtta                              <  1:379925‑M1/50‑19 (MQ=255)
                 ttttATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcgtta                              >  1:2539647‑M1/1‑31 (MQ=255)
                  tttATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaa                             >  1:888088‑M1/1‑30 (MQ=255)
                   ttATATTCACTGAAAATATTTTTAATCTTgaatattgtcgggcgttaaaa                           >  1:1044471‑M1/1‑29 (MQ=255)
                   ttATATTCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaaa                            <  1:684352‑M1/49‑21 (MQ=255)
                     atatTCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaaaata                         >  1:1279329‑M1/1‑27 (MQ=255)
                     atatTCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaaaat                          <  1:582494‑M1/49‑23 (MQ=255)
                      tatTCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaaaatat                        >  1:520838‑M1/1‑26 (MQ=255)
                        ttCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaaaatatt                       <  1:1538132‑M1/49‑26 (MQ=255)
                         tCACTGAAAATATTTTTAATCTTgaatattatcgggcgttaaaatatta                      >  1:2558072‑M1/1‑23 (MQ=255)
                           aCTGAAAATATTTTTAATCTTgaatattatcgggcgttaaaatattacat                   >  1:1597887‑M1/1‑21 (MQ=255)
                              gAAAATATTTTTAATCTTgaatattatcgggcgttaaaatattacattgt                >  1:1353864‑M1/1‑18 (MQ=255)
                               aaaaTATTTTTAATCTTgaatattatcgggcgttaaaatattacattgtg               >  1:1179766‑M1/1‑17 (MQ=255)
                                  atatTTTTAATCTTgaatattatcgggcgttaaaatattacattgtggtt            >  1:854411‑M1/1‑14 (MQ=255)
                                  atatTTTTAATCTTgaatattatcgggcgttaaaatattacattgtggtt            <  1:1064033‑M1/50‑37 (MQ=255)
                                   tatTTTTAATCTTgaatattatcgggcgttaaaatattacattgtggttt           <  1:2303259‑M1/50‑38 (MQ=255)
                                    atTTTTAATCTTgaatattatcgggcgttaaaatattacattgtggtttt          >  1:1101872‑M1/1‑12 (MQ=255)
                                     tttttAATCTTgaatattatcgggcgttaaaatattacattgtggttttc         <  1:1647042‑M1/50‑40 (MQ=255)
                                     tttttAATCTTgaatattatcgggcgttaaaatattacattgtggtttt          >  1:964065‑M1/1‑11 (MQ=255)
                                     tttttAATCTTgaatattatcgggcgttaaaatattacattgtggtttt          <  1:1650184‑M1/49‑39 (MQ=255)
                                       tttAATCTTgaatattatcgggcgttaaaatattacattgtggttttcaa       >  1:734475‑M1/1‑9 (MQ=255)
                                       tttAATCTTgaatattatcgggcgttaaaatattacattgtggttttca        <  1:654456‑M1/49‑41 (MQ=255)
                                        ttAATCTTgaatattatcgggcgttaaaatattacattgtggttttcaa       <  1:2456554‑M1/49‑42 (MQ=255)
                                         tAATCTTgaatattatcgggcgttaaaatattacattgtggttttcaat      <  1:1493844‑M1/49‑43 (MQ=255)
                                         tAATCTTgaatattatcgggcgttaaaatattacattgtggttttcaat      <  1:1061396‑M1/49‑43 (MQ=255)
                                          aaTCTTgaatattatcgggcgttaaaatattacattgtggttttcaatgc    >  1:2622717‑M1/1‑6 (MQ=255)
                                           aTCTTgaatattatcgggcgttaaaatattacattgtggttttcaatgcg   <  1:2137610‑M1/50‑46 (MQ=255)
                                            tCTTgaatattatcgggcgttaaaatattacattgtggttttcaatgcgt  >  1:2050696‑M1/1‑4 (MQ=255)
                                             cTTgaatattatcgggcgttaaaatattacattgtggttttcaatgcgt  >  1:2268651‑M1/1‑3 (MQ=255)
                                               |                                              
AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGC  >  NC_000913/389850‑389943

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