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 27 [0] [0] 27 [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)26 (0.870) 25/98 0.1 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

AAAAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGT  >  NC_000913/3367594‑3367692
                                                 |                                                 
tgctcagacacattacattcactgtatttataacaacatcaatgaatact                                                   <  1:503848‑M2/1‑1 (MQ=255)
tgctcagacacattacattcactgtatttataacaacatcaatgaatact                                                   <  1:2044462‑M2/1‑1 (MQ=255)
        cacattacattcactgtatttataacaacatcaatgaatacTAACAATGc                                           >  1:139273‑M2/42‑50 (MQ=255)
             tacattcactgtatttataacaacatcaatgaatacTAACAATGCCGGa                                       >  1:1789693‑M2/37‑49 (MQ=255)
                attcactgtatttataacaacatcaatgaatacTAACAATGCCGGACAg                                    >  1:2196503‑M2/34‑49 (MQ=255)
                    actgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTa                               <  1:2002007‑M2/21‑1 (MQ=255)
                      tgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTAc                              >  1:786401‑M2/28‑49 (MQ=255)
                       gtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTa                            >  1:1481748‑M2/27‑50 (MQ=255)
                        tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACt                             <  1:2180795‑M2/23‑1 (MQ=255)
                         atttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           >  1:253035‑M2/25‑49 (MQ=255)
                             ataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAgg                       >  1:746085‑M2/21‑49 (MQ=255)
                              taacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGc                      <  1:1679007‑M2/30‑1 (MQ=255)
                                acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAtt                   >  1:1114906‑M2/18‑50 (MQ=255)
                                acaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCAt                    <  1:1647091‑M2/32‑1 (MQ=255)
                                  aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCg                 <  1:562597‑M2/35‑1 (MQ=255)
                                   acatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGt                >  1:1346098‑M2/15‑50 (MQ=255)
                                    catcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGtt               <  1:772090‑M2/37‑1 (MQ=255)
                                    catcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGtt               >  1:2137766‑M2/14‑50 (MQ=255)
                                     atcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGttt              <  1:2010278‑M2/38‑1 (MQ=255)
                                        aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta            <  1:1228408‑M2/40‑1 (MQ=255)
                                        aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta            >  1:21812‑M2/10‑49 (MQ=255)
                                           gaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAc         >  1:605388‑M2/7‑49 (MQ=255)
                                            aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACg        <  1:1835675‑M2/44‑1 (MQ=255)
                                            aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa       >  1:2189576‑M2/6‑50 (MQ=255)
                                              tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     <  1:215186‑M2/47‑1 (MQ=255)
                                                cTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCg    <  1:650764‑M2/48‑1 (MQ=255)
                                                 tAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGt  >  1:234892/1‑50 (MQ=255)
                                                 |                                                 
AAAAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGT  >  NC_000913/3367594‑3367692

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