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 31 [0] [0] 31 [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)30 (0.800) 27/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

AAATCGAACTTGATGGGAAAGCTGTGATAAGCGTAGTCTACTTGTGATAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGT  >  NC_000913/3367596‑3367692
                                               |                                                 
ctcagacacattacattcactgtatttataacaacatcaatgaatacTa                                                  <  1:2012194‑M2/2‑1 (MQ=255)
   agacacattacattcactgtatttataacaacatcaatgaatacTAACaa                                              >  1:645335‑M2/45‑50 (MQ=255)
     acacattacattcactgtatttataacaacatcaatgaatacTAACAATg                                            >  1:1947993‑M2/43‑50 (MQ=255)
        cattacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                         >  1:1944578‑M2/40‑50 (MQ=255)
        cattacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                         <  1:1775662‑M2/11‑1 (MQ=255)
         attacattcactgtatttataacaacatcaatgaatacTAACAATGCCg                                         >  1:843914‑M2/39‑49 (MQ=255)
          ttacattcactgtatttataacaacatcaatgaatacTAACAATGCCGGa                                       >  1:372409‑M2/38‑50 (MQ=255)
              attcactgtatttataacaacatcaatgaatacTAACAATGCCGGACAgg                                   <  1:1580386‑M2/17‑1 (MQ=255)
               ttcactgtatttataacaacatcaatgaatacTAACAATGCCGGACAg                                    >  1:1856207‑M2/33‑48 (MQ=255)
                 cactgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAg                                 <  1:1211660‑M2/19‑1 (MQ=255)
                   ctgtatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTAc                              >  1:2555435‑M2/29‑50 (MQ=255)
                      tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           <  1:160004‑M2/25‑1 (MQ=255)
                      tatttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           <  1:1084348‑M2/25‑1 (MQ=255)
                       atttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTAt                           >  1:2630049‑M2/25‑49 (MQ=255)
                         ttataacaacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAg                        >  1:2242563‑M2/23‑50 (MQ=255)
                               caacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTc                  <  1:270902‑M2/34‑1 (MQ=255)
                                aacatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCg                 <  1:1387714‑M2/35‑1 (MQ=255)
                                 acatcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGt                >  1:720124‑M2/15‑50 (MQ=255)
                                   atcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGttt              >  1:2230451‑M2/13‑50 (MQ=255)
                                    tcaatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGttt              >  1:1070831‑M2/12‑49 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTtat           <  1:439168‑M2/41‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTtat           <  1:124169‑M2/41‑1 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta            >  1:1516432‑M2/10‑49 (MQ=255)
                                      aatgaatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTta            <  1:1032792‑M2/40‑1 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa       >  1:2560957‑M2/6‑50 (MQ=255)
                                          aatacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGa       <  1:1276636‑M2/45‑1 (MQ=255)
                                           atacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGAt      <  1:579886‑M2/46‑1 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     >  1:2679970‑M2/4‑50 (MQ=255)
                                            tacTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATc     >  1:876858‑M2/4‑50 (MQ=255)
                                              cTAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCgg   >  1:1458051‑M2/2‑50 (MQ=255)
                                               tAACAATGCCGGACAGGAGTACTATCAGGCATTCGTTTTATACGATCGGt  >  1:721591/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: