Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 167,401 Δ5,915 bp fhuA[fhuB] fhuA, fhuC, fhuD, [fhuB]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 167401 173315 5915 36 [0] [0] 36 fhuA–[fhuB] fhuA,fhuC,fhuD,[fhuB]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1674000 (0.000)36 (0.590) 34/492 0.4 100% intergenic (+136/‑84) mrcB/fhuA peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
?NC_000913 173316 = 0 (0.000)coding (1855/1983 nt) fhuB iron(III) hydroxamate ABC transporter membrane subunit

GATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGG  >  NC_000913/167154‑167650
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
gATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttac                                                                                                                                                                                                                                                        >  2:297552‑M1/1‑247 (MQ=255)
                   agcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactg                                                                                                                                                                                                                                     >  1:310459‑M1/1‑228 (MQ=255)
                       aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAAGTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgt                                                                                                                                                                                                                                 <  1:52786‑M1/251‑28 (MQ=255)
                         aGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgg                                                                                                                                                                                                                               <  1:461289‑M1/251‑30 (MQ=255)
                            cGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcg                                                                                                                                                                                                                            >  1:299173‑M1/1‑219 (MQ=255)
                                 cAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGATGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTAACGCTTGCGAACCttactgctgg                                                                                                                                                                                                                                                  >  2:309510‑M1/1‑214 (MQ=255)
                                 cAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTGCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctgg                                                                                                                                                                                                                                                  <  1:309510‑M1/224‑11 (MQ=255)
                                           cTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttcc                                                                                                                                                                                                             >  1:165839‑M1/1‑204 (MQ=255)
                                                 aaGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattcc                                                                                                                                                                                                        >  1:198604‑M1/1‑198 (MQ=255)
                                                                gCGACGGTGTAGCCGGTGGGATCAAGGATATGTTTTGTAATAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggct                                                                                                                                                                                        <  2:54952‑M1/251‑69 (MQ=255)
                                                                     ggTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgt                                                                                                                                                                                   <  1:173557‑M1/251‑74 (MQ=255)
                                                                           gCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCGGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacct                                                                                                                                                                             <  2:430157‑M1/251‑80 (MQ=255)
                                                                                   gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggc                                                                                                                                                                     <  1:139351‑M1/251‑88 (MQ=255)
                                                                                       aaGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattcca                                                                                                                                                                                                       >  1:414744‑M1/1‑160 (MQ=255)
                                                                                       aaGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggggtgggcggatggtgctttttccattcca                                                                                                                                                                                                       <  2:414744‑M1/213‑54 (MQ=255)
                                                                                                 ttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTCGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatatttta                                                                                                                                                        <  1:262567‑M1/250‑101 (MQ=255)
                                                                                                          aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgtt                                                                                                                                              <  1:247675‑M1/251‑111 (MQ=255)
                                                                                                                aTCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggcggctgtcaacctttatcggcgcgccatattttatctatttgttgagaa                                                                                                                                         >  1:436715‑M1/1‑135 (MQ=255)
                                                                                                                        gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagcc                                                                                                                                <  1:297552‑M1/251‑125 (MQ=255)
                                                                                                                               aCCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatt                                                                                                                         <  2:116325‑M1/251‑132 (MQ=255)
                                                                                                                                ccGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcccccaattttatctatttgttgagaaagcagagcccttatttt                                                                                                                        >  2:100273‑M1/1‑119 (MQ=255)
                                                                                                                                        cGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatat                                                                                                                <  2:190816‑M1/251‑141 (MQ=255)
                                                                                                                                                 aTCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttagctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgggccatattttatctatttgttttgaaagcaaagccgttcatttttccatataacggccct                                                                                                       >  2:11430‑M1/1‑102 (MQ=255)
                                                                                                                                                               cATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacgg                                                                                         >  1:140880‑M1/1‑88 (MQ=255)
                                                                                                                                                                        tAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagaagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgt                                                                                >  1:19559‑M1/1‑79 (MQ=255)
                                                                                                                                                                                     ctttctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccag                                                                                                                                                                                                      >  1:9382‑M1/1‑66 (MQ=255)
                                                                                                                                                                                     ctttctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccag                                                                                                                                                                                                      <  2:9382‑M1/120‑55 (MQ=255)
                                                                                                                                                                                              cATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttttctctttttttagaaagcaaagccgttaatttttgcctataacggccctgtacgcctggccggataaaacgtttaccccaccaccgattcccc                                                           >  2:527‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                tCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGCTTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgccgtcaacccttatcggcgcgccatattttatctaattgttgagaaaggagagccgttaatttttccctataacggccctggaaggctgggcggattaagcgtttacgggacaacggctgcccggtg                                                        >  1:139606‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                 cTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgccaaactttatcggcgcgccctattttatcttttttttggggaagcagagacgttaatttttccatataacggccccggacgcctggccgggcaaagcgtttacgcgacaacgggtggcccttgc                                                       >  2:60704‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                    gttgttTATTAACCCTTCAGGAACGCTCAGATTGCGTCCCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgac                                                    <  2:166351‑M1/251‑201 (MQ=255)
                                                                                                                                                                                                    gttgttTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgac                                                    <  1:219899‑M1/251‑201 (MQ=255)
                                                                                                                                                                                                    gttgttTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgac                                                    <  1:151404‑M1/251‑201 (MQ=255)
                                                                                                                                                                                                                                        tACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccg                 >  1:298139‑M1/1‑15 (MQ=255)
                                                                                                                                                                                                                                                gCGAACCttactgctggtgttcgctgactggtgtgggcggacggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgtag         <  1:236091‑M1/250‑244 (MQ=255)
                                                                                                                                                                                                                                                      cttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaactttatcggagcgccatattttatctttttgttgggaaaggagagccgttaattttttcatataacggccctgtacgcctggaagggtaaagcgtttacgcgacaaaggatgcccgaagtcacgctggcgcgtcttatcgggcctacaaacggcccccaattgtagggcggct  >  1:123393‑M1/1‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
GATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGG  >  NC_000913/167154‑167650

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