Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 167,484 Δ5,964 bp [fhuA]fhuB [fhuA], fhuC, fhuD, fhuB

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 167484 173447 5964 87 [0] [0] 89 [fhuA]–fhuB [fhuA],fhuC,fhuD,fhuB

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1674830 (0.000)86 (0.970) 60/98 0.0 100% intergenic (+219/‑1) mrcB/fhuA peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
?NC_000913 173448 = 0 (0.000)intergenic (+4/+154) fhuB/hemL iron(III) hydroxamate ABC transporter membrane subunit/glutamate‑1‑semialdehyde aminotransferase

GCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173401‑173497
                                               |                                                 
atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 <  1:4768015‑M2/3‑1 (MQ=255)
 tcattttcgtttacgttatcattcactttacatcagagatataccaTTc                                                 >  1:2334731‑M2/47‑49 (MQ=255)
 tcattctcgtttacgttatcattcactttacatcagagatataccaTTcc                                                <  1:3116526‑M2/4‑1 (MQ=255)
  cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:5436018‑M2/46‑50 (MQ=255)
  cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:1937188‑M2/46‑50 (MQ=255)
   attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                              <  1:6279451‑M2/6‑1 (MQ=255)
   attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:1530408‑M2/45‑49 (MQ=255)
   attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:5620281‑M2/45‑49 (MQ=255)
    ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             >  1:6468501‑M2/44‑50 (MQ=255)
    ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             <  1:1787824‑M2/7‑1 (MQ=255)
     tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             >  1:4371801‑M2/43‑49 (MQ=255)
      ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            <  1:240291‑M2/8‑1 (MQ=255)
      ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            <  1:5637385‑M2/8‑1 (MQ=255)
       tcgtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          >  1:1234301‑M2/41‑50 (MQ=255)
        cgtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         >  1:3026072‑M2/40‑50 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:4951468‑M2/9‑1 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:3262987‑M2/10‑1 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:384555‑M2/10‑1 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:4121099‑M2/11‑1 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:1641624‑M2/11‑1 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                        >  1:4811846‑M2/39‑50 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                        >  1:3140927‑M2/39‑50 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                        >  1:4096801‑M2/39‑50 (MQ=255)
         gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                        >  1:751393‑M2/39‑50 (MQ=255)
          tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:4618272‑M2/13‑1 (MQ=255)
          tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:6236263‑M2/13‑1 (MQ=255)
           ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      >  1:261649‑M2/37‑50 (MQ=255)
            tacgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                     >  1:5368906‑M2/36‑50 (MQ=255)
            tacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      <  1:171702‑M2/14‑1 (MQ=255)
             acgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                    >  1:3635760‑M2/35‑50 (MQ=255)
              cgttatcattcactttacatcagagatataccaTTCCATATAACggccat                                   <  1:1756825‑M2/17‑3 (MQ=255)
              cgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                    <  1:1225734‑M2/16‑1 (MQ=255)
              cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   >  1:4850019‑M2/34‑50 (MQ=255)
               gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:6192869‑M2/33‑50 (MQ=255)
               gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  <  1:3503051‑M2/18‑1 (MQ=255)
               gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:5921280‑M2/33‑50 (MQ=255)
               gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:2212647‑M2/33‑50 (MQ=255)
               gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:5438630‑M2/33‑50 (MQ=255)
               gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:290581‑M2/33‑50 (MQ=255)
                ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:4763004‑M2/32‑49 (MQ=255)
                ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:806314‑M2/32‑49 (MQ=255)
                ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:1060682‑M2/32‑49 (MQ=255)
                ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                 <  1:1748883‑M2/19‑1 (MQ=255)
                 tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                >  1:1115003‑M2/31‑50 (MQ=255)
                 tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                >  1:385876‑M2/31‑50 (MQ=255)
                   tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              <  1:1529002‑M2/22‑1 (MQ=255)
                   tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              <  1:5358575‑M2/22‑1 (MQ=255)
                   tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              <  1:1371049‑M2/22‑1 (MQ=255)
                   tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              <  1:492257‑M2/22‑1 (MQ=255)
                    cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             >  1:159528‑M2/28‑50 (MQ=255)
                     attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGcc                            <  1:839437‑M2/24‑1 (MQ=255)
                      ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           >  1:3192562‑M2/26‑50 (MQ=255)
                       tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:5259681‑M2/25‑50 (MQ=255)
                         actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTgg                         <  1:73020‑M2/27‑1 (MQ=255)
                         actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:5041545‑M2/23‑50 (MQ=255)
                           tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:6527705‑M2/21‑50 (MQ=255)
                            ttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:2669724‑M2/20‑49 (MQ=255)
                             tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACgg                     <  1:2985665‑M2/31‑1 (MQ=255)
                             tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:3038368‑M2/19‑50 (MQ=255)
                              acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    <  1:3938490‑M2/32‑1 (MQ=255)
                              acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:4469896‑M2/18‑49 (MQ=255)
                              acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:63588‑M2/18‑49 (MQ=255)
                               catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   >  1:6738434‑M2/17‑49 (MQ=255)
                                atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                 <  1:2335337‑M2/35‑1 (MQ=255)
                                atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:4372877‑M2/16‑49 (MQ=255)
                                 tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa                <  1:4533346‑M2/36‑1 (MQ=255)
                                 tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa                <  1:4325960‑M2/36‑1 (MQ=255)
                                 tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:788777‑M2/15‑48 (MQ=255)
                                 tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:330856‑M2/15‑48 (MQ=255)
                                  cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg               >  1:3955030‑M2/14‑50 (MQ=255)
                                   agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:1018700‑M2/38‑1 (MQ=255)
                                   agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:4179551‑M2/38‑1 (MQ=255)
                                     agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt            >  1:1809810‑M2/11‑50 (MQ=255)
                                      gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:4096135‑M2/10‑50 (MQ=255)
                                        tataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa         <  1:2485609‑M2/43‑1 (MQ=255)
                                         ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        >  1:6042292‑M2/7‑50 (MQ=255)
                                          taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcg       <  1:2287107‑M2/45‑1 (MQ=255)
                                          taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcg       >  1:5505813‑M2/6‑50 (MQ=255)
                                           accaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      >  1:5664293‑M2/5‑50 (MQ=255)
                                           accaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      >  1:1200224‑M2/5‑50 (MQ=255)
                                            ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgcg     <  1:3056564‑M2/47‑1 (MQ=255)
                                            ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      >  1:5499498‑M2/4‑49 (MQ=255)
                                            ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      <  1:554126‑M2/46‑1 (MQ=255)
                                            ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      >  1:4529739‑M2/4‑49 (MQ=255)
                                             caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:5188124‑M2/3‑50 (MQ=255)
                                              aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:924025‑M2/2‑49 (MQ=255)
                                               ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc   >  1:423439/1‑49 (MQ=255)
                                               ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc   >  1:5260410/1‑49 (MQ=255)
                                               ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:356188/50‑1 (MQ=255)
                                               |                                                 
GCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173401‑173497

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