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 29 [0] [0] 29 [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)27 (0.900) 23/98 0.1 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

TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167438‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACG  >  NC_000913/173448‑173492
                                                                                           
TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                            >  1:602813/1‑49
 CATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                          >  1:1965786/1‑50
   TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                        >  1:1920772/1‑50
        GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                   >  1:2034387/1‑50
        GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                    <  1:1830658/49‑1
        GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATA                                      <  1:676686/47‑1
        GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATA                                      <  1:1230760/47‑1
         TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                  >  1:1456448/1‑50
         TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                  >  1:1959950/1‑50
           TACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGC                                 <  1:625345/49‑1
            ACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCC                               >  1:2201642/1‑50
               TTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                             >  1:540361/1‑49
                TATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTA                           >  1:897479/1‑50
                  TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                          <  1:1831087/49‑1
                  TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                          <  1:1362911/49‑1
                     TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCT                      >  1:1705883/1‑50
                     TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCT                      <  1:1500503/50‑1
                            TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA               >  1:704119/1‑50
                             ACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA               >  1:288724/1‑49
                               ATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAA            <  1:2094072/50‑1
                                  AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC         <  1:1812992/50‑1
                                  AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC         >  1:1149498/1‑50
                                  AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG          <  1:2121491/49‑1
                                   GAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCG        >  1:1449230/1‑50
                                   GAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCG        >  1:1130813/1‑50
                                      ATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT      >  1:668100/1‑49
                                         TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACG  >  1:395912/1‑50
                                                                                           
TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167438‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACG  >  NC_000913/173448‑173492

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