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 133 [0] [0] 130 [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)130 (1.200) 72/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

Only 100 of 130 total aligned reads displayed.

CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173399‑173495
                                                 |                                               
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                 >  1:542346‑M2/50‑50 (MQ=255)
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                 >  1:4131072‑M2/50‑50 (MQ=255)
taatcattctcgtttacgttatcattcactttacatcagagatataccat                                                 <  1:7762542‑M2/1‑1 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                               <  1:6256128‑M2/3‑1 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                               >  1:6539371‑M2/48‑50 (MQ=255)
  atcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                               >  1:6628879‑M2/48‑50 (MQ=255)
   tcattctcgtttacgttatcattcactttacatcagagatataccaTTc                                               >  1:3882242‑M2/47‑49 (MQ=255)
    cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             >  1:5659477‑M2/46‑50 (MQ=255)
    cattctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             >  1:5883070‑M2/46‑50 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                            <  1:4627445‑M2/6‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                            <  1:4566428‑M2/6‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCat                                            <  1:7817210‑M2/6‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             <  1:793641‑M2/5‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             >  1:1755382‑M2/45‑49 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             <  1:4269606‑M2/5‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             <  1:3971121‑M2/5‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             <  1:4842962‑M2/5‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             <  1:1163038‑M2/5‑1 (MQ=255)
     attctcgtttacgttatcattcactttacatcagagatataccaTTCCa                                             <  1:7538853‑M2/5‑1 (MQ=255)
       tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                           >  1:1671358‑M2/43‑49 (MQ=255)
       tctcgtttacgttatcattcactttacatcagagatataccaTTCCata                                           >  1:6491837‑M2/43‑49 (MQ=255)
        ctcgtttacgttatcattcactttacatcagagatataccaTTCCatat                                          >  1:4640910‑M2/42‑49 (MQ=255)
         tcgtttacgttgtcattccctttacatcagagatataccaTTCCAtata                                         <  1:7032633‑M2/9‑1 (MQ=255)
          cgtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                       >  1:7267945‑M2/40‑50 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                         <  1:5207831‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCAtata                                         <  1:6475400‑M2/9‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                        <  1:7341902‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATaa                                        <  1:211972‑M2/10‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAAc                                       <  1:6925119‑M2/11‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                      <  1:3375160‑M2/12‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                      >  1:1619449‑M2/39‑50 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                      <  1:8218994‑M2/12‑1 (MQ=255)
           gtttacgttatcattcactttacatcagagatataccaTTCCATATAACg                                      <  1:4044347‑M2/12‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                     >  1:2901807‑M2/38‑50 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                     >  1:6097828‑M2/38‑50 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                     <  1:7048745‑M2/13‑1 (MQ=255)
            tttacgttatcattcactttacatcagagatataccaTTCCATATAACgg                                     >  1:6797880‑M2/38‑50 (MQ=255)
             ttacgttatcattcactttacatcagagatataccaTTCCATATAACGGc                                    >  1:5613322‑M2/37‑50 (MQ=255)
              tacgttatcattcactttacatcagagatataccaTTCCATATAACGGcc                                   >  1:6105055‑M2/36‑50 (MQ=255)
               acgttatcattcactttacatcagagatataccaTTCCATATAACGGccc                                  >  1:8309228‑M2/35‑50 (MQ=255)
                cgttatcattcactttacatcagagatataccaTTCCATATAACGGCCCt                                 >  1:5393449‑M2/34‑50 (MQ=255)
                 gttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                >  1:3967451‑M2/33‑50 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                >  1:4230489‑M2/32‑49 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                >  1:4991892‑M2/32‑49 (MQ=255)
                  ttatcattcactttacatcagagatataccaTTCCATATAACGGCCCTg                                >  1:2748650‑M2/32‑49 (MQ=255)
                   tatcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTa                              >  1:345730‑M2/31‑50 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                             <  1:2923175‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                             <  1:4836331‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                             <  1:2057880‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTAc                             <  1:3947050‑M2/21‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                            >  1:1674468‑M2/29‑50 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                            <  1:3808483‑M2/22‑1 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                            >  1:4345242‑M2/29‑50 (MQ=255)
                     tcattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACg                            <  1:3006055‑M2/22‑1 (MQ=255)
                      cattcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                           >  1:5494626‑M2/28‑50 (MQ=255)
                       attcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGc                           <  1:5301833‑M2/23‑1 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                         <  1:3733728‑M2/25‑1 (MQ=255)
                        ttcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCt                         >  1:4362876‑M2/26‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                        >  1:182341‑M2/25‑50 (MQ=255)
                         tcactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                        >  1:339372‑M2/25‑50 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                        <  1:1530563‑M2/26‑1 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                        <  1:2443087‑M2/26‑1 (MQ=255)
                          cactttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTg                        <  1:7744410‑M2/26‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTgg                       <  1:909815‑M2/27‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                      <  1:1039928‑M2/28‑1 (MQ=255)
                           actttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGa                      <  1:6123184‑M2/28‑1 (MQ=255)
                            ctttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGAc                     >  1:1291136‑M2/22‑50 (MQ=255)
                             tttacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                    >  1:7698796‑M2/21‑50 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                    <  1:7097595‑M2/30‑1 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                    <  1:4791497‑M2/30‑1 (MQ=255)
                               tacatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACg                    <  1:1886815‑M2/30‑1 (MQ=255)
                                acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                  >  1:1514646‑M2/18‑49 (MQ=255)
                                acatcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGa                  >  1:3587646‑M2/18‑49 (MQ=255)
                                 catcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                <  1:5488171‑M2/34‑1 (MQ=255)
                                  atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaa               <  1:2954402‑M2/35‑1 (MQ=255)
                                  atcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                >  1:3016836‑M2/16‑49 (MQ=255)
                                   tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa              <  1:5971560‑M2/36‑1 (MQ=255)
                                   tcagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATa                >  1:6508550‑M2/15‑48 (MQ=255)
                                    cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg             >  1:3561066‑M2/14‑50 (MQ=255)
                                    cagagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAg             >  1:5792309‑M2/14‑50 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc            <  1:3567085‑M2/38‑1 (MQ=255)
                                     agagatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGc            <  1:7715823‑M2/38‑1 (MQ=255)
                                       agatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGt          >  1:5265931‑M2/11‑50 (MQ=255)
                                        gatataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt         >  1:4387074‑M2/10‑50 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt        <  1:2789405‑M2/42‑1 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt         >  1:6784084‑M2/9‑49 (MQ=255)
                                         atataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt         >  1:6424647‑M2/9‑49 (MQ=255)
                                         atataccaTGCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGttt        <  1:8027499‑M2/42‑1 (MQ=255)
                                          tataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa       >  1:2502030‑M2/8‑50 (MQ=255)
                                          tataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTa       <  1:7867843‑M2/43‑1 (MQ=255)
                                           ataccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc      >  1:7495760‑M2/7‑50 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcg     >  1:5429825‑M2/6‑50 (MQ=255)
                                            taccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc      <  1:214296‑M2/44‑1 (MQ=255)
                                             accaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc    >  1:1052160‑M2/5‑50 (MQ=255)
                                             accaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc    >  1:3170195‑M2/5‑50 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgcg   <  1:3325065‑M2/47‑1 (MQ=255)
                                              ccaTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAcgc    <  1:6321077‑M2/46‑1 (MQ=255)
                                               caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa  <  1:7935855‑M2/48‑1 (MQ=255)
                                               caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa  >  1:191516‑M2/3‑50 (MQ=255)
                                               caTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa  >  1:1639662‑M2/3‑50 (MQ=255)
                                                 |                                               
CGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA  >  NC_000913/173399‑173495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: