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 64 [0] [0] 65 [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)64 (1.100) 48/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

CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173399‑173497
                                                 |                                                 
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                   >  1:2326662‑M2/50‑50 (MQ=255)
 aatcattctcgtttacgttatcattcactttacatcagagatataccatt                                                  >  1:4354446‑M2/49‑50 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccatt                                                  <  1:518325‑M2/2‑1 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 >  1:3303512‑M2/48‑50 (MQ=255)
   tcattctcgtttacgttatcattcactttacatcagagatataccaTTcc                                                <  1:2427534‑M2/4‑1 (MQ=255)
   tcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 >  1:416501‑M2/47‑49 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                              <  1:754369‑M2/6‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:2256534‑M2/45‑49 (MQ=255)
      ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             <  1:4172418‑M2/7‑1 (MQ=255)
      ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             <  1:602018‑M2/7‑1 (MQ=255)
       tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             >  1:1284605‑M2/43‑49 (MQ=255)
        ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            >  1:908936‑M2/42‑49 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:2581213‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:3878154‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:4258308‑M2/11‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:282205‑M2/13‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:2363588‑M2/13‑1 (MQ=255)
             ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      >  1:3545715‑M2/37‑50 (MQ=255)
              tacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      <  1:1689121‑M2/14‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   >  1:2517136‑M2/34‑50 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   <  1:2808297‑M2/17‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   <  1:2587017‑M2/17‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   >  1:4350371‑M2/34‑50 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   >  1:1771974‑M2/34‑50 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   >  1:172888‑M2/34‑50 (MQ=255)
                 gttatcattcactttacatcagagatataccaTTCCATATAACGGCTCTg                                  >  1:780816‑M2/33‑50 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                 <  1:2963731‑M2/19‑1 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                <  1:2531060‑M2/20‑1 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                <  1:1371678‑M2/20‑1 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               >  1:2758392‑M2/30‑50 (MQ=255)
                      cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             >  1:983164‑M2/28‑50 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           >  1:2542078‑M2/26‑50 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATACAACGGCCCTGTACGCCt                           <  1:3230889‑M2/25‑1 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          <  1:1020898‑M2/26‑1 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:3704520‑M2/24‑49 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:2308068‑M2/23‑50 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:3995312‑M2/23‑50 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:3632921‑M2/23‑50 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:1617847‑M2/23‑50 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:3264027‑M2/23‑50 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:3013906‑M2/23‑50 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                       <  1:1540443‑M2/29‑1 (MQ=255)
                             tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:4214041‑M2/21‑50 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      <  1:640433‑M2/30‑1 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:3219914‑M2/19‑50 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   <  1:1458291‑M2/33‑1 (MQ=255)
                                  atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:3632277‑M2/16‑49 (MQ=255)
                                  atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:2615974‑M2/16‑49 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              >  1:2629194‑M2/13‑50 (MQ=255)
                                       agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt            >  1:3425758‑M2/11‑50 (MQ=255)
                                        gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:2914531‑M2/10‑50 (MQ=255)
                                        gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           <  1:1462112‑M2/41‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt          <  1:2118341‑M2/42‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt          <  1:1309924‑M2/42‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:3308765‑M2/9‑49 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:2450591‑M2/9‑49 (MQ=255)
                                           ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        >  1:249053‑M2/7‑50 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        <  1:3332571‑M2/44‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      <  1:3525574‑M2/46‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      >  1:3404279‑M2/4‑49 (MQ=255)
                                               caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    <  1:4312295‑M2/48‑1 (MQ=255)
                                               caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:1716493‑M2/3‑50 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:1595134‑M2/2‑49 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc   <  1:1239746‑M2/49‑1 (MQ=255)
                                                 ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:857679/50‑1 (MQ=255)
                                                 |                                                 
CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173399‑173497

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