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 64 [0] [0] 65 [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)64 (1.100) 48/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

TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167435‑167532
                                                |                                                 
tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  >  1:2326662‑M1/1‑49 (MQ=255)
 aaTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 >  1:4354446‑M1/1‑48 (MQ=255)
  aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:3303512‑M1/1‑47 (MQ=255)
  aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 <  1:518325‑M1/49‑3 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcc                                               <  1:2427534‑M1/50‑5 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:416501‑M1/1‑46 (MQ=255)
     aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                             <  1:754369‑M1/50‑7 (MQ=255)
     aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:2256534‑M1/1‑44 (MQ=255)
      ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            <  1:4172418‑M1/50‑8 (MQ=255)
      ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            <  1:602018‑M1/50‑8 (MQ=255)
       tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:1284605‑M1/1‑42 (MQ=255)
        ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                           >  1:908936‑M1/1‑41 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                        <  1:4258308‑M1/49‑12 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                          <  1:2581213‑M1/47‑10 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                          <  1:3878154‑M1/47‑10 (MQ=255)
            tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:2363588‑M1/50‑14 (MQ=255)
            tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:282205‑M1/50‑14 (MQ=255)
             ttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     >  1:3545715‑M1/1‑36 (MQ=255)
              tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     <  1:1689121‑M1/49‑15 (MQ=255)
                cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:1771974‑M1/1‑33 (MQ=255)
                cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:172888‑M1/1‑33 (MQ=255)
                cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:2517136‑M1/1‑33 (MQ=255)
                cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:4350371‑M1/1‑33 (MQ=255)
                cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  <  1:2587017‑M1/50‑18 (MQ=255)
                cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  <  1:2808297‑M1/50‑18 (MQ=255)
                 gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggctctg                                 >  1:780816‑M1/1‑32 (MQ=255)
                  ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                                <  1:2963731‑M1/50‑20 (MQ=255)
                    aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              >  1:2758392‑M1/1‑29 (MQ=255)
                    aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:1371678‑M1/49‑21 (MQ=255)
                    aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:2531060‑M1/49‑21 (MQ=255)
                      cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            >  1:983164‑M1/1‑27 (MQ=255)
                        ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                          >  1:2542078‑M1/1‑25 (MQ=255)
                        ttCACTTTACATCAGAGATATACCAttccatacaacggccctgtacgcct                          <  1:3230889‑M1/50‑26 (MQ=255)
                          cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         <  1:1020898‑M1/49‑27 (MQ=255)
                          cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:3704520‑M1/1‑23 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:2308068‑M1/1‑22 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:3995312‑M1/1‑22 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:3632921‑M1/1‑22 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:1617847‑M1/1‑22 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:3264027‑M1/1‑22 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:3013906‑M1/1‑22 (MQ=255)
                            cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      <  1:1540443‑M1/50‑30 (MQ=255)
                             tttACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     >  1:4214041‑M1/1‑20 (MQ=255)
                               tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:3219914‑M1/1‑18 (MQ=255)
                               tACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                     <  1:640433‑M1/48‑31 (MQ=255)
                                 cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  <  1:1458291‑M1/49‑34 (MQ=255)
                                  aTCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 >  1:3632277‑M1/1‑15 (MQ=255)
                                  aTCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 >  1:2615974‑M1/1‑15 (MQ=255)
                                     agagATATACCAttccatataacggccctgtacgcctggacggataaagc             >  1:2629194‑M1/1‑12 (MQ=255)
                                       agaTATACCAttccatataacggccctgtacgcctggacggataaagcgt           >  1:3425758‑M1/1‑10 (MQ=255)
                                        gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:2914531‑M1/1‑9 (MQ=255)
                                        gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt          <  1:1462112‑M1/50‑42 (MQ=255)
                                         atatACCAttccatataacggccctgtacgcctggacggataaagcgttt         <  1:2118341‑M1/50‑43 (MQ=255)
                                         atatACCAttccatataacggccctgtacgcctggacggataaagcgttt         <  1:1309924‑M1/50‑43 (MQ=255)
                                         atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:3308765‑M1/1‑8 (MQ=255)
                                         atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:2450591‑M1/1‑8 (MQ=255)
                                           ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac       >  1:249053‑M1/1‑6 (MQ=255)
                                            taCCAttccatataacggccctgtacgcctggacggataaagcgtttac       <  1:3332571‑M1/49‑45 (MQ=255)
                                              ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc     <  1:3525574‑M1/49‑47 (MQ=255)
                                              ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc     >  1:3404279‑M1/1‑3 (MQ=255)
                                               cAttccatataacggccctgtacgcctggacggataaagcgtttacgcga   <  1:4312295‑M1/50‑49 (MQ=255)
                                               cAttccatataacggccctgtacgcctggacggataaagcgtttacgcga   >  1:1716493‑M1/1‑2 (MQ=255)
                                                attccatataacggccctgtacgcctggacggataaagcgtttacgcgac  <  1:1239746‑M1/50‑50 (MQ=255)
                                                attccatataacggccctgtacgcctggacggataaagcgtttacgcga   >  1:1595134‑M1/1‑1 (MQ=255)
                                                |                                                 
TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167435‑167532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ 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.