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] 36 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)34 (1.240) 27/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

AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/389850‑389897
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  NC_000913/399806‑399854
                                                                                                 
AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTG                                                  <  1:90904/49‑1
AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTG                                                  <  1:261165/49‑1
AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTG                                                  >  1:1405603/1‑49
 ATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAA                                                <  1:2060538/50‑1
 ATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAA                                                <  1:950141/50‑1
 ATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAA                                                <  1:1823132/50‑1
     AAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATT                                            <  1:1651000/50‑1
     AAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATAT                                             >  1:759046/1‑49
     AAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATAT                                             >  1:2258438/1‑49
      AACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTA                                           >  1:1704536/1‑50
      AACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTA                                           <  1:2013799/50‑1
         AGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCG                                        >  1:527297/1‑50
          GCAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGG                                       <  1:1708058/50‑1
          GCAATAATTTTATATTCACTGAAAATATTTTAAATCTTGAATATTATCGG                                       >  1:833417/1‑50
           CAATAATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGG                                      >  1:1184683/1‑50
                ATTTTATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTT                                  >  1:1466147/1‑49
                    TATATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAA                              >  1:813666/1‑49
                      TATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATAT                           <  1:1278072/50‑1
                      TATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATAT                           <  1:871054/50‑1
                      TATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATA                            <  1:549801/49‑1
                       ATTCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATAT                           >  1:1720946/1‑49
                        TTCACTGAAAATATTTTAAATCGTGAATATTATCGGGCGTTAAAATATT                          <  1:1186144/49‑1
                         TCACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTA                         >  1:870254/1‑49
                          CACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTAC                        <  1:380664/49‑1
                           ACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACAT                      <  1:1300480/50‑1
                           ACTGAAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACA                       <  1:701571/49‑1
                               AAAATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTG                  <  1:1940556/50‑1
                                  ATATTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTT               <  1:1182154/50‑1
                                     TTTTTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTT             >  1:1262517/1‑49
                                        TTAATCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAAT         >  1:2089061/1‑50
                                            TCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT     >  1:722589/1‑50
                                            TCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT     >  1:779127/1‑50
                                            TCTTGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGT     >  1:1716422/1‑50
                                               TGAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  1:1053075/1‑50
                                                                                                 
AATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/389850‑389897
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATATTATCGGGCGTTAAAATATTACATTGTGGTTTTCAATGCGTTAT  >  NC_000913/399806‑399854

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.