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.

ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167434‑167532
                                                 |                                                 
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                   >  1:1220510/1‑50 (MQ=255)
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                   >  1:2692847/1‑50 (MQ=255)
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                   <  1:2066923/50‑1 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  <  1:5273522‑M1/50‑2 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  <  1:3011230‑M1/50‑2 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  >  1:4372660‑M1/1‑49 (MQ=255)
  aaTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 >  1:5170171‑M1/1‑48 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                <  1:219308‑M1/50‑4 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 <  1:4132391‑M1/49‑3 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 <  1:2024161‑M1/49‑3 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:5500322‑M1/1‑46 (MQ=255)
     cATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:1573091‑M1/1‑45 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                             <  1:4229602‑M1/50‑7 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              <  1:3688135‑M1/49‑6 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            <  1:4144104‑M1/50‑8 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                           <  1:1255410‑M1/49‑9 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                           <  1:1473173‑M1/49‑9 (MQ=255)
          tcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         >  1:3481580‑M1/1‑40 (MQ=255)
           cGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                        >  1:4324047‑M1/1‑39 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       <  1:1604640‑M1/50‑13 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                        <  1:2583365‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         <  1:3484710‑M1/48‑11 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         <  1:3562800‑M1/48‑11 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         <  1:862453‑M1/48‑11 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         <  1:2132791‑M1/48‑11 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                          <  1:75249‑M1/47‑10 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                          <  1:2479374‑M1/47‑10 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:4393878‑M1/50‑14 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:4843874‑M1/50‑14 (MQ=255)
               tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     <  1:3400158‑M1/49‑15 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccc                                   >  1:2054650‑M1/1‑34 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggcc                                    <  1:4458556‑M1/49‑16 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggcc                                    <  1:3369551‑M1/49‑16 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccca                                  <  1:3070769‑M1/50‑18 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccc                                   <  1:5174857‑M1/49‑17 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 <  1:1276466‑M1/50‑19 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 <  1:1386422‑M1/50‑19 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:5360421‑M1/1‑32 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                                <  1:5232726‑M1/50‑20 (MQ=255)
                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               >  1:3787089‑M1/1‑30 (MQ=255)
                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               >  1:5117182‑M1/1‑30 (MQ=255)
                     aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:2431827‑M1/49‑21 (MQ=255)
                     aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:3523582‑M1/49‑21 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                             <  1:4874575‑M1/50‑23 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              <  1:267558‑M1/49‑22 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              <  1:5175177‑M1/49‑22 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              <  1:5130667‑M1/49‑22 (MQ=255)
                       cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            >  1:4890039‑M1/1‑27 (MQ=255)
                       cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            >  1:843490‑M1/1‑27 (MQ=255)
                       cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            >  1:2802477‑M1/1‑27 (MQ=255)
                        aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                           >  1:5506133‑M1/1‑26 (MQ=255)
                        aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            <  1:2833261‑M1/49‑24 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                          <  1:3622549‑M1/50‑26 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                          <  1:164817‑M1/50‑26 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                          >  1:5351621‑M1/1‑25 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:1127502‑M1/1‑24 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:453691‑M1/1‑24 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:1598301‑M1/1‑24 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:2035694‑M1/1‑24 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:2286345‑M1/1‑24 (MQ=255)
                           cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:5188277‑M1/1‑23 (MQ=255)
                            aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:2099717‑M1/1‑22 (MQ=255)
                            aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgg                        <  1:1894734‑M1/49‑28 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccttgtacgcctggac                      <  1:2691549‑M1/50‑30 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      >  1:931826‑M1/1‑21 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     >  1:647719‑M1/1‑20 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     >  1:4363084‑M1/1‑20 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     >  1:845343‑M1/1‑20 (MQ=255)
                               ttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     >  1:4148383‑M1/1‑19 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     <  1:4427677‑M1/48‑31 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:2729854‑M1/1‑17 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:377605‑M1/1‑17 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 <  1:3801345‑M1/50‑35 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 <  1:4740203‑M1/50‑35 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  <  1:5387219‑M1/49‑34 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  <  1:3728163‑M1/49‑34 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  <  1:1230493‑M1/49‑34 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  >  1:3966844‑M1/1‑16 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  <  1:3551856‑M1/49‑34 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa               <  1:4195997‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 >  1:720626‑M1/1‑14 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 >  1:2642648‑M1/1‑14 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc             >  1:5398403‑M1/1‑12 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc             <  1:1086409‑M1/50‑39 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc             <  1:5429674‑M1/50‑39 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag              <  1:1077399‑M1/49‑38 (MQ=255)
                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagcg            <  1:546105‑M1/50‑40 (MQ=255)
                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt           <  1:3681323‑M1/50‑41 (MQ=255)
                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt           >  1:3508175‑M1/1‑10 (MQ=255)
                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt           <  1:3857659‑M1/50‑41 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt          <  1:260675‑M1/50‑42 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgttt         <  1:4690558‑M1/50‑43 (MQ=255)
                                            ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac       >  1:159384‑M1/1‑6 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac       <  1:4079234‑M1/49‑45 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac       <  1:1643162‑M1/49‑45 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac       <  1:4350281‑M1/49‑45 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc     <  1:5197259‑M1/49‑47 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc     <  1:1525619‑M1/49‑47 (MQ=255)
                                                 attccatataacggccctgtacgcctggtcggataaagcgtttacgcga   >  1:4865842‑M1/1‑1 (MQ=255)
                                                 attccatataacggccctgtacgcctggacggataaagcgtttacgcgac  <  1:2290302‑M1/50‑50 (MQ=255)
                                                 |                                                 
ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167434‑167532

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: