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 152 [0] [0] 148 [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)146 (1.350) 78/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

Only 100 of 148 total aligned reads displayed.

CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173399‑173497
                                                 |                                                 
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                   >  1:2675686‑M2/50‑50 (MQ=255)
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                   <  1:207655‑M2/1‑1 (MQ=255)
 aatcattctcgtttacgttatcattcactttacatcagagatataccatt                                                  >  1:1531033‑M2/49‑50 (MQ=255)
   tcattctcgtttacgttatcattcactttacatcagagatataccaTTcc                                                <  1:4354398‑M2/4‑1 (MQ=255)
   tcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 >  1:1822697‑M2/47‑49 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               <  1:341419‑M2/5‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:5788123‑M2/45‑49 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:2174412‑M2/45‑49 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               <  1:1923821‑M2/5‑1 (MQ=255)
      ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             >  1:1559182‑M2/44‑50 (MQ=255)
        ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            <  1:5182498‑M2/8‑1 (MQ=255)
        ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            <  1:2962283‑M2/8‑1 (MQ=255)
         tcgtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          >  1:5018696‑M2/41‑50 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:5400616‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:4398824‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:4857258‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:4773704‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:2867581‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:1043064‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:3722760‑M2/11‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:1353260‑M2/11‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:1221497‑M2/13‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:7580326‑M2/13‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:8164464‑M2/13‑1 (MQ=255)
             ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      >  1:7114214‑M2/37‑50 (MQ=255)
             ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      >  1:1919404‑M2/37‑50 (MQ=255)
              tacgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                     >  1:7488070‑M2/36‑50 (MQ=255)
               acgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                     <  1:2194527‑M2/15‑1 (MQ=255)
               acgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                     <  1:4137741‑M2/15‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                    <  1:5738846‑M2/16‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   >  1:2679753‑M2/34‑50 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   <  1:6468597‑M2/17‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                   <  1:1709177‑M2/17‑1 (MQ=255)
                 gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:6754898‑M2/33‑50 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                 <  1:3717046‑M2/19‑1 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                 <  1:3936061‑M2/19‑1 (MQ=255)
                   tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                >  1:8017759‑M2/31‑50 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                <  1:7783314‑M2/20‑1 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                <  1:5155852‑M2/20‑1 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               >  1:4686136‑M2/30‑50 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              >  1:3541699‑M2/29‑50 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              <  1:6042523‑M2/22‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATAGAACGGCCCTGTACg                              <  1:6700339‑M2/22‑1 (MQ=255)
                      cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             >  1:4450750‑M2/28‑50 (MQ=255)
                       attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGcc                            <  1:2940720‑M2/24‑1 (MQ=255)
                       attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGcc                            >  1:6907988‑M2/27‑50 (MQ=255)
                       attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             <  1:4176033‑M2/23‑1 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           >  1:751772‑M2/26‑50 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           <  1:7996186‑M2/25‑1 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           <  1:5567926‑M2/25‑1 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:7585188‑M2/25‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:2047855‑M2/25‑50 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:4637380‑M2/24‑49 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          <  1:466811‑M2/26‑1 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          <  1:3641153‑M2/26‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTgg                         <  1:3335931‑M2/27‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTgg                         <  1:8243316‑M2/27‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:3137689‑M2/23‑50 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:3227965‑M2/23‑50 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                       <  1:7677263‑M2/29‑1 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                       <  1:3044118‑M2/29‑1 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                       >  1:4002343‑M2/22‑50 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                       <  1:5333152‑M2/29‑1 (MQ=255)
                             tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:1060601‑M2/21‑50 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACgg                     <  1:6748005‑M2/31‑1 (MQ=255)
                                acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    <  1:1378993‑M2/32‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   <  1:1510383‑M2/33‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  <  1:8066580‑M2/34‑1 (MQ=255)
                                  atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                 <  1:6594900‑M2/35‑1 (MQ=255)
                                  atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:1483608‑M2/16‑49 (MQ=255)
                                   tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:568461‑M2/15‑48 (MQ=255)
                                   tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:668274‑M2/15‑48 (MQ=255)
                                    cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg               >  1:6856422‑M2/14‑50 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg               <  1:6573767‑M2/37‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:4651498‑M2/38‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:7966562‑M2/38‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              >  1:3079031‑M2/13‑50 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:3500098‑M2/38‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:4295479‑M2/38‑1 (MQ=255)
                                       agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt            >  1:5044118‑M2/11‑50 (MQ=255)
                                        gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:7106138‑M2/10‑50 (MQ=255)
                                        gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           <  1:3451374‑M2/41‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt          <  1:459991‑M2/42‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt          <  1:2021024‑M2/42‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:5572177‑M2/9‑49 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:2988707‑M2/9‑49 (MQ=255)
                                          tataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa         <  1:5663410‑M2/43‑1 (MQ=255)
                                           ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        >  1:85155‑M2/7‑50 (MQ=255)
                                           ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        >  1:6135415‑M2/7‑50 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcg       >  1:1387292‑M2/6‑50 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        <  1:1061175‑M2/44‑1 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        <  1:27853‑M2/44‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgcg     <  1:2403609‑M2/47‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      >  1:3322045‑M2/4‑49 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      <  1:5498315‑M2/46‑1 (MQ=255)
                                               caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:7971913‑M2/3‑50 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:3617622‑M2/2‑49 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa    >  1:5289564‑M2/2‑49 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc   <  1:6206097‑M2/49‑1 (MQ=255)
                                                 ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:626618/50‑1 (MQ=255)
                                                 |                                                 
CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173399‑173497

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: 
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.