Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,362,162 Δ5,481 bp [yhcA][yhcF] [yhcA], yhcD, yhcE, insH‑10, yhcE, [yhcF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3362162 3367642 5481 40 [0] [0] 40 [yhcA]–[yhcF] [yhcA],yhcD,yhcE,insH‑10,yhcE,[yhcF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 33621610 (0.000)40 (1.470) 34/96 0.0 100% coding (50/675 nt) yhcA putative fimbrial chaperone YhcA
?NC_000913 3367643 = 0 (0.000)intergenic (+1/‑184) yhcF/yhcG DUF1120 domain‑containing protein YhcF/DUF1016 domain‑containing protein YhcG

AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGA  >  NC_000913/3367596‑3367687
                                               |                                            
ctcagacacattacattcactgtatttataacaacatcaatgaatact                                              >  1:109642‑M2/48‑48 (MQ=255)
ctcagacacattacattcactgtatttataacaacatcaatgaatact                                              <  1:1603‑M2/1‑1 (MQ=255)
ctcagacacattacattcactgtatttataacaacatcaatgaatacTa                                             <  1:1506430‑M2/2‑1 (MQ=255)
   agacacattacattcactgtatttataacaacatcaatgaatacTAACa                                          <  1:727667‑M2/5‑1 (MQ=255)
      cacattacattcactgtatttataacaacatcaatgaatacTAACAATg                                       >  1:1106343‑M2/42‑49 (MQ=255)
      cacattacattcactgtatttataacaacatcaatgaatacTAACAATg                                       <  1:1185162‑M2/8‑1 (MQ=255)
      cacattacattcactgtatttataacaacatcaatgaatacTAACAATg                                       <  1:1869011‑M2/8‑1 (MQ=255)
       acattacattcactgtatttataacaacatcaatgaatacTAACAATGc                                      <  1:1093386‑M2/9‑1 (MQ=255)
       acattacattcactgtatttataacaacatcaatgaatacTAACAATGc                                      >  1:1112936‑M2/41‑49 (MQ=255)
        cattacattcactgtatttataacaacatcaatgaatacTAACAATGcc                                     <  1:510645‑M2/10‑1 (MQ=255)
         attacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                    <  1:792271‑M2/11‑1 (MQ=255)
            acattcactgtatttataacaacatcaatgaatacTAACAATGCCGGa                                  >  1:578288‑M2/36‑48 (MQ=255)
                 cactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGa                             <  1:876585‑M2/18‑1 (MQ=255)
                 cactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGa                             <  1:850683‑M2/18‑1 (MQ=255)
                  actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGt                           <  1:554313‑M2/20‑1 (MQ=255)
                    tgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTAc                         >  1:1563278‑M2/28‑49 (MQ=255)
                      tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTa                       >  1:434984‑M2/26‑49 (MQ=255)
                       atttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                      >  1:500103‑M2/25‑49 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                      >  1:467566‑M2/24‑48 (MQ=255)
                         ttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                      >  1:531882‑M2/23‑47 (MQ=255)
                          tataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAg                   >  1:376273‑M2/22‑49 (MQ=255)
                              acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAt               <  1:419176‑M2/32‑1 (MQ=255)
                              acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAt               <  1:44243‑M2/32‑1 (MQ=255)
                                aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTc             >  1:1922132‑M2/16‑49 (MQ=255)
                                aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTc             >  1:679768‑M2/16‑49 (MQ=255)
                                 acatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCg            <  1:1262376‑M2/35‑1 (MQ=255)
                                 acatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCg            >  1:2076963‑M2/15‑49 (MQ=255)
                                  catcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGt           <  1:1063066‑M2/36‑1 (MQ=255)
                                    tcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGttt         <  1:19672‑M2/38‑1 (MQ=255)
                                     caatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGttt         >  1:1799096‑M2/11‑48 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGtttt        <  1:417326‑M2/39‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGtttt        <  1:210332‑M2/39‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta       <  1:1927846‑M2/40‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta       <  1:687884‑M2/40‑1 (MQ=255)
                                       atgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta       >  1:1561877‑M2/9‑48 (MQ=255)
                                        tgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTtata     >  1:1138877‑M2/8‑49 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc    <  1:51157‑M2/43‑1 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc    >  1:778247‑M2/6‑48 (MQ=255)
                                           atacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa  >  1:1372890‑M2/5‑49 (MQ=255)
                                           atacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa  <  1:1563680‑M2/45‑1 (MQ=255)
                                               |                                            
AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGA  >  NC_000913/3367596‑3367687

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: 
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.