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 113 [0] [0] 113 [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)110 (1.510) 69/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 113 total aligned reads displayed.

CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173399‑173497
                                                 |                                                 
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                   <  1:3011230‑M2/1‑1 (MQ=255)
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                   <  1:5273522‑M2/1‑1 (MQ=255)
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                   >  1:4372660‑M2/50‑50 (MQ=255)
 aatcattctcgtttacgttatcattcactttacatcagagatataccatt                                                  >  1:5170171‑M2/49‑50 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccatt                                                  <  1:4132391‑M2/2‑1 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccatt                                                  <  1:2024161‑M2/2‑1 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 <  1:219308‑M2/3‑1 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 <  1:1707606‑M2/3‑1 (MQ=255)
   tcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                 >  1:5500322‑M2/47‑49 (MQ=255)
    cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               >  1:1573091‑M2/46‑50 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                              <  1:4229602‑M2/6‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                               <  1:3688135‑M2/5‑1 (MQ=255)
      ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                             <  1:4144104‑M2/7‑1 (MQ=255)
        ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            <  1:1473173‑M2/8‑1 (MQ=255)
        ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                            <  1:1255410‑M2/8‑1 (MQ=255)
         tcgtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          >  1:3481580‑M2/41‑50 (MQ=255)
          cgtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         >  1:4324047‑M2/40‑50 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:2479374‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                           <  1:75249‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:862453‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:3562800‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:3484710‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                          <  1:2132791‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                         <  1:2583365‑M2/11‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                        <  1:1604640‑M2/12‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:4843874‑M2/13‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                       <  1:4393878‑M2/13‑1 (MQ=255)
              tacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                      <  1:3400158‑M2/14‑1 (MQ=255)
               acgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                    >  1:2054650‑M2/35‑50 (MQ=255)
               acgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                     <  1:3369551‑M2/15‑1 (MQ=255)
               acgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                     <  1:4458556‑M2/15‑1 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGgccca                                   <  1:3070769‑M2/17‑2 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                    <  1:5174857‑M2/16‑1 (MQ=255)
                 gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:5360421‑M2/33‑50 (MQ=255)
                 gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  <  1:1386422‑M2/18‑1 (MQ=255)
                 gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  <  1:1276466‑M2/18‑1 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                  >  1:2673571‑M2/32‑49 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                 <  1:5232726‑M2/19‑1 (MQ=255)
                   tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                >  1:3787089‑M2/31‑50 (MQ=255)
                   tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                >  1:5117182‑M2/31‑50 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                <  1:2431827‑M2/20‑1 (MQ=255)
                    atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                                <  1:3523582‑M2/20‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               <  1:5130667‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               <  1:267558‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                               <  1:5175177‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                              <  1:4874575‑M2/22‑1 (MQ=255)
                      cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             >  1:2802477‑M2/28‑50 (MQ=255)
                      cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             >  1:4890039‑M2/28‑50 (MQ=255)
                      cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             >  1:843490‑M2/28‑50 (MQ=255)
                       attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGcc                            >  1:5506133‑M2/27‑50 (MQ=255)
                       attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                             <  1:2833261‑M2/23‑1 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           <  1:3622549‑M2/25‑1 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                           >  1:5351621‑M2/26‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:1127502‑M2/25‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:453691‑M2/25‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:1598301‑M2/25‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:2035694‑M2/25‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:2286345‑M2/25‑50 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                          >  1:5188277‑M2/24‑49 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTgg                         <  1:1894734‑M2/27‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                        >  1:2099717‑M2/23‑50 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCTTGTACGCCTGGAc                       <  1:2691549‑M2/29‑1 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                       >  1:931826‑M2/22‑50 (MQ=255)
                             tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:4363084‑M2/21‑50 (MQ=255)
                             tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:647719‑M2/21‑50 (MQ=255)
                             tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:845343‑M2/21‑50 (MQ=255)
                              ttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      >  1:4148383‑M2/20‑49 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                      <  1:4427677‑M2/30‑1 (MQ=255)
                                acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:377605‑M2/18‑49 (MQ=255)
                                acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                    >  1:2729854‑M2/18‑49 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   <  1:5387219‑M2/33‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   <  1:3728163‑M2/33‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   <  1:1230493‑M2/33‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   >  1:3966844‑M2/17‑49 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                   <  1:3551856‑M2/33‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  <  1:3801345‑M2/34‑1 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  <  1:4740203‑M2/34‑1 (MQ=255)
                                   tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa                <  1:4195997‑M2/36‑1 (MQ=255)
                                   tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                  >  1:720626‑M2/15‑48 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg               <  1:1077399‑M2/37‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:1086409‑M2/38‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              <  1:5429674‑M2/38‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc              >  1:5398403‑M2/13‑50 (MQ=255)
                                      gagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCg             <  1:546105‑M2/39‑1 (MQ=255)
                                       agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt            <  1:3857659‑M2/40‑1 (MQ=255)
                                       agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt            <  1:3681323‑M2/40‑1 (MQ=255)
                                       agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt            >  1:3508175‑M2/11‑50 (MQ=255)
                                        gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           <  1:260675‑M2/41‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt          <  1:4690558‑M2/42‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt           >  1:2248943‑M2/9‑49 (MQ=255)
                                           ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        >  1:159384‑M2/7‑50 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        <  1:1643162‑M2/44‑1 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        <  1:4079234‑M2/44‑1 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc        <  1:4350281‑M2/44‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      <  1:5197259‑M2/46‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc      <  1:1525619‑M2/46‑1 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGTCGGATAAAGCGTTTACGCGa    >  1:4865842‑M2/2‑49 (MQ=255)
                                                aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc   <  1:2290302‑M2/49‑1 (MQ=255)
                                                 ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:1002674/50‑1 (MQ=255)
                                                 ttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACa  <  1:15348/50‑1 (MQ=255)
                                                 |                                                 
CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACA  >  NC_000913/173399‑173497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: