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 46 [0] [0] 49 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)43 (1.200) 36/94 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

AAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAA  >  NC_000913/738914‑739012
                                                 |                                                 
aaaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                   <  1:3209‑M2/1‑1 (MQ=255)
  aggcaactggtgatttttattgattgctgagttttgaatacattgaaCGt                                                 >  1:235091‑M2/48‑50 (MQ=255)
  aggcaactggtgatttttattgattgctgagttttgaatacattgaaCGt                                                 <  1:1422795‑M2/3‑1 (MQ=255)
   ggcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                                >  1:949859‑M2/47‑50 (MQ=255)
      aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTac                                             <  1:2627906‑M2/7‑1 (MQ=255)
      aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTac                                             <  1:9350‑M2/7‑1 (MQ=255)
      aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                              >  1:2684662‑M2/44‑49 (MQ=255)
          ggtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                          >  1:440659‑M2/40‑49 (MQ=255)
           gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTa                                        <  1:2019460‑M2/12‑1 (MQ=255)
              atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACt                                      <  1:1821979‑M2/14‑1 (MQ=255)
              atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTc                                     <  1:2553967‑M2/15‑1 (MQ=255)
               tttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTcc                                    <  1:2153829‑M2/16‑1 (MQ=255)
                ttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                   >  1:2771215‑M2/34‑50 (MQ=255)
                 tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                   <  1:1938749‑M2/17‑1 (MQ=255)
                  ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                  >  1:60169‑M2/32‑49 (MQ=255)
                   tattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGt                                 >  1:2227030‑M2/31‑49 (MQ=255)
                    gttgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTa                                >  1:1132121‑M2/30‑49 (MQ=255)
                    attgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTa                                >  1:1132159‑M2/30‑49 (MQ=255)
                     ttgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTaaa                              >  1:330268‑M2/29‑50 (MQ=255)
                      tgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                             >  1:1159064‑M2/28‑50 (MQ=255)
                       gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            >  1:483336‑M2/27‑50 (MQ=255)
                       gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            >  1:746476‑M2/27‑50 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            <  1:1096557‑M2/24‑1 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCg                           <  1:344108‑M2/25‑1 (MQ=255)
                          tgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:208693‑M2/24‑50 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:1277250‑M2/23‑49 (MQ=255)
                            ctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTc                        <  1:2274860‑M2/28‑1 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:2204364‑M2/19‑48 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:2532660‑M2/19‑48 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:1443197‑M2/19‑48 (MQ=255)
                                 ttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                  <  1:2034891‑M2/34‑1 (MQ=255)
                                    tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg                >  1:1490835‑M2/14‑49 (MQ=255)
                                     gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt               <  1:2335080‑M2/37‑1 (MQ=255)
                                     gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt               <  1:2710346‑M2/37‑1 (MQ=255)
                                     gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt               >  1:860857‑M2/13‑49 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa             <  1:1378613‑M2/39‑1 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            <  1:793991‑M2/40‑1 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:411499‑M2/11‑50 (MQ=255)
                                          cattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGaa         <  1:2792608‑M2/43‑1 (MQ=255)
                                           attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt        >  1:2195235‑M2/7‑50 (MQ=255)
                                            ttgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc       >  1:39883‑M2/6‑50 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     >  1:2328192‑M2/4‑50 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     >  1:2618018‑M2/4‑50 (MQ=255)
                                               aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     <  1:1039526‑M2/47‑1 (MQ=255)
                                                aCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTAc    >  1:102902‑M2/2‑49 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  >  1:2349670/1‑50 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  >  1:2001709/1‑50 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  >  1:1792811/1‑50 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  >  1:1551969/1‑50 (MQ=255)
                                                 |                                                 
AAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAA  >  NC_000913/738914‑739012

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: