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

TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167435‑167531
                                                |                                                
taatCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCa                                                  <  1:72441/49‑1 (MQ=255)
  aTCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAtt                                                >  1:419337‑M1/1‑47 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               <  1:607357‑M1/49‑4 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:1386509‑M1/1‑46 (MQ=255)
   tCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttc                                               >  1:65337‑M1/1‑46 (MQ=255)
     aTTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttcca                                             >  1:1494114‑M1/1‑44 (MQ=255)
      ttCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccat                                            <  1:1180724‑M1/49‑7 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatataac                                       >  1:1986716‑M1/1‑38 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:1347945‑M1/46‑9 (MQ=255)
           gTTTACGTTATCATTCACTTTACATCAGAGATATACCAttccatat                                          <  1:739880‑M1/46‑9 (MQ=255)
                    aTCATTCACTTTACATCAGAGATATACCAttccatataacggccctgt                               <  1:1824309‑M1/48‑20 (MQ=255)
                     tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             <  1:559667‑M1/49‑22 (MQ=255)
                     tCATTCACTTTACATCAGAGATATACCAttccatataacggccctgtac                             >  1:1372167‑M1/1‑28 (MQ=255)
                        ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                          <  1:592260‑M1/49‑25 (MQ=255)
                        ttCACTTTACATCAGAGATATACCAttccatataacggccctgtacgcc                          >  1:482449‑M1/1‑25 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgg                       <  1:1078028‑M1/49‑28 (MQ=255)
                           aCTTTACATCAGAGATATACCAttccatataacggccctgtacgcctg                        <  1:643648‑M1/48‑27 (MQ=255)
                            cTTTACATCAGAGATATACCAttccatataacggccctgtacgcctgga                      <  1:1980633‑M1/49‑29 (MQ=255)
                             tttACATCAGAGATATACCAttccatataacggccctgtacgcctggac                     >  1:322353‑M1/1‑20 (MQ=255)
                             tttACATCAGAGATATACCAttccatataacggccctgtacgcctggac                     >  1:1559976‑M1/1‑20 (MQ=255)
                               tACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    <  1:952302‑M1/48‑31 (MQ=255)
                               tACATCAGAGATATACCAttccatataacggccctgtacgcctggacg                    <  1:311307‑M1/48‑31 (MQ=255)
                                  aTCAGAGATATACCAttccatataacggccctgtacgcctggacggata                <  1:528841‑M1/49‑35 (MQ=255)
                                     agagATATACCAttccatataacggccctgtacgcctggacggataaag             <  1:1371244‑M1/49‑38 (MQ=255)
                                      gagaTATACCAttccatataacggccctgtacgcctggacggataaagc            <  1:1037018‑M1/49‑39 (MQ=255)
                                        gaTATACCAttccatataacggccctgtacgcctggacggataaagcgt          >  1:825634‑M1/1‑9 (MQ=255)
                                         atatACCAttccatataacggccctgtacgcctggacggataaagcgtt         <  1:2038601‑M1/49‑42 (MQ=255)
                                            taCCAttccatataacggccctgtacgcctggacggataaagcgtttac      <  1:1894517‑M1/49‑45 (MQ=255)
                                               cAttccatataacggccctgtacgcctggacggataaagcgtttacgc    >  1:1221614‑M1/1‑2 (MQ=255)
                                                attccatataacggccctgtacgcctggacggataaagcgtttacgcga  >  1:953484‑M1/1‑1 (MQ=255)
                                                |                                                
TAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATC  >  NC_000913/167435‑167531

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: