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 207 [0] [0] 204 [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)201 (1.350) 83/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 207 total aligned reads displayed.

ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167434‑167531
                                                 |                                                
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                  >  2:775888/1‑50 (MQ=255)
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                  >  1:1640030/1‑50 (MQ=255)
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                  >  1:87582/1‑50 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                 <  1:2846577‑M1/50‑2 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:4802997‑M1/1‑47 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:1442782‑M1/1‑47 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               <  1:5287009‑M1/50‑4 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                <  2:1030567‑M1/49‑3 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:4717482‑M1/1‑46 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  2:488507‑M1/1‑46 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                            <  2:4251645‑M1/50‑7 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                             <  2:3909306‑M1/49‑6 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           >  2:1357028‑M1/1‑43 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           >  2:4731388‑M1/1‑43 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          >  1:2731628‑M1/1‑41 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          >  2:3441477‑M1/1‑41 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          >  1:4066345‑M1/1‑41 (MQ=255)
           cGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       >  2:531070‑M1/1‑39 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      <  2:3179638‑M1/50‑13 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      >  2:20453‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      <  1:5630799‑M1/50‑13 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                      >  2:821544‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       <  2:4570809‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       <  1:2024319‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       <  1:2869023‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                         <  1:2985791‑M1/47‑10 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                         <  2:3645766‑M1/47‑10 (MQ=255)
            gTTTACGTTATCATTCACTATACATCAGAGATATACCAttccatata                                         <  1:5437114‑M1/47‑10 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     >  2:565815‑M1/1‑37 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     >  1:3435433‑M1/1‑37 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     >  2:117539‑M1/1‑37 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     <  1:2172545‑M1/50‑14 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     >  1:4400787‑M1/1‑37 (MQ=255)
               tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                    <  2:123690‑M1/49‑15 (MQ=255)
               tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                    <  1:843591‑M1/49‑15 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccc                                  >  1:1982683‑M1/1‑34 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  2:142902‑M1/1‑32 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:3718571‑M1/1‑32 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                               <  1:3628483‑M1/50‑20 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:2592363‑M1/1‑31 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:952034‑M1/1‑31 (MQ=255)
                     aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             >  2:2821274‑M1/1‑29 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                            <  1:2086260‑M1/50‑23 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             <  1:689862‑M1/49‑22 (MQ=255)
                       cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                           >  2:261169‑M1/1‑27 (MQ=255)
                        aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                          >  2:4913010‑M1/1‑26 (MQ=255)
                        aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                          <  1:5100835‑M1/50‑25 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                         >  2:1446823‑M1/1‑25 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAtaccatataacggccctgtacgcct                         <  1:4497439‑M1/50‑26 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        >  1:5476446‑M1/1‑24 (MQ=255)
                          tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        >  1:4589559‑M1/1‑24 (MQ=255)
                           cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        >  1:147537‑M1/1‑23 (MQ=255)
                            aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                      >  1:748843‑M1/1‑22 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                     >  2:2091715‑M1/1‑21 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                     <  1:1766255‑M1/50‑30 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                     <  1:4964921‑M1/50‑30 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    >  1:4600945‑M1/1‑20 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    >  2:1736553‑M1/1‑20 (MQ=255)
                              tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    >  1:3320975‑M1/1‑20 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:1879149‑M1/1‑18 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:2294396‑M1/1‑18 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  <  1:3918101‑M1/49‑33 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:1167283‑M1/1‑17 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:2053939‑M1/1‑17 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  2:5109974‑M1/1‑17 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggata                <  2:5151737‑M1/50‑35 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                 >  1:3670737‑M1/1‑16 (MQ=255)
                                   aTCAGAGATATACCAttccatataacggccctgtacgcctggacggata                >  1:591045‑M1/1‑15 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa              <  1:1220533‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa              <  1:2688949‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa              <  2:1935587‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                >  1:3973764‑M1/1‑14 (MQ=255)
                                      agagATATTCCAttccatataacggccctgtacgcctggacggataaag             <  2:5610180‑M1/49‑38 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc            >  2:3063512‑M1/1‑12 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc            >  1:3052383‑M1/1‑12 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc            >  2:1542875‑M1/1‑12 (MQ=255)
                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagcg           >  2:5427482‑M1/1‑11 (MQ=255)
                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagcg           >  1:3203310‑M1/1‑11 (MQ=255)
                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagcg           <  2:4083254‑M1/50‑40 (MQ=255)
                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt          >  1:3496215‑M1/1‑10 (MQ=255)
                                        agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt          <  2:4716352‑M1/50‑41 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt         >  2:5686878‑M1/1‑9 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt         <  2:4872514‑M1/50‑42 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt         >  1:5011940‑M1/1‑9 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgttt        <  2:3767803‑M1/50‑43 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgttt        <  1:1035695‑M1/50‑43 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta       >  2:1696949‑M1/1‑7 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta       <  2:3012529‑M1/50‑44 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta       >  1:1530455‑M1/1‑7 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta       <  2:2329607‑M1/50‑44 (MQ=255)
                                            ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac      >  1:1091284‑M1/1‑6 (MQ=255)
                                            ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac      >  1:3154927‑M1/1‑6 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg     <  1:2973794‑M1/50‑46 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac      <  2:5338615‑M1/49‑45 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac      <  1:1719012‑M1/49‑45 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac      <  2:3314546‑M1/49‑45 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac      <  2:2096846‑M1/49‑45 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:5406540‑M1/1‑3 (MQ=255)
                                                cAttccatataacggccctgtacgcctggacggataaagcgtttacgcga  <  2:690592‑M1/50‑49 (MQ=255)
                                                 attccatataacggccctgtacgcctggacggataaagcgtttacgcga  >  1:5178045‑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 < 31 ≤ 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.