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 28 [0] [0] 28 [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)28 (0.920) 21/96 0.2 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

ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167437‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495
                                                                                               
ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAT                                                 <  1:278372/48‑1
 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:1478946/1‑49
  CATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:1656571/1‑48
   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                             >  1:196521/1‑49
   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCC                                              <  1:110094/48‑1
   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCC                                              <  1:1998065/48‑1
      CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                           >  1:6024/1‑48
      CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                           >  1:1519965/1‑48
         GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                       <  1:981887/49‑1
          TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                      >  1:197760/1‑49
          TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                      <  1:701249/49‑1
          TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                      >  1:395559/1‑49
           TTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                     >  1:1345451/1‑49
               GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCT                                 >  1:1512570/1‑49
                      TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCC                          >  1:661138/1‑49
                        CACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                        >  1:1828001/1‑49
                        CACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                        >  1:795279/1‑49
                          CTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                      >  1:1963080/1‑49
                          CTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                      >  1:1562157/1‑49
                            TTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACG                    >  1:1755587/1‑49
                                    GAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC            >  1:1898198/1‑49
                                      GATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGT          >  1:2117195/1‑49
                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT         >  1:662248/1‑49
                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT         >  1:1991481/1‑49
                                            CCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGC    >  1:1308787/1‑49
                                             CATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGC    >  1:252599/1‑48
                                              ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  1:1348976/1‑49
                                              ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  1:1903749/1‑49
                                                                                               
ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167437‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495

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