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 64 [0] [0] 64 [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)62 (1.070) 41/98 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

ATGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACCTTAGCCACAGGGATGGTTCCTGAAACATCAGTATTGCTGGTTGATGAA  >  NC_000913/3362112‑3362210
                                                 |                                                 
aTGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATAc                                                   <  1:3752767/50‑1 (MQ=255)
 tGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATAc                                                   >  1:2816956/1‑49 (MQ=255)
   cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaa                                                <  1:2673649‑M1/50‑4 (MQ=255)
   cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaa                                                <  1:2912391‑M1/50‑4 (MQ=255)
   cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACta                                                 <  1:4177918‑M1/49‑3 (MQ=255)
   cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACta                                                 <  1:2969129‑M1/49‑3 (MQ=255)
   cTCAGACACATTACATTCACTGTATTTATAACAACAGCAATGAATACtaa                                                <  1:169669‑M1/50‑4 (MQ=255)
     cAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaac                                               >  1:3332467‑M1/1‑45 (MQ=255)
     cAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaac                                               >  1:4020902‑M1/1‑45 (MQ=255)
     cAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaac                                               >  1:3420306‑M1/1‑45 (MQ=255)
      aGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaa                                             >  1:473026‑M1/1‑44 (MQ=255)
      aGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaa                                             <  1:2041149‑M1/50‑7 (MQ=255)
      aGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaaca                                              <  1:3814664‑M1/49‑6 (MQ=255)
      aGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaac                                               <  1:2609363‑M1/48‑5 (MQ=255)
        acacATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaatg                                           <  1:3092134‑M1/50‑9 (MQ=255)
        acacATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaatg                                           <  1:3251509‑M1/50‑9 (MQ=255)
          acattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgcc                                         >  1:3513696‑M1/1‑40 (MQ=255)
          acattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgcc                                         >  1:410448‑M1/1‑40 (MQ=255)
          acattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgcc                                         >  1:4354906‑M1/1‑40 (MQ=255)
           cattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccg                                        <  1:4218926‑M1/50‑12 (MQ=255)
            attacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccg                                        >  1:399889‑M1/1‑38 (MQ=255)
             ttacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                      >  1:870207‑M1/1‑37 (MQ=255)
             ttacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                      >  1:1560597‑M1/1‑37 (MQ=255)
              tacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                      >  1:1247793‑M1/1‑36 (MQ=255)
                 attCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                   >  1:3554832‑M1/1‑33 (MQ=255)
                 attCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                   >  1:1189216‑M1/1‑33 (MQ=255)
                  ttCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacagga                                 <  1:2788590‑M1/50‑19 (MQ=255)
                    cACTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggag                                <  1:3521466‑M1/49‑20 (MQ=255)
                    cACTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggag                                >  1:1109935‑M1/1‑30 (MQ=255)
                     aCTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagt                               <  1:1952870‑M1/49‑21 (MQ=255)
                     aCTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagt                               <  1:3214893‑M1/49‑21 (MQ=255)
                      cTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                             <  1:4126740‑M1/50‑23 (MQ=255)
                       tGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                             >  1:936140‑M1/1‑27 (MQ=255)
                         tATTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                          >  1:1480062‑M1/1‑25 (MQ=255)
                         tATTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                          >  1:2647500‑M1/1‑25 (MQ=255)
                         tATTTATAACAACATCAATGAATACtaacaatgccggacaggagtacta                           <  1:1023356‑M1/49‑25 (MQ=255)
                          aTTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                          >  1:3657230‑M1/1‑24 (MQ=255)
                           tttATAACAACATCAATGAATACtaacaatgccggacaggagtactatca                        <  1:129936‑M1/50‑28 (MQ=255)
                           tttATAACAACATCAATGAATACtaacaatgccggacaggagtactatc                         <  1:1321195‑M1/49‑27 (MQ=255)
                           tttATAACAACATCAATGAATACtaacaatgccggacaggagtactatc                         <  1:780794‑M1/49‑27 (MQ=255)
                           tttATAACAACATCAATGAATACtaacaatgccggacaggagtactat                          >  1:611461‑M1/1‑23 (MQ=255)
                               taACAACATCAATGAATACtaacaatgccggacaggagtactatcaggca                    >  1:2249783‑M1/1‑19 (MQ=255)
                               taACAACATCAATGAATACtaacaatgccggacaggagtactatcaggca                    <  1:3855917‑M1/50‑32 (MQ=255)
                                 acaacaTCAATGAATACtaacaatgccggacaggagtactatcaggcatt                  <  1:926270‑M1/50‑34 (MQ=255)
                                  caacaTCAATGAATACtaacaatgccggacaggagtactatcaggcattc                 >  1:1310650‑M1/1‑16 (MQ=255)
                                   aacaTCAATGAATACtaacaatgccggacaggagtactatcaggcattcg                <  1:2726843‑M1/50‑36 (MQ=255)
                                     caTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgtt              <  1:821743‑M1/50‑38 (MQ=255)
                                     caTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgtt              <  1:1336119‑M1/50‑38 (MQ=255)
                                         aaTGAATACtaacaatgccggacaggagtactatcaggcattcgttttat          <  1:1769857‑M1/50‑42 (MQ=255)
                                          aTGAATACtaacaatgccggacaggagtactatcaggcattcgttttata         >  1:181403‑M1/1‑8 (MQ=255)
                                           tGAATACtaacaatgccggacaggagtactatcaggcattcgttttatac        >  1:3790406‑M1/1‑7 (MQ=255)
                                           tGAATACtaacaatgccggacaggagtactatcaggcattcgttttatac        >  1:2682376‑M1/1‑7 (MQ=255)
                                           tGAATACtaacaatgccggacaggagtactatcaggcattcgttttatac        >  1:2071450‑M1/1‑7 (MQ=255)
                                            gAATACtaacaatgccggacaggagtactatcaggcattcgttttatac        >  1:3018214‑M1/1‑6 (MQ=255)
                                             aaTACtaacaatgccggacaggagtactatcaggcattcgttttatacga      <  1:438805‑M1/50‑46 (MQ=255)
                                             aaTACtaacaatgccggacaggagtactatcaggcattcgttttatacga      <  1:1997982‑M1/50‑46 (MQ=255)
                                              aTACtaacaatgccggacaggagtactatcaggcattcgttttatacgat     <  1:1473504‑M1/50‑47 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc    <  1:3103873‑M1/50‑48 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc    >  1:1425854‑M1/1‑3 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc    <  1:1979195‑M1/50‑48 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc    >  1:877551‑M1/1‑3 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc    <  1:918808‑M1/50‑48 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc    <  1:2707415‑M1/50‑48 (MQ=255)
                                                 ctaacaatgccggacaggagtactatcaggcattcgttttatacgatcgg  >  1:4365525‑M1/1‑1 (MQ=255)
                                                 |                                                 
ATGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACCTTAGCCACAGGGATGGTTCCTGAAACATCAGTATTGCTGGTTGATGAA  >  NC_000913/3362112‑3362210

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