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 88 [0] [0] 88 [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 (1.480) 56/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

ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167434‑167532
                                                 |                                                 
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                   >  1:4077896/1‑50 (MQ=255)
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                   >  1:1196243/1‑50 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  <  1:2632117‑M1/50‑2 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  >  1:3133389‑M1/1‑49 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  <  1:3806219‑M1/50‑2 (MQ=255)
 tAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAt                                                  <  1:1699518‑M1/50‑2 (MQ=255)
  aaTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 >  1:4421285‑M1/1‑48 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                <  1:3657047‑M1/50‑4 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:3286654‑M1/1‑47 (MQ=255)
   aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 <  1:4065300‑M1/49‑3 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcc                                               <  1:1247218‑M1/50‑5 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:2946353‑M1/1‑46 (MQ=255)
     cATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:2797197‑M1/1‑45 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                             <  1:2833939‑M1/50‑7 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                             <  1:2112503‑M1/50‑7 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:1760616‑M1/1‑44 (MQ=255)
      aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:1293871‑M1/1‑44 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            <  1:487151‑M1/50‑8 (MQ=255)
        tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:1977743‑M1/1‑42 (MQ=255)
        tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:320241‑M1/1‑42 (MQ=255)
         ctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                           <  1:290576‑M1/49‑9 (MQ=255)
          tcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         >  1:2750086‑M1/1‑40 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       >  1:499282‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       >  1:764059‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       <  1:2919730‑M1/50‑13 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                          <  1:2756197‑M1/47‑10 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:4144374‑M1/50‑14 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:3426013‑M1/50‑14 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      >  1:1124148‑M1/1‑37 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:1246043‑M1/50‑14 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      >  1:1974038‑M1/1‑37 (MQ=255)
              ttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     >  1:2420353‑M1/1‑36 (MQ=255)
              ttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     >  1:1372120‑M1/1‑36 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggcc                                    <  1:3115756‑M1/49‑16 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  <  1:1075346‑M1/50‑18 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:3327564‑M1/1‑33 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:2665141‑M1/1‑33 (MQ=255)
                 cGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccct                                  >  1:1907514‑M1/1‑33 (MQ=255)
                  gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:268435‑M1/1‑32 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:2447899‑M1/1‑31 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:3263420‑M1/1‑31 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:916993‑M1/1‑31 (MQ=255)
                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               >  1:898485‑M1/1‑30 (MQ=255)
                     aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:938979‑M1/49‑21 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              <  1:649159‑M1/49‑22 (MQ=255)
                       cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            >  1:2865991‑M1/1‑27 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                          <  1:3285811‑M1/50‑26 (MQ=255)
                           cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         >  1:3776723‑M1/1‑23 (MQ=255)
                           cACTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                         <  1:428098‑M1/49‑27 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      >  1:1499854‑M1/1‑21 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      >  1:526372‑M1/1‑21 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      <  1:3278272‑M1/50‑30 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      >  1:2433644‑M1/1‑21 (MQ=255)
                             cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      >  1:2622235‑M1/1‑21 (MQ=255)
                               ttACATCAGAGATATTCCAttccatataacggccctgtacgcctggacg                     >  1:2723134‑M1/1‑19 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:1081765‑M1/1‑18 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:3752122‑M1/1‑17 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  >  1:38509‑M1/1‑16 (MQ=255)
                                  cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  >  1:2693769‑M1/1‑16 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa               <  1:3408867‑M1/50‑37 (MQ=255)
                                    tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 >  1:748380‑M1/1‑14 (MQ=255)
                                     cAGAGATATACCAttccatataacggccctgtacgcctggacggataaag              >  1:1209689‑M1/1‑13 (MQ=255)
                                     cAGAGATATACCAttccatataacggccctgtacgcctggacggataaag              >  1:150999‑M1/1‑13 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc             <  1:569382‑M1/50‑39 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag              <  1:3328043‑M1/49‑38 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:174626‑M1/1‑9 (MQ=255)
                                         gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt          <  1:3893‑M1/50‑42 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgttt         <  1:3293869‑M1/50‑43 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgttt         <  1:148149‑M1/50‑43 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:2507065‑M1/1‑8 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:2065908‑M1/1‑8 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:458287‑M1/1‑8 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:268136‑M1/1‑8 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:2059928‑M1/1‑8 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:1694753‑M1/1‑8 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:3493280‑M1/1‑8 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta        <  1:2188568‑M1/50‑44 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta        <  1:4424489‑M1/50‑44 (MQ=255)
                                           tataCCAttccatataacggccctgtacgcctggacggataaagcgttta        <  1:833227‑M1/50‑44 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg      <  1:2506346‑M1/50‑46 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg      <  1:1969598‑M1/50‑46 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg      >  1:524468‑M1/1‑5 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttac       <  1:3932032‑M1/49‑45 (MQ=255)
                                              aCCAttccatataacggccctgtacgcctggacggataaagcgtttacgc     >  1:569058‑M1/1‑4 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgcg    <  1:293992‑M1/50‑48 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgcg    <  1:3349746‑M1/50‑48 (MQ=255)
                                               ccAttccatataacggccctgtacgcctggacggataaagcgtttacgc     >  1:2060280‑M1/1‑3 (MQ=255)
                                                 attccatataacggccctgtacgcctggacggataaagcgtttacgcgac  <  1:3571790‑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 < 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.