Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 729,365 Δ9,598 bp rhsCybfD rhsC, ybfB, ybfO, ybfC, ybfQ, ybfL, ybfD

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 729365 738962 9598 49 [0] [0] 50 rhsC–ybfD rhsC,ybfB,ybfO,ybfC,ybfQ,ybfL,ybfD

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 7293640 (0.000)46 (1.300) 37/92 0.0 100% intergenic (+24/‑219) ybfA/rhsC DUF2517 domain‑containing protein YbfA/rhs element protein RhsC
?NC_000913 738963 = 0 (0.000)intergenic (+110/‑38) ybfD/ybgA H repeat‑associated putative transposase YbfD/DUF1722 domain‑containing protein YbgA

AATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTAC  >  NC_000913/738915‑739010
                                                |                                               
aaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                 >  1:1770693‑M2/49‑49 (MQ=255)
aaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                 <  1:1508005‑M2/1‑1 (MQ=255)
 aggcaactggtgatttttattgattgctgagttttgaatacattgaaCg                                                <  1:812557‑M2/2‑1 (MQ=255)
  ggcaactggtgatttttattgattgctgagttttgaatacattgaaCGt                                               <  1:1503679‑M2/3‑1 (MQ=255)
   gcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                              >  1:121126‑M2/46‑49 (MQ=255)
   gcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                              <  1:519182‑M2/4‑1 (MQ=255)
   gcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                              <  1:1329966‑M2/4‑1 (MQ=255)
    caactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                              >  1:562728‑M2/45‑48 (MQ=255)
    caactggtgatttttattgattgctgagttttgaatacattgaaCGTCt                                             <  1:761033‑M2/5‑1 (MQ=255)
     aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                            <  1:318727‑M2/6‑1 (MQ=255)
      actggtgatttttattgattgctgagttttgaatacattgaaCGTCTac                                           <  1:2415771‑M2/7‑1 (MQ=255)
          gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                        >  1:1322685‑M2/39‑48 (MQ=255)
          gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                        <  1:1629999‑M2/10‑1 (MQ=255)
          gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                        <  1:572666‑M2/10‑1 (MQ=255)
           tgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTa                                      >  1:2254372‑M2/38‑49 (MQ=255)
             atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACt                                    >  1:138610‑M2/36‑49 (MQ=255)
              tttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTc                                   <  1:2451407‑M2/15‑1 (MQ=255)
                tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTcc                                  <  1:1693499‑M2/16‑1 (MQ=255)
                tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                 >  1:1387922‑M2/33‑49 (MQ=255)
                 ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                >  1:555056‑M2/32‑49 (MQ=255)
                 ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                >  1:1377445‑M2/32‑49 (MQ=255)
                     tgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTaaa                            >  1:1299034‑M2/28‑49 (MQ=255)
                      gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                           >  1:1595296‑M2/27‑49 (MQ=255)
                      gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                           >  1:1057720‑M2/27‑49 (MQ=255)
                      gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                           >  1:417708‑M2/27‑49 (MQ=255)
                       attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                          <  1:2019996‑M2/24‑1 (MQ=255)
                            tgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                       >  1:613845‑M2/21‑47 (MQ=255)
                            tgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                       >  1:341492‑M2/21‑47 (MQ=255)
                             gagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                    >  1:153102‑M2/20‑49 (MQ=255)
                              agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCa                     <  1:543978‑M2/29‑1 (MQ=255)
                              agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAgg                   <  1:521865‑M2/31‑1 (MQ=255)
                              agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAgg                   <  1:2623761‑M2/31‑1 (MQ=255)
                              agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                    <  1:1375‑M2/30‑1 (MQ=255)
                                  ttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAg                 >  1:1472792‑M2/15‑47 (MQ=255)
                                  ttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                <  1:2090300‑M2/34‑1 (MQ=255)
                                   tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg              >  1:2058173‑M2/14‑49 (MQ=255)
                                   tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg              >  1:1292584‑M2/14‑49 (MQ=255)
                                    gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg              >  1:2369837‑M2/13‑48 (MQ=255)
                                     aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt             <  1:2405849‑M2/37‑1 (MQ=255)
                                     aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTa            <  1:433774‑M2/38‑1 (MQ=255)
                                       tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt          >  1:2314971‑M2/10‑49 (MQ=255)
                                          attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGaa       >  1:342464‑M2/7‑49 (MQ=255)
                                           ttgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt      >  1:264643‑M2/6‑49 (MQ=255)
                                           ttgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt      <  1:1890853‑M2/44‑1 (MQ=255)
                                             gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc     >  1:1552389‑M2/4‑48 (MQ=255)
                                             gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc     >  1:1445268‑M2/4‑48 (MQ=255)
                                              aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa   >  1:2381779‑M2/3‑49 (MQ=255)
                                              aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa   >  1:1983011‑M2/3‑49 (MQ=255)
                                              aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa   >  1:1863830‑M2/3‑49 (MQ=255)
                                                cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTAc  >  1:830703/1‑48 (MQ=255)
                                                |                                               
AATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTAC  >  NC_000913/738915‑739010

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