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] 64 [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.250) 45/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

AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167436‑167532
                                               |                                                 
aaTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 >  1:1886620‑M1/1‑48 (MQ=255)
 aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:3215566‑M1/1‑47 (MQ=255)
 aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                                >  1:286952‑M1/1‑47 (MQ=255)
 aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                 <  1:2882949‑M1/49‑3 (MQ=255)
  tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcc                                               <  1:945832‑M1/50‑5 (MQ=255)
   cATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:3252916‑M1/1‑45 (MQ=255)
    aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                             <  1:3022732‑M1/50‑7 (MQ=255)
    aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                             <  1:1741252‑M1/50‑7 (MQ=255)
    aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:1183108‑M1/1‑44 (MQ=255)
    aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                              >  1:3714962‑M1/1‑44 (MQ=255)
     ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            <  1:3808919‑M1/50‑8 (MQ=255)
      tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:1671240‑M1/1‑42 (MQ=255)
      tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:1829548‑M1/1‑42 (MQ=255)
      tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:1077604‑M1/1‑42 (MQ=255)
      tctcGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                            >  1:79758‑M1/1‑42 (MQ=255)
         cGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                        >  1:534355‑M1/1‑39 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       >  1:475410‑M1/1‑38 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       >  1:1128352‑M1/1‑38 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                       >  1:3054250‑M1/1‑38 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                        <  1:1876185‑M1/49‑12 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                        <  1:3344349‑M1/49‑12 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataa                                         <  1:2673693‑M1/48‑11 (MQ=255)
          gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                          <  1:1140108‑M1/47‑10 (MQ=255)
           tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                      <  1:1618163‑M1/50‑14 (MQ=255)
            ttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     >  1:1491876‑M1/1‑36 (MQ=255)
            ttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                     >  1:2394462‑M1/1‑36 (MQ=255)
                gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:774440‑M1/1‑32 (MQ=255)
                gTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                                 >  1:3948232‑M1/1‑32 (MQ=255)
                 ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                                <  1:3015332‑M1/50‑20 (MQ=255)
                 ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                                <  1:287236‑M1/50‑20 (MQ=255)
                  tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               >  1:466084‑M1/1‑30 (MQ=255)
                   aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:2601374‑M1/49‑21 (MQ=255)
                   aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                               <  1:855508‑M1/49‑21 (MQ=255)
                    tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              <  1:2820800‑M1/49‑22 (MQ=255)
                    tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                              <  1:1163000‑M1/49‑22 (MQ=255)
                     cATTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            >  1:1482842‑M1/1‑27 (MQ=255)
                      aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                           <  1:3129294‑M1/50‑25 (MQ=255)
                      aTTCACTTTACATCAGAGATATACCAttccatataacggccctgtacgc                            <  1:2187262‑M1/49‑24 (MQ=255)
                       ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                          <  1:3091542‑M1/50‑26 (MQ=255)
                          aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                       >  1:3833322‑M1/1‑22 (MQ=255)
                           cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctggac                      >  1:3277865‑M1/1‑21 (MQ=255)
                              tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:3688788‑M1/1‑18 (MQ=255)
                              tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   >  1:367625‑M1/1‑18 (MQ=255)
                               aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga                   <  1:209944‑M1/49‑33 (MQ=255)
                                cATCAGAGATATACCAttccatataacggccctgtacgcctggacggat                  <  1:486110‑M1/49‑34 (MQ=255)
                                 aTCAGAGATATACCAttccatataacggccctgtacgcctggacggataa                <  1:3631194‑M1/50‑36 (MQ=255)
                                  tCAGAGATATACCAttccatataacggccctgtacgcctggacggataaa               <  1:2132317‑M1/50‑37 (MQ=255)
                                  tCAGAGATATACCAttccatataacggccctgtacgcctggacggata                 >  1:2479453‑M1/1‑14 (MQ=255)
                                   cAGAGATATACCAttccatataacggccctgtacgcctggacggataaag              >  1:2075094‑M1/1‑13 (MQ=255)
                                    agagATATACCAttccatataacggccctgtacgcctggacggataaagc             <  1:209109‑M1/50‑39 (MQ=255)
                                    agagATATACCAttccatataacggccctgtacgcctggacggataaagc             >  1:3501566‑M1/1‑12 (MQ=255)
                                    agagATATACCAttccatataacggccctgtacgcctggacggataaagc             >  1:1489163‑M1/1‑12 (MQ=255)
                                    agagATATACCAttccatataacggccctgtacgcctggacggataaag              <  1:891859‑M1/49‑38 (MQ=255)
                                       gaTATACCAttccatataacggccctgtacgcctggacggataaagcgtt          >  1:2162361‑M1/1‑9 (MQ=255)
                                         tataCCAttccatataacggccctgtacgcctggacggataaagcgttta        <  1:2024857‑M1/50‑44 (MQ=255)
                                         tataCCAttccatataacggccctgtacgcctggacggataaagcgttta        <  1:1460072‑M1/50‑44 (MQ=255)
                                          ataCCAttccatataacggccctgtacgcctggacggataaagcgtttac       >  1:1397453‑M1/1‑6 (MQ=255)
                                           taCCAttccatataacggccctgtacgcctggacggataaagcgtttac       <  1:1197866‑M1/49‑45 (MQ=255)
                                             ccAttccatataacggccctgtacgcctggacggataaagcgtttacgcg    <  1:3606573‑M1/50‑48 (MQ=255)
                                             ccAttccatataacggccctgtacgcctggacggataaagcgtttacgcg    <  1:1137311‑M1/50‑48 (MQ=255)
                                              cAttccatataacggccctgtacgcctggacggataaagcgtttacgcga   <  1:467731‑M1/50‑49 (MQ=255)
                                               attccatataacggccctgtacgcctggacggataaagcgtttacgcgac  <  1:258886‑M1/50‑50 (MQ=255)
                                               attccatataacggccctgtacgcctggacggataaagcgtttacgcga   >  1:3777316‑M1/1‑1 (MQ=255)
                                               attccatataacggccctgtacgcctggacggataaagcgtttacgcga   >  1:2431989‑M1/1‑1 (MQ=255)
                                               |                                                 
AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCG  >  NC_000913/167436‑167532

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