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 34 [0] [0] 33 [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)33 (1.220) 25/96 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

ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167437‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC  >  NC_000913/173448‑173491
                                                                                           
ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATT                                            <  1:406868/49‑1
 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                           <  1:1186255/49‑1
 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                           >  1:331022/1‑49
   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                         >  1:404836/1‑49
   ATTCTCGTTTACGTTAACATTCACTTTACATCAGAGATATACCATACC                                          <  1:1354244/48‑1
      CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                       <  1:1205441/48‑1
      CTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                       <  1:1658203/48‑1
         GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                   >  1:826444/1‑49
          TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                  <  1:598017/49‑1
          TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                  >  1:445685/1‑49
          TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                  <  1:188798/49‑1
               GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCT                             <  1:1257853/49‑1
               GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCT                             <  1:647258/49‑1
                 TATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                            >  1:995966/1‑48
                   TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                         <  1:183268/49‑1
                   TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTA                          <  1:402410/48‑1
                      TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCC                      <  1:1587674/49‑1
                         ACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                    >  1:296447/1‑48
                          CTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                  >  1:1678083/1‑49
                              ACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA              >  1:1944758/1‑49
                              ACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA              >  1:1287392/1‑49
                               CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGAT             <  1:72292/49‑1
                                   AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG         >  1:537066/1‑49
                                   AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG         >  1:1384172/1‑49
                                   AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG         >  1:745666/1‑49
                                   AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG         >  1:1041533/1‑49
                                    GAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC        >  1:1743411/1‑49
                                     AGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCG       >  1:334918/1‑49
                                      GATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGT      >  1:1854702/1‑49
                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT     >  1:284888/1‑49
                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT     <  1:1515884/49‑1
                                          TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC  >  1:930974/1‑49
                                          TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTA   <  1:1469368/48‑1
                                                                                           
ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167437‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC  >  NC_000913/173448‑173491

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.