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 36 [0] [0] 35 [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)34 (1.240) 29/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

AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGG  >  NC_000913/3367596‑3367691
                                               |                                                
ctcagacacattacattcactgtatttataacaacatcaatgaatact                                                  <  1:2116787‑M2/1‑1 (MQ=255)
ctcagacacattacattcactgtatttataacaacatcaatgaatacTa                                                 <  1:1729716‑M2/2‑1 (MQ=255)
 tcagacacattacattcactgtatttataacaacatcaatgaatacTaa                                                >  1:98383‑M2/47‑49 (MQ=255)
   agacacattacattcactgtatttttaacaacatcaatgaatacTAAc                                               >  1:83285‑M2/45‑48 (MQ=255)
    gacacattacattcactgtatttataacaacatcaatgaatacTAACaa                                             >  1:1279040‑M2/44‑49 (MQ=255)
     acacattacattcactgtatttataacaacatcaatgaatacTAACAAt                                            >  1:1352671‑M2/43‑49 (MQ=255)
     acacattacattcactgtatttataacaacatcaatgaatacTAACAAt                                            <  1:18623‑M2/7‑1 (MQ=255)
       acattacattcactgtatttataacaacatcaatgaatacTAACAATGc                                          >  1:176914‑M2/41‑49 (MQ=255)
        cattacattcactgtatttataacaacatcaatgaatacTAACAATGcc                                         <  1:535500‑M2/10‑1 (MQ=255)
        cattacattcactgtatttataacaacatcaatgaatacTAACAATGcc                                         >  1:1599791‑M2/40‑49 (MQ=255)
          ttacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                        >  1:1780363‑M2/38‑48 (MQ=255)
             cattcactgtatttataacaacatcaatgaatacTAACAATGCCGGAc                                     <  1:310782‑M2/14‑1 (MQ=255)
                 cactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAg                                >  1:2070135‑M2/31‑49 (MQ=255)
                  actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAg                                <  1:74196‑M2/19‑1 (MQ=255)
                  actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAg                                >  1:699049‑M2/30‑48 (MQ=255)
                  actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGt                               <  1:1030144‑M2/20‑1 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                          <  1:1379234‑M2/25‑1 (MQ=255)
                        tttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                          <  1:1400409‑M2/25‑1 (MQ=255)
                           ataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAg                       <  1:100551‑M2/28‑1 (MQ=255)
                           ataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAg                       <  1:795055‑M2/28‑1 (MQ=255)
                            taacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAgg                      <  1:240750‑M2/29‑1 (MQ=255)
                              acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCa                    <  1:1395704‑M2/31‑1 (MQ=255)
                              acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAt                   >  1:1503366‑M2/18‑49 (MQ=255)
                              acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAt                   <  1:1132501‑M2/32‑1 (MQ=255)
                                aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTc                 <  1:1414474‑M2/34‑1 (MQ=255)
                                aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTc                 <  1:1869511‑M2/34‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta           >  1:2116175‑M2/10‑49 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta           >  1:1090846‑M2/10‑49 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc        >  1:776288‑M2/6‑48 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc        <  1:803034‑M2/43‑1 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACg       <  1:1950150‑M2/44‑1 (MQ=255)
                                           atacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa      <  1:1143031‑M2/45‑1 (MQ=255)
                                              cTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc    <  1:2069957‑M2/47‑1 (MQ=255)
                                              cTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc    <  1:1965393‑M2/47‑1 (MQ=255)
                                               tAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCgg  >  1:210598/1‑49 (MQ=255)
                                               |                                                
AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGG  >  NC_000913/3367596‑3367691

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.