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 49 [0] [0] 50 [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)49 (1.310) 37/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

TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167435‑167531
                                                |                                                
tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                 <  1:314773‑M1/50‑2 (MQ=255)
  aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               <  1:2395423‑M1/50‑4 (MQ=255)
  aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:1275352‑M1/1‑47 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcc                                              <  1:1929200‑M1/50‑5 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:1692686‑M1/1‑46 (MQ=255)
      ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           >  1:693083‑M1/1‑43 (MQ=255)
      ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           <  1:745008‑M1/50‑8 (MQ=255)
      ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           <  1:710117‑M1/50‑8 (MQ=255)
       tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                           >  1:1190388‑M1/1‑42 (MQ=255)
        ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:1472260‑M1/49‑9 (MQ=255)
        ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:560261‑M1/49‑9 (MQ=255)
        ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:734138‑M1/49‑9 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       <  1:1155770‑M1/49‑12 (MQ=255)
            tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     <  1:326160‑M1/50‑14 (MQ=255)
            tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     <  1:1434060‑M1/50‑14 (MQ=255)
            tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     <  1:2324366‑M1/50‑14 (MQ=255)
            tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                     >  1:2178806‑M1/1‑37 (MQ=255)
                 gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:1424584‑M1/1‑32 (MQ=255)
                 gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:79612‑M1/1‑32 (MQ=255)
                 gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                <  1:2529079‑M1/50‑19 (MQ=255)
                  ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                               <  1:386472‑M1/50‑20 (MQ=255)
                  ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                               <  1:1359119‑M1/50‑20 (MQ=255)
                  ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                >  1:208619‑M1/1‑31 (MQ=255)
                    aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             >  1:822311‑M1/1‑29 (MQ=255)
                    aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             >  1:1688284‑M1/1‑29 (MQ=255)
                     tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacg                            <  1:1943769‑M1/50‑23 (MQ=255)
                     tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             <  1:854417‑M1/49‑22 (MQ=255)
                     tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             <  1:2082984‑M1/49‑22 (MQ=255)
                      cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                           >  1:410145‑M1/1‑27 (MQ=255)
                       aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                          >  1:360876‑M1/1‑26 (MQ=255)
                        ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                         >  1:1372035‑M1/1‑25 (MQ=255)
                         tCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        >  1:1828355‑M1/1‑24 (MQ=255)
                          cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        >  1:283749‑M1/1‑23 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                      <  1:495500‑M1/50‑29 (MQ=255)
                             tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    >  1:2319112‑M1/1‑20 (MQ=255)
                               tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                  >  1:36191‑M1/1‑18 (MQ=255)
                                 cATCAGAGATATACCAttccatataacggccctgtacgcctggacggata                <  1:606347‑M1/50‑35 (MQ=255)
                                 cATCAGAGATATACCAttccatataacggccctgtacgcctggacggata                <  1:159698‑M1/50‑35 (MQ=255)
                                 cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                 <  1:2218473‑M1/49‑34 (MQ=255)
                                 cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                 <  1:2195786‑M1/49‑34 (MQ=255)
                                  aTCAGAGATATACCAttccatataacggccctgtacgcctggacggata                >  1:598383‑M1/1‑15 (MQ=255)
                                    cAGAGATATACCAttccatataacggccctgtacgcctggacggataaag             >  1:2017384‑M1/1‑13 (MQ=255)
                                    cAGAGATATACCAttccatataacggccctgtacgcctggacggataaag             >  1:797677‑M1/1‑13 (MQ=255)
                                         atatACCAttccatataacggccctgtacgcctggacggataaagcgtt         >  1:422850‑M1/1‑8 (MQ=255)
                                           ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac      >  1:1091676‑M1/1‑6 (MQ=255)
                                            taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg     >  1:1933297‑M1/1‑5 (MQ=255)
                                            taCCAttccatataacggccctgtacgcctggacggataaagcgtttac      <  1:1167749‑M1/49‑45 (MQ=255)
                                             aCCAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:769645‑M1/1‑4 (MQ=255)
                                               cAttccatataacggccctgtacgcctggacggataaagcgtttacgcga  >  1:1132325‑M1/1‑2 (MQ=255)
                                                |                                                
TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167435‑167531

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.