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 59 [0] [0] 60 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)56 (1.140) 45/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:1033014‑M2/50‑50 (MQ=255)
aaaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                   <  1:186367‑M2/1‑1 (MQ=255)
aaaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                   <  1:1827057‑M2/1‑1 (MQ=255)
aaaggcaactggtgatttttattgattgctgagttttgaatacattgaac                                                   <  1:176229‑M2/1‑1 (MQ=255)
 aaggcaactggtgatttttattgattgctgagttttgaatacattgaaCg                                                  <  1:46555‑M2/2‑1 (MQ=255)
 aaggcaactggtgatttttattgattgctgagttttgaatacattgaaCg                                                  <  1:3724549‑M2/2‑1 (MQ=255)
 aaggcaactggtgatttttattgattgctgagttttgaatacattgaaCg                                                  >  1:1845428‑M2/49‑50 (MQ=255)
   ggcaactggtgatttttattgattgctgagttttgaatacattgaaCGTc                                                >  1:3663729‑M2/47‑50 (MQ=255)
     caactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                              >  1:1063291‑M2/45‑50 (MQ=255)
     caactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                              <  1:2836324‑M2/6‑1 (MQ=255)
      aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTa                                              >  1:2753503‑M2/44‑49 (MQ=255)
       actggtgatttttattgattgctgagttttgaatacattgaaCGTCTaca                                            <  1:2447206‑M2/8‑1 (MQ=255)
           gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                          <  1:3096076‑M2/10‑1 (MQ=255)
           gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTa                                        >  1:296186‑M2/39‑50 (MQ=255)
           gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTa                                        <  1:1992147‑M2/12‑1 (MQ=255)
            tgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTAc                                       <  1:1841649‑M2/13‑1 (MQ=255)
            tgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTAc                                       >  1:2096080‑M2/38‑50 (MQ=255)
            tgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTAc                                       >  1:1864597‑M2/38‑50 (MQ=255)
             gatttttattgattgctgagttttgaatacattgaaCGTCTACACTTACt                                      >  1:802594‑M2/37‑50 (MQ=255)
              atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACt                                      <  1:2567851‑M2/14‑1 (MQ=255)
              atttttattgattgctgagttttgaatacattgaaCGTATACACTTACt                                      <  1:803314‑M2/14‑1 (MQ=255)
               tttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTcc                                    <  1:3406233‑M2/16‑1 (MQ=255)
                ttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                   >  1:1940783‑M2/34‑50 (MQ=255)
                 tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                   <  1:1964273‑M2/17‑1 (MQ=255)
                  ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                  >  1:1968164‑M2/32‑49 (MQ=255)
                     ttgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTaaa                              <  1:1987762‑M2/22‑1 (MQ=255)
                      tgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                             >  1:1170364‑M2/28‑50 (MQ=255)
                       gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            >  1:2725170‑M2/27‑50 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            <  1:2174887‑M2/24‑1 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            <  1:2250915‑M2/24‑1 (MQ=255)
                        attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                            <  1:430629‑M2/24‑1 (MQ=255)
                         ttgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGc                          <  1:326110‑M2/26‑1 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         <  1:1551875‑M2/27‑1 (MQ=255)
                           gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                         <  1:1788159‑M2/27‑1 (MQ=255)
                            ctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCa                       <  1:1846717‑M2/29‑1 (MQ=255)
                             tgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:1535726‑M2/21‑50 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                      >  1:3181309‑M2/19‑48 (MQ=255)
                               agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGa                    <  1:1945729‑M2/32‑1 (MQ=255)
                                gttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAg                   <  1:1058687‑M2/33‑1 (MQ=255)
                                gttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAg                   <  1:1994716‑M2/33‑1 (MQ=255)
                                   ttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCa                 <  1:3483223‑M2/35‑1 (MQ=255)
                                   ttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg                >  1:1290607‑M2/15‑50 (MQ=255)
                                    tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg                <  1:704167‑M2/36‑1 (MQ=255)
                                     gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt               <  1:1375775‑M2/37‑1 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa             >  1:2748647‑M2/12‑50 (MQ=255)
                                      aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa             <  1:3899367‑M2/39‑1 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            <  1:2615846‑M2/40‑1 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            <  1:1544185‑M2/40‑1 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:607303‑M2/11‑50 (MQ=255)
                                       atacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt            >  1:139843‑M2/11‑50 (MQ=255)
                                        tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATg           <  1:19046‑M2/41‑1 (MQ=255)
                                          cattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGaa         <  1:2766673‑M2/43‑1 (MQ=255)
                                          cattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGa          <  1:3598226‑M2/42‑1 (MQ=255)
                                           attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt        <  1:2487008‑M2/44‑1 (MQ=255)
                                           attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt        >  1:9014‑M2/7‑50 (MQ=255)
                                            ttgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc       >  1:1030512‑M2/6‑50 (MQ=255)
                                               aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     <  1:2477384‑M2/47‑1 (MQ=255)
                                               aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     <  1:2704398‑M2/47‑1 (MQ=255)
                                               aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa     <  1:947937‑M2/47‑1 (MQ=255)
                                                 cGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACaa  <  1:3913303/50‑1 (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: