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

TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167435‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495
                                                                                                 
TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAT                                                 <  1:314773/50‑1
  ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               <  1:2395423/50‑1
  ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:1275352/1‑50
   TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCC                                              <  1:1929200/50‑1
   TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:1692686/1‑49
      TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                           >  1:693083/1‑50
      TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                           <  1:745008/50‑1
      TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                           <  1:710117/50‑1
       TCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                           >  1:1190388/1‑49
        CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:1472260/49‑1
        CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:560261/49‑1
        CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:734138/49‑1
           GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                       <  1:1155770/49‑1
            TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                     <  1:326160/50‑1
            TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                     <  1:1434060/50‑1
            TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                     <  1:2324366/50‑1
            TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                     >  1:2178806/1‑50
                 GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                >  1:1424584/1‑50
                 GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                >  1:79612/1‑50
                 GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                <  1:2529079/50‑1
                  TTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGT                               <  1:386472/50‑1
                  TTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGT                               <  1:1359119/50‑1
                  TTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                >  1:208619/1‑49
                    ATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                             >  1:822311/1‑50
                    ATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                             >  1:1688284/1‑50
                     TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACG                            <  1:1943769/50‑1
                     TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                             <  1:854417/49‑1
                     TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                             <  1:2082984/49‑1
                      CATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGC                           >  1:410145/1‑50
                       ATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCC                          >  1:360876/1‑50
                        TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCT                         >  1:1372035/1‑50
                         TCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                        >  1:1828355/1‑50
                          CACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                        >  1:283749/1‑49
                           ACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                      <  1:495500/50‑1
                             TTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACG                    >  1:2319112/1‑50
                               TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                  >  1:36191/1‑50
                                 CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                <  1:606347/50‑1
                                 CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                <  1:159698/50‑1
                                 CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGAT                 <  1:2218473/49‑1
                                 CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGAT                 <  1:2195786/49‑1
                                  ATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                >  1:598383/1‑49
                                    CAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG             >  1:2017384/1‑50
                                    CAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG             >  1:797677/1‑50
                                         ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT         >  1:422850/1‑49
                                           ATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC      >  1:1091676/1‑50
                                            TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACG     >  1:1933297/1‑50
                                            TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC      <  1:1167749/49‑1
                                             ACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGC    >  1:769645/1‑50
                                               CATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  1:1132325/1‑50
                                                                                                 
TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167435‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495

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.