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

GGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173400‑173496
                                                |                                                
aatcattctcgtttacgttatcattcactttacatcagagatataccatt                                                 >  1:1886620‑M2/49‑50 (MQ=255)
 atcattctcgtttacgttatcattcactttacatcagagatataccatt                                                 <  1:2882949‑M2/2‑1 (MQ=255)
 atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                >  1:3215566‑M2/48‑50 (MQ=255)
 atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                                >  1:286952‑M2/48‑50 (MQ=255)
  tcattctcgtttacgttatcattcactttacatcagagatataccaTTcc                                               <  1:945832‑M2/4‑1 (MQ=255)
   cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                              >  1:3252916‑M2/46‑50 (MQ=255)
    attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                             <  1:3022732‑M2/6‑1 (MQ=255)
    attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                             <  1:1741252‑M2/6‑1 (MQ=255)
    attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                              >  1:1183108‑M2/45‑49 (MQ=255)
    attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                              >  1:3714962‑M2/45‑49 (MQ=255)
     ttctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            <  1:3808919‑M2/7‑1 (MQ=255)
      tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            >  1:1671240‑M2/43‑49 (MQ=255)
      tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            >  1:1829548‑M2/43‑49 (MQ=255)
      tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            >  1:1077604‑M2/43‑49 (MQ=255)
      tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                            >  1:79758‑M2/43‑49 (MQ=255)
         cgtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                        >  1:534355‑M2/40‑50 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                          <  1:1140108‑M2/9‑1 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                         <  1:2673693‑M2/10‑1 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                        <  1:1876185‑M2/11‑1 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                        <  1:3344349‑M2/11‑1 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                       >  1:475410‑M2/39‑50 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                       >  1:1128352‑M2/39‑50 (MQ=255)
          gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                       >  1:3054250‑M2/39‑50 (MQ=255)
           tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                      <  1:1618163‑M2/13‑1 (MQ=255)
            ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                     >  1:1491876‑M2/37‑50 (MQ=255)
            ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                     >  1:2394462‑M2/37‑50 (MQ=255)
                gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                 >  1:774440‑M2/33‑50 (MQ=255)
                gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                 >  1:3948232‑M2/33‑50 (MQ=255)
                 ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                <  1:3015332‑M2/19‑1 (MQ=255)
                 ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGt                                <  1:287236‑M2/19‑1 (MQ=255)
                  tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                               >  1:466084‑M2/31‑50 (MQ=255)
                   atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                               <  1:2601374‑M2/20‑1 (MQ=255)
                   atcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                               <  1:855508‑M2/20‑1 (MQ=255)
                    tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                              <  1:2820800‑M2/21‑1 (MQ=255)
                    tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                              <  1:1163000‑M2/21‑1 (MQ=255)
                     cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                            >  1:1482842‑M2/28‑50 (MQ=255)
                      attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGcc                           <  1:3129294‑M2/24‑1 (MQ=255)
                      attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                            <  1:2187262‑M2/23‑1 (MQ=255)
                       ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                          <  1:3091542‑M2/25‑1 (MQ=255)
                          actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                       >  1:3833322‑M2/23‑50 (MQ=255)
                           ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                      >  1:3277865‑M2/22‑50 (MQ=255)
                              tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                   >  1:3688788‑M2/19‑50 (MQ=255)
                              tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                   >  1:367625‑M2/19‑50 (MQ=255)
                               acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                   <  1:209944‑M2/32‑1 (MQ=255)
                                catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGAt                  <  1:486110‑M2/33‑1 (MQ=255)
                                 atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                <  1:3631194‑M2/35‑1 (MQ=255)
                                  tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa               <  1:2132317‑M2/36‑1 (MQ=255)
                                  tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                 >  1:2479453‑M2/15‑48 (MQ=255)
                                   cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg              >  1:2075094‑M2/14‑50 (MQ=255)
                                    agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg              <  1:891859‑M2/37‑1 (MQ=255)
                                    agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc             <  1:209109‑M2/38‑1 (MQ=255)
                                    agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc             >  1:3501566‑M2/13‑50 (MQ=255)
                                    agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc             >  1:1489163‑M2/13‑50 (MQ=255)
                                       gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt          >  1:2162361‑M2/10‑50 (MQ=255)
                                         tataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa        <  1:2024857‑M2/43‑1 (MQ=255)
                                         tataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa        <  1:1460072‑M2/43‑1 (MQ=255)
                                          ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc       >  1:1397453‑M2/7‑50 (MQ=255)
                                           taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc       <  1:1197866‑M2/44‑1 (MQ=255)
                                             ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgcg    <  1:3606573‑M2/47‑1 (MQ=255)
                                             ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgcg    <  1:1137311‑M2/47‑1 (MQ=255)
                                              caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa   <  1:467731‑M2/48‑1 (MQ=255)
                                               aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa   >  1:3777316‑M2/2‑49 (MQ=255)
                                               aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa   >  1:2431989‑M2/2‑49 (MQ=255)
                                               aTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAc  <  1:258886‑M2/49‑1 (MQ=255)
                                                |                                                
GGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173400‑173496

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.