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

ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167434‑167531
                                                 |                                                
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                  >  1:594753/1‑50 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               <  1:4768015‑M1/50‑4 (MQ=255)
    tCATTTTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:2334731‑M1/1‑46 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcc                                              <  1:3116526‑M1/50‑5 (MQ=255)
     cATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                             >  1:1937188‑M1/1‑45 (MQ=255)
     cATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                             >  1:5436018‑M1/1‑45 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                            <  1:6279451‑M1/50‑7 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                             >  1:1530408‑M1/1‑44 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                             >  1:5620281‑M1/1‑44 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           <  1:1787824‑M1/50‑8 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           >  1:6468501‑M1/1‑43 (MQ=255)
        tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           >  1:4371801‑M1/1‑42 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:240291‑M1/49‑9 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:5637385‑M1/49‑9 (MQ=255)
          tcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                        >  1:1234301‑M1/1‑40 (MQ=255)
           cGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       >  1:3026072‑M1/1‑39 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      >  1:3140927‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      >  1:4811846‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      >  1:751393‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      >  1:4096801‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       <  1:4121099‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       <  1:1641624‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                        <  1:384555‑M1/48‑11 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                        <  1:3262987‑M1/48‑11 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                         <  1:4951468‑M1/47‑10 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     <  1:4618272‑M1/50‑14 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     <  1:6236263‑M1/50‑14 (MQ=255)
              ttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                    >  1:261649‑M1/1‑36 (MQ=255)
               tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggcc                                   >  1:5368906‑M1/1‑35 (MQ=255)
               tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                    <  1:171702‑M1/49‑15 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccc                                  >  1:3635760‑M1/1‑34 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                 >  1:4850019‑M1/1‑33 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccc                                  <  1:1225734‑M1/49‑17 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccat                                 <  1:1756825‑M1/50‑18 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:2212647‑M1/1‑32 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                <  1:3503051‑M1/50‑19 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:6192869‑M1/1‑32 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:290581‑M1/1‑32 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:5438630‑M1/1‑32 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:5921280‑M1/1‑32 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                               <  1:1748883‑M1/50‑20 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:4763004‑M1/1‑31 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:806314‑M1/1‑31 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:1060682‑M1/1‑31 (MQ=255)
                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                              >  1:1115003‑M1/1‑30 (MQ=255)
                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                              >  1:385876‑M1/1‑30 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                            <  1:1371049‑M1/50‑23 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                            <  1:1529002‑M1/50‑23 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                            <  1:492257‑M1/50‑23 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                            <  1:5358575‑M1/50‑23 (MQ=255)
                       cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                           >  1:159528‑M1/1‑27 (MQ=255)
                        aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                          <  1:839437‑M1/50‑25 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                         >  1:3192562‑M1/1‑25 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        >  1:5259681‑M1/1‑24 (MQ=255)
                            aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                      >  1:5041545‑M1/1‑22 (MQ=255)
                            aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgg                       <  1:73020‑M1/49‑28 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    >  1:6527705‑M1/1‑20 (MQ=255)
                               ttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    >  1:2669724‑M1/1‑19 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:3038368‑M1/1‑18 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgg                   <  1:2985665‑M1/49‑32 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:63588‑M1/1‑17 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:4469896‑M1/1‑17 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  <  1:3938490‑M1/49‑33 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                 >  1:6738434‑M1/1‑16 (MQ=255)
                                   aTCAGAGATATACCAttccatataacggccctgtacgcctggacggataa               <  1:2335337‑M1/50‑36 (MQ=255)
                                   aTCAGAGATATACCAttccatataacggccctgtacgcctggacggata                >  1:4372877‑M1/1‑15 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa              <  1:4533346‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa              <  1:4325960‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                >  1:788777‑M1/1‑14 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                >  1:330856‑M1/1‑14 (MQ=255)
                                     cAGAGATATACCAttccatataacggccctgtacgcctggacggataaag             >  1:3955030‑M1/1‑13 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc            <  1:4179551‑M1/50‑39 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc            <  1:1018700‑M1/50‑39 (MQ=255)
                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt          >  1:1809810‑M1/1‑10 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt         >  1:4096135‑M1/1‑9 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta       <  1:2485609‑M1/50‑44 (MQ=255)
                                            ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac      >  1:6042292‑M1/1‑6 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg     <  1:2287107‑M1/50‑46 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg     >  1:5505813‑M1/1‑5 (MQ=255)
                                              aCCAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:1200224‑M1/1‑4 (MQ=255)
                                              aCCAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:5664293‑M1/1‑4 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgcg   <  1:3056564‑M1/50‑48 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:5499498‑M1/1‑3 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc    <  1:554126‑M1/49‑47 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:4529739‑M1/1‑3 (MQ=255)
                                                cAttccatataacggccctgtacgcctggacggataaagcgtttacgcga  >  1:5188124‑M1/1‑2 (MQ=255)
                                                 attccatataacggccctgtacgcctggacggataaagcgtttacgcga  >  1:924025‑M1/1‑1 (MQ=255)
                                                 |                                                
ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167434‑167531

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