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 64 [0] [0] 64 [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)64 (1.250) 45/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

AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167436‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173448‑173496
                                                                                                 
AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATT                                                 >  1:1886620/1‑50
 ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                                >  1:3215566/1‑50
 ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTC                                                >  1:286952/1‑50
 ATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATT                                                 <  1:2882949/49‑1
  TCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCC                                               <  1:945832/50‑1
   CATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                              >  1:3252916/1‑50
    ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCAT                                             <  1:3022732/50‑1
    ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCAT                                             <  1:1741252/50‑1
    ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                              >  1:1183108/1‑49
    ATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCA                                              >  1:3714962/1‑49
     TTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                            <  1:3808919/50‑1
      TCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                            >  1:1671240/1‑49
      TCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                            >  1:1829548/1‑49
      TCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                            >  1:1077604/1‑49
      TCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATA                                            >  1:79758/1‑49
         CGTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                        >  1:534355/1‑50
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                       >  1:475410/1‑50
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                       >  1:1128352/1‑50
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACG                                       >  1:3054250/1‑50
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                        <  1:1876185/49‑1
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAAC                                        <  1:3344349/49‑1
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAA                                         <  1:2673693/48‑1
          GTTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATA                                          <  1:1140108/47‑1
           TTTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGG                                      <  1:1618163/50‑1
            TTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGC                                     >  1:1491876/1‑50
            TTACGTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGC                                     >  1:2394462/1‑50
                GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                 >  1:774440/1‑50
                GTTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTG                                 >  1:3948232/1‑50
                 TTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGT                                <  1:3015332/50‑1
                 TTATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGT                                <  1:287236/50‑1
                  TATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTA                               >  1:466084/1‑50
                   ATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTA                               <  1:2601374/49‑1
                   ATCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTA                               <  1:855508/49‑1
                    TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                              <  1:2820800/49‑1
                    TCATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTAC                              <  1:1163000/49‑1
                     CATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGC                            >  1:1482842/1‑50
                      ATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCC                           <  1:3129294/50‑1
                      ATTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGC                            <  1:2187262/49‑1
                       TTCACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCT                          <  1:3091542/50‑1
                          ACTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGA                       >  1:3833322/1‑50
                           CTTTACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGAC                      >  1:3277865/1‑50
                              TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                   >  1:3688788/1‑50
                              TACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                   >  1:367625/1‑50
                               ACATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGA                   <  1:209944/49‑1
                                CATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGAT                  <  1:486110/49‑1
                                 ATCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAA                <  1:3631194/50‑1
                                  TCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAA               <  1:2132317/50‑1
                                  TCAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATA                 >  1:2479453/1‑48
                                   CAGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG              >  1:2075094/1‑50
                                    AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC             <  1:209109/50‑1
                                    AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC             >  1:3501566/1‑50
                                    AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGC             >  1:1489163/1‑50
                                    AGAGATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAG              <  1:891859/49‑1
                                       GATATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTT          >  1:2162361/1‑50
                                         TATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTA        <  1:2024857/50‑1
                                         TATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTA        <  1:1460072/50‑1
                                          ATACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC       >  1:1397453/1‑50
                                           TACCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAC       <  1:1197866/49‑1
                                             CCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCG    <  1:3606573/50‑1
                                             CCATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCG    <  1:1137311/50‑1
                                              CATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA   <  1:467731/50‑1
                                               ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  <  1:258886/50‑1
                                               ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA   >  1:3777316/1‑49
                                               ATTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA   >  1:2431989/1‑49
                                                                                                 
AATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167436‑167483
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGAC  >  NC_000913/173448‑173496

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