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 85 [0] [0] 86 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)82 (1.180) 50/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:3935146‑M2/50‑50 (MQ=255)
aaaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                   >  1:1590181‑M2/50‑50 (MQ=255)
 aaggcaactggtgatttttattgattgctgagttttgaatacattgaaCg                                                  >  1:3297801‑M2/49‑50 (MQ=255)
   ggcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                                <  1:794045‑M2/4‑1 (MQ=255)
   ggcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                                <  1:1805867‑M2/4‑1 (MQ=255)
    gcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                                >  1:3558234‑M2/46‑49 (MQ=255)
    gcaactggtgatttttattgattgctgagttttgaatacattgaaCGTCt                                               <  1:2981670‑M2/5‑1 (MQ=255)
     caactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                              >  1:3243129‑M2/45‑50 (MQ=255)
      aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                              >  1:3225904‑M2/44‑49 (MQ=255)
       actggtgatttttattgattgctgagttttgaatacattgaaCGTCTaca                                            >  1:854314‑M2/43‑50 (MQ=255)
       actggtgatttttattgattgctgagttttgaatacattgaaCGTCTaca                                            <  1:254076‑M2/8‑1 (MQ=255)
       actggtgatttttattgattgctgagttttgaatacattgaaCGTCTaca                                            >  1:169501‑M2/43‑50 (MQ=255)
         tggtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                          >  1:5142691‑M2/41‑50 (MQ=255)
         tggtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                          >  1:672708‑M2/41‑50 (MQ=255)
          ggtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                          >  1:4946723‑M2/40‑49 (MQ=255)
           gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                          <  1:3429164‑M2/10‑1 (MQ=255)
           gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTa                                        >  1:468579‑M2/39‑50 (MQ=255)
              atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACt                                      <  1:1727660‑M2/14‑1 (MQ=255)
              atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTc                                     <  1:1605975‑M2/15‑1 (MQ=255)
               tttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTcc                                    <  1:5129308‑M2/16‑1 (MQ=255)
               tttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTcc                                    <  1:3897867‑M2/16‑1 (MQ=255)
                 tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                   <  1:4209250‑M2/17‑1 (MQ=255)
                 tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                  >  1:1607332‑M2/33‑50 (MQ=255)
                  ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                  >  1:1285429‑M2/32‑49 (MQ=255)
                  ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                  <  1:1850161‑M2/18‑1 (MQ=255)
                   tattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGt                                 >  1:2894173‑M2/31‑49 (MQ=255)
                    attgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTa                                >  1:2941996‑M2/30‑49 (MQ=255)
                    attgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTa                                >  1:4048979‑M2/30‑49 (MQ=255)
                     ttgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTaaa                              >  1:360831‑M2/29‑50 (MQ=255)
                     ttgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTaaa                              >  1:3716729‑M2/29‑50 (MQ=255)
                     ttgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTaaa                              >  1:3741419‑M2/29‑50 (MQ=255)
                       gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            >  1:4304500‑M2/27‑50 (MQ=255)
                       gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            >  1:130491‑M2/27‑50 (MQ=255)
                       gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                             <  1:1819889‑M2/23‑1 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            <  1:1716979‑M2/24‑1 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCg                           >  1:443914‑M2/26‑50 (MQ=255)
                         ttgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGc                          >  1:2503294‑M2/25‑50 (MQ=255)
                          tgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:3762562‑M2/24‑50 (MQ=255)
                          tgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:529641‑M2/24‑50 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:4527380‑M2/23‑49 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:1773779‑M2/23‑49 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:1309690‑M2/23‑49 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         >  1:3564708‑M2/23‑49 (MQ=255)
                            ctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCa                       <  1:4655960‑M2/29‑1 (MQ=255)
                              gagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:4327367‑M2/20‑49 (MQ=255)
                              gagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:2186934‑M2/20‑49 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:5354318‑M2/19‑48 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:976070‑M2/19‑48 (MQ=255)
                                gttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAg                   >  1:1550639‑M2/18‑50 (MQ=255)
                                gttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGAAg                   >  1:3176433‑M2/18‑50 (MQ=255)
                                 ttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAg                   >  1:1041069‑M2/17‑49 (MQ=255)
                                 ttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                  <  1:734753‑M2/34‑1 (MQ=255)
                                 ttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                  <  1:1113953‑M2/34‑1 (MQ=255)
                                 ttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                  <  1:5542156‑M2/34‑1 (MQ=255)
                                 ttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGc                  <  1:4256031‑M2/34‑1 (MQ=255)
                                  tttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAg                   >  1:4775656‑M2/16‑48 (MQ=255)
                                   ttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg                >  1:1832265‑M2/15‑50 (MQ=255)
                                    tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg                >  1:3207806‑M2/14‑49 (MQ=255)
                                    tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg                >  1:1052549‑M2/14‑49 (MQ=255)
                                     gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt               <  1:3958538‑M2/37‑1 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa             <  1:1917584‑M2/39‑1 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa             >  1:3520671‑M2/12‑50 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa             >  1:4795400‑M2/12‑50 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTa              <  1:149900‑M2/38‑1 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTa              <  1:2473429‑M2/38‑1 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTTAGGAGCAGTAAt            >  1:1793025‑M2/11‑50 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:2210619‑M2/11‑50 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:673476‑M2/11‑50 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:3022900‑M2/11‑50 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:5390637‑M2/11‑50 (MQ=255)
                                        tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:2880112‑M2/10‑49 (MQ=255)
                                        tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:5075532‑M2/10‑49 (MQ=255)
                                        tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATg           <  1:4661899‑M2/41‑1 (MQ=255)
                                           attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt        >  1:1705209‑M2/7‑50 (MQ=255)
                                           attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt        >  1:1235592‑M2/7‑50 (MQ=255)
                                            ttgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc       <  1:4698215‑M2/45‑1 (MQ=255)
                                            ttgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc       <  1:3096536‑M2/45‑1 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCt      <  1:1277347‑M2/46‑1 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCt      <  1:4545103‑M2/46‑1 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCt      <  1:10436‑M2/46‑1 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     >  1:948167‑M2/4‑50 (MQ=255)
                                              gaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     >  1:2056909‑M2/4‑50 (MQ=255)
                                               aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     <  1:5166254‑M2/47‑1 (MQ=255)
                                               aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTAc    <  1:5360350‑M2/48‑1 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  <  1:2097613/50‑1 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  <  1:382809/50‑1 (MQ=255)
                                                 |                                                 
AAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACAA  >  NC_000913/738914‑739012

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: