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

AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGT  >  NC_000913/3367596‑3367692
                                               |                                                 
ctcagacacattacattcactgtatttataacaacatcaatgaatacTaa                                                 <  1:2912391‑M2/3‑1 (MQ=255)
ctcagacacattacattcactgtatttataacaacatcaatgaatacTaa                                                 <  1:2673649‑M2/3‑1 (MQ=255)
ctcagacacattacattcactgtatttataacaacatcaatgaatacTa                                                  <  1:2969129‑M2/2‑1 (MQ=255)
ctcagacacattacattcactgtatttataacaacatcaatgaatacTa                                                  <  1:4177918‑M2/2‑1 (MQ=255)
ctcagacacattacattcactgtatttataacaacagcaatgaatacTaa                                                 <  1:169669‑M2/3‑1 (MQ=255)
  cagacacattacattcactgtatttataacaacatcaatgaatacTAAc                                                >  1:3332467‑M2/46‑49 (MQ=255)
  cagacacattacattcactgtatttataacaacatcaatgaatacTAAc                                                >  1:3420306‑M2/46‑49 (MQ=255)
  cagacacattacattcactgtatttataacaacatcaatgaatacTAAc                                                >  1:4020902‑M2/46‑49 (MQ=255)
   agacacattacattcactgtatttataacaacatcaatgaatacTAAc                                                <  1:2609363‑M2/4‑1 (MQ=255)
   agacacattacattcactgtatttataacaacatcaatgaatacTAACaa                                              >  1:473026‑M2/45‑50 (MQ=255)
   agacacattacattcactgtatttataacaacatcaatgaatacTAACaa                                              <  1:2041149‑M2/6‑1 (MQ=255)
   agacacattacattcactgtatttataacaacatcaatgaatacTAACa                                               <  1:3814664‑M2/5‑1 (MQ=255)
     acacattacattcactgtatttataacaacatcaatgaatacTAACAATg                                            <  1:3092134‑M2/8‑1 (MQ=255)
     acacattacattcactgtatttataacaacatcaatgaatacTAACAATg                                            <  1:3251509‑M2/8‑1 (MQ=255)
       acattacattcactgtatttataacaacatcaatgaatacTAACAATGcc                                          >  1:410448‑M2/41‑50 (MQ=255)
       acattacattcactgtatttataacaacatcaatgaatacTAACAATGcc                                          >  1:4354906‑M2/41‑50 (MQ=255)
       acattacattcactgtatttataacaacatcaatgaatacTAACAATGcc                                          >  1:3513696‑M2/41‑50 (MQ=255)
        cattacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                         <  1:4218926‑M2/11‑1 (MQ=255)
         attacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                         >  1:399889‑M2/39‑49 (MQ=255)
          ttacattcactgtatttataacaacatcaatgaatacTAACAATGCCGGa                                       >  1:870207‑M2/38‑50 (MQ=255)
          ttacattcactgtatttataacaacatcaatgaatacTAACAATGCCGGa                                       >  1:1560597‑M2/38‑50 (MQ=255)
           tacattcactgtatttataacaacatcaatgaatacTAACAATGCCGGa                                       >  1:1247793‑M2/37‑49 (MQ=255)
              attcactgtatttataacaacatcaatgaatacTAACAATGCCGGACAg                                    >  1:3554832‑M2/34‑49 (MQ=255)
              attcactgtatttataacaacatcaatgaatacTAACAATGCCGGACAg                                    >  1:1189216‑M2/34‑49 (MQ=255)
               ttcactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGa                                  <  1:2788590‑M2/18‑1 (MQ=255)
                 cactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAg                                 <  1:3521466‑M2/19‑1 (MQ=255)
                 cactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAg                                 >  1:1109935‑M2/31‑49 (MQ=255)
                  actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGt                                <  1:3214893‑M2/20‑1 (MQ=255)
                  actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGt                                <  1:1952870‑M2/20‑1 (MQ=255)
                   ctgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTAc                              <  1:4126740‑M2/22‑1 (MQ=255)
                    tgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTAc                              >  1:936140‑M2/28‑49 (MQ=255)
                      tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTa                            <  1:1023356‑M2/24‑1 (MQ=255)
                      tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           >  1:2647500‑M2/26‑50 (MQ=255)
                      tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           >  1:1480062‑M2/26‑50 (MQ=255)
                       atttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           >  1:3657230‑M2/25‑49 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           >  1:611461‑M2/24‑48 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATc                          <  1:780794‑M2/26‑1 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATc                          <  1:1321195‑M2/26‑1 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCa                         <  1:129936‑M2/27‑1 (MQ=255)
                            taacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCa                     <  1:3855917‑M2/31‑1 (MQ=255)
                            taacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCa                     >  1:2249783‑M2/20‑50 (MQ=255)
                              acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAtt                   <  1:926270‑M2/33‑1 (MQ=255)
                               caacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTc                  >  1:1310650‑M2/17‑50 (MQ=255)
                                aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCg                 <  1:2726843‑M2/35‑1 (MQ=255)
                                  catcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGtt               <  1:821743‑M2/37‑1 (MQ=255)
                                  catcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGtt               <  1:1336119‑M2/37‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTtat           <  1:1769857‑M2/41‑1 (MQ=255)
                                       atgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTtata          >  1:181403‑M2/9‑50 (MQ=255)
                                        tgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc         >  1:2071450‑M2/8‑50 (MQ=255)
                                        tgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc         >  1:2682376‑M2/8‑50 (MQ=255)
                                        tgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc         >  1:3790406‑M2/8‑50 (MQ=255)
                                         gaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc         >  1:3018214‑M2/7‑49 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa       <  1:438805‑M2/45‑1 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa       <  1:1997982‑M2/45‑1 (MQ=255)
                                           atacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGAt      <  1:1473504‑M2/46‑1 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     <  1:2707415‑M2/47‑1 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     <  1:918808‑M2/47‑1 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     >  1:877551‑M2/4‑50 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     <  1:3103873‑M2/47‑1 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     >  1:1425854‑M2/4‑50 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     <  1:1979195‑M2/47‑1 (MQ=255)
                                              cTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCgg   >  1:4365525‑M2/2‑50 (MQ=255)
                                               tAACCATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGt  >  1:3738484/1‑50 (MQ=255)
                                               tAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGt  >  1:3038195/1‑50 (MQ=255)
                                               |                                                 
AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGT  >  NC_000913/3367596‑3367692

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: