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 30 [0] [0] 31 [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)29 (1.060) 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
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173448‑173495
                                                                                               
ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATT                                                >  1:419337/1‑49
 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:65337/1‑49
 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               <  1:607357/49‑1
 TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                               >  1:1386509/1‑49
   ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                             >  1:1494114/1‑49
    TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCAT                                            <  1:1180724/49‑1
         GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                       >  1:1986716/1‑49
         GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:739880/46‑1
         GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATAT                                          <  1:1347945/46‑1
                  ATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGT                               <  1:1824309/48‑1
                   TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                             >  1:1372167/1‑49
                   TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                             <  1:559667/49‑1
                      TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCC                          <  1:592260/49‑1
                      TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCC                          >  1:482449/1‑49
                         ACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGG                       <  1:1078028/49‑1
                         ACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTG                        <  1:643648/48‑1
                          CTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                      <  1:1980633/49‑1
                           TTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGAC                     >  1:1559976/1‑49
                           TTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGAC                     >  1:322353/1‑49
                             TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACG                    <  1:311307/48‑1
                             TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACG                    <  1:952302/48‑1
                                ATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                <  1:528841/49‑1
                                   AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG             <  1:1371244/49‑1
                                    GAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC            <  1:1037018/49‑1
                                      GATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGT          >  1:825634/1‑49
                                       ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT         <  1:2038601/49‑1
                                          TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC      <  1:1894517/49‑1
                                             CATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGC    >  1:1221614/1‑48
                                              ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  1:953484/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.