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

ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAA  >  NC_000913/167434‑167528
                                                 |                                             
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                               <  1:1164070/50‑1 (MQ=255)
ataatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                               <  1:1428199/50‑1 (MQ=255)
    tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                            >  1:602813‑M1/1‑46 (MQ=255)
     cATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                          >  1:1965786‑M1/1‑45 (MQ=255)
       ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccata                                        >  1:1920772‑M1/1‑43 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacg                                   >  1:2034387‑M1/1‑38 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                    <  1:1830658‑M1/49‑12 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                      <  1:676686‑M1/47‑10 (MQ=255)
            gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatata                                      <  1:1230760‑M1/47‑10 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                  >  1:1456448‑M1/1‑37 (MQ=255)
             tttACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacgg                                  >  1:1959950‑M1/1‑37 (MQ=255)
               tACGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggc                                 <  1:625345‑M1/49‑15 (MQ=255)
                aCGTTATCATTCACTTTACATCAGAGATATACCAttccatataacggccc                               >  1:2201642‑M1/1‑34 (MQ=255)
                   ttATCATTCACTTTACATCAGAGATATACCAttccatataacggccctg                             >  1:540361‑M1/1‑31 (MQ=255)
                    tATCATTCACTTTACATCAGAGATATACCAttccatataacggccctgta                           >  1:897479‑M1/1‑30 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                          <  1:1831087‑M1/49‑22 (MQ=255)
                      tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                          <  1:1362911‑M1/49‑22 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                      >  1:1705883‑M1/1‑25 (MQ=255)
                         ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcct                      <  1:1500503‑M1/50‑26 (MQ=255)
                                tACATCAGAGATATACCAttccatataacggccctgtacgcctggacgga               >  1:704119‑M1/1‑18 (MQ=255)
                                 aCATCAGAGATATACCAttccatataacggccctgtacgcctggacgga               >  1:288724‑M1/1‑17 (MQ=255)
                                   aTCAGAGATATACCAttccatataacggccctgtacgcctggacggataa            <  1:2094072‑M1/50‑36 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc         <  1:1812992‑M1/50‑39 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaagc         >  1:1149498‑M1/1‑12 (MQ=255)
                                      agagATATACCAttccatataacggccctgtacgcctggacggataaag          <  1:2121491‑M1/49‑38 (MQ=255)
                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagcg        >  1:1449230‑M1/1‑11 (MQ=255)
                                       gagaTATACCAttccatataacggccctgtacgcctggacggataaagcg        >  1:1130813‑M1/1‑11 (MQ=255)
                                          atatACCAttccatataacggccctgtacgcctggacggataaagcgtt      >  1:668100‑M1/1‑8 (MQ=255)
                                             taCCAttccatataacggccctgtacgcctggacggataaagcgtttacg  >  1:395912‑M1/1‑5 (MQ=255)
                                                 |                                             
ATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAA  >  NC_000913/167434‑167528

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: