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 34 [0] [0] 32 [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)32 (1.050) 28/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

CTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3362115‑3362161
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCG  >  NC_000913/3367643‑3367690
                                                                                               
CTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACTA                                                <  1:1963335/49‑1
 TCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACTAA                                               >  1:1908761/1‑49
       ACATTACATTCACTGTATTTATAACAACATCAATGAATACTAACAATGCC                                        >  1:1725184/1‑50
             CATTCACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAG                                  >  1:543281/1‑50
              ATTCACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAGG                                 <  1:2426936/50‑1
              ATTCACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAG                                  >  1:574252/1‑49
                 CACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAG                               >  1:1683108/1‑49
                 CACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAG                               >  1:280054/1‑49
                  ACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTA                             >  1:2616212/1‑50
                  ACTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGT                              <  1:245349/49‑1
                   CTGTATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTAC                            >  1:2376960/1‑50
                      TATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTAT                         <  1:258506/50‑1
                      TATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTAT                         >  1:1969456/1‑50
                      TATTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACT                           <  1:2131655/48‑1
                        TTTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCA                       <  1:698925/50‑1
                         TTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAG                      >  1:642645/1‑50
                         TTATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAG                      >  1:2405955/1‑50
                          TATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGG                     >  1:2580252/1‑50
                          TATAACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAG                      <  1:1116706/49‑1
                             AACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGGCAT                  >  1:894686/1‑50
                             AACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGGCAT                  >  1:2315704/1‑50
                             AACAACATCAATGAATACTAACAATGCCGGACAGGAGTAATATCAGGCAT                  >  1:2298846/1‑50
                              ACAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGGCATT                 <  1:1752853/50‑1
                               CAACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGGCATTC                <  1:621084/50‑1
                                 ACATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGT              >  1:683715/1‑50
                                   ATCAATGAATACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTT            <  1:1189887/50‑1
                                        TGAATACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAC       <  1:2518270/50‑1
                                        TGAATACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAC       >  1:57209/1‑50
                                         GAATACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATAC       >  1:567284/1‑49
                                          AATACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGA     >  1:690082/1‑50
                                            TACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATC   <  1:2210369/50‑1
                                             ACTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCG  >  1:2108581/1‑50
                                                                                               
CTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3362115‑3362161
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCG  >  NC_000913/3367643‑3367690

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