Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,618,992 Δ6,688 bp rhsB[yhhJ] rhsB, yhhH, yhhI, [yhhJ]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3618992 3625679 6688 98 [0] [0] 98 rhsB–[yhhJ] rhsB,yhhH,yhhI,[yhhJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 36189910 (0.000)96 (1.000) 54/96 0.0 100% intergenic (+2/‑201) nikR/rhsB DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
?NC_000913 3625680 = NA (NA)coding (1124/1125 nt) yhhJ ABC transporter family protein YhhJ

TTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625631‑3625728
                                                 |                                                
gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtt                                                  <  1:1382025‑M2/1‑1 (MQ=255)
gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtt                                                  <  1:3584954‑M2/1‑1 (MQ=255)
 cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 <  1:5739958‑M2/2‑1 (MQ=255)
 cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 <  1:6651514‑M2/2‑1 (MQ=255)
  ggcgtgcggcacgggcattcgcagtgcttgccgaaggaagatcgagtTAc                                                >  1:3153899‑M2/48‑50 (MQ=255)
   gcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACg                                               <  1:957899‑M2/4‑1 (MQ=255)
   gcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACg                                               <  1:5787138‑M2/4‑1 (MQ=255)
    cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              <  1:3127678‑M2/5‑1 (MQ=255)
    cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:1960890‑M2/46‑50 (MQ=255)
      tgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCa                                            <  1:2205404‑M2/7‑1 (MQ=255)
      tgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCa                                            >  1:3626251‑M2/44‑50 (MQ=255)
         ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         >  1:5672588‑M2/41‑50 (MQ=255)
         ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         >  1:4879190‑M2/41‑50 (MQ=255)
            acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:3885345‑M2/13‑1 (MQ=255)
            acgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTcc                                      <  1:1398999‑M2/13‑1 (MQ=255)
             cgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     <  1:6658675‑M2/14‑1 (MQ=255)
              gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     >  1:1720827‑M2/36‑49 (MQ=255)
              gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     >  1:924121‑M2/36‑49 (MQ=255)
              gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     >  1:2967439‑M2/36‑49 (MQ=255)
               ggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCaa                                   <  1:2587518‑M2/16‑1 (MQ=255)
               ggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCaa                                   >  1:2660885‑M2/35‑50 (MQ=255)
                gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  >  1:3326072‑M2/34‑50 (MQ=255)
                gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  >  1:1105860‑M2/34‑50 (MQ=255)
                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 >  1:4746283‑M2/33‑50 (MQ=255)
                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:1409834‑M2/18‑1 (MQ=255)
                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:3372942‑M2/18‑1 (MQ=255)
                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:5740650‑M2/32‑50 (MQ=255)
                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:6365253‑M2/32‑50 (MQ=255)
                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                >  1:2511106‑M2/32‑50 (MQ=255)
                   tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               <  1:5176293‑M2/20‑1 (MQ=255)
                   tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               >  1:4897427‑M2/31‑50 (MQ=255)
                   tttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGt                               >  1:7291964‑M2/31‑50 (MQ=255)
                     tgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTc                              <  1:6678105‑M2/21‑1 (MQ=255)
                     tgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             >  1:1422193‑M2/29‑50 (MQ=255)
                     tgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             >  1:6291947‑M2/29‑50 (MQ=255)
                      gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCtt                            >  1:6306388‑M2/28‑50 (MQ=255)
                      gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:2828335‑M2/22‑1 (MQ=255)
                      gcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             <  1:3231805‑M2/22‑1 (MQ=255)
                       cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCtt                            <  1:2024382‑M2/23‑1 (MQ=255)
                       cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCtt                            <  1:3284709‑M2/23‑1 (MQ=255)
                       cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCtt                            <  1:4015366‑M2/23‑1 (MQ=255)
                       cagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTc                           >  1:1847795‑M2/27‑50 (MQ=255)
                        agtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTcc                          <  1:926555‑M2/25‑1 (MQ=255)
                        agtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTcc                          <  1:394405‑M2/25‑1 (MQ=255)
                         gtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCt                         <  1:1771078‑M2/26‑1 (MQ=255)
                         gtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCt                         >  1:5876240‑M2/25‑50 (MQ=255)
                         gtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCt                         <  1:1901435‑M2/26‑1 (MQ=255)
                          tgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTg                        <  1:1610975‑M2/27‑1 (MQ=255)
                          tgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTg                        <  1:6850805‑M2/27‑1 (MQ=255)
                          tgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTg                        <  1:3280029‑M2/27‑1 (MQ=255)
                            cttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGaa                      <  1:5367246‑M2/29‑1 (MQ=255)
                             ttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAAt                     >  1:29717‑M2/21‑50 (MQ=255)
                               gccgaaggaagattgagtTACGTCATTGTCCCAATCGTCTTCCTGAATc                    <  1:6538710‑M2/31‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAAt                     <  1:3094234‑M2/30‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:5165210‑M2/31‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:5023540‑M2/31‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:4481314‑M2/31‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    <  1:735266‑M2/31‑1 (MQ=255)
                               gccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:2894387‑M2/19‑50 (MQ=255)
                                ccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCg                   >  1:1738716‑M2/18‑49 (MQ=255)
                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:5725437‑M2/34‑1 (MQ=255)
                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:1272420‑M2/34‑1 (MQ=255)
                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:781894‑M2/34‑1 (MQ=255)
                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:3079571‑M2/34‑1 (MQ=255)
                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:3549107‑M2/34‑1 (MQ=255)
                                  gaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                <  1:3353942‑M2/35‑1 (MQ=255)
                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagc               <  1:3032281‑M2/36‑1 (MQ=255)
                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:7422162‑M2/15‑49 (MQ=255)
                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:4402394‑M2/15‑49 (MQ=255)
                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:7086374‑M2/14‑48 (MQ=255)
                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:844508‑M2/14‑48 (MQ=255)
                                     ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:7383827‑M2/13‑50 (MQ=255)
                                     ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:6113564‑M2/13‑50 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            <  1:3390574‑M2/39‑1 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            >  1:7446703‑M2/12‑50 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            >  1:3052954‑M2/12‑50 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            <  1:2610611‑M2/39‑1 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            >  1:3793537‑M2/12‑50 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            <  1:4311896‑M2/39‑1 (MQ=255)
                                      gaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            <  1:5477733‑M2/39‑1 (MQ=255)
                                      gaagattgagaTACGCCATTGTCCCAATCGTCTTCCTGAATCGCagcagc            >  1:3151847‑M2/12‑50 (MQ=255)
                                       aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           <  1:7297295‑M2/40‑1 (MQ=255)
                                       aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           <  1:7348488‑M2/40‑1 (MQ=255)
                                       aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           >  1:5136586‑M2/11‑50 (MQ=255)
                                          attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        <  1:3607446‑M2/43‑1 (MQ=255)
                                          attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:2012968‑M2/8‑50 (MQ=255)
                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       >  1:6386926‑M2/7‑50 (MQ=255)
                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       >  1:385792‑M2/7‑50 (MQ=255)
                                           ttgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAAt       >  1:2911926‑M2/7‑50 (MQ=255)
                                            tgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      >  1:1338426‑M2/6‑50 (MQ=255)
                                             gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATgg     >  1:1286170‑M2/5‑50 (MQ=255)
                                              agtTACGCCGTTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:6818509‑M2/47‑1 (MQ=255)
                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATgg     <  1:7042863‑M2/46‑1 (MQ=255)
                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:3947349‑M2/47‑1 (MQ=255)
                                              agtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    <  1:85017‑M2/47‑1 (MQ=255)
                                               gtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGt    >  1:4974381‑M2/3‑49 (MQ=255)
                                                 tACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  >  1:4762414/1‑49 (MQ=255)
                                                 tACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  >  1:4904311/1‑49 (MQ=255)
                                                 |                                                
TTGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625631‑3625728

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