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

ATGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACCTTAGCCACAGGGATGGTTCCTGAAACATCAGTATTGCTGGTTGATGA  >  NC_000913/3362112‑3362209
                                                 |                                                
aTGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATAc                                                  <  1:2176887/50‑1 (MQ=255)
 tGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACt                                                 <  1:503848‑M1/50‑2 (MQ=255)
 tGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACt                                                 <  1:2044462‑M1/50‑2 (MQ=255)
         cacaTTACATTCACTGTATTTATAACAACATCAATGAATACtaacaatgc                                         >  1:139273‑M1/1‑41 (MQ=255)
              tacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                     >  1:1789693‑M1/1‑36 (MQ=255)
                 attCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                  >  1:2196503‑M1/1‑33 (MQ=255)
                     aCTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagta                             <  1:2002007‑M1/50‑22 (MQ=255)
                       tGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                            >  1:786401‑M1/1‑27 (MQ=255)
                        gTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtacta                          >  1:1481748‑M1/1‑26 (MQ=255)
                         tATTTATAACAACATCAATGAATACtaacaatgccggacaggagtact                           <  1:2180795‑M1/48‑24 (MQ=255)
                          aTTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                         >  1:253035‑M1/1‑24 (MQ=255)
                              ataACAACATCAATGAATACtaacaatgccggacaggagtactatcagg                     >  1:746085‑M1/1‑20 (MQ=255)
                               taACAACATCAATGAATACtaacaatgccggacaggagtactatcaggc                    <  1:1679007‑M1/49‑31 (MQ=255)
                                 acaacaTCAATGAATACtaacaatgccggacaggagtactatcaggcatt                 >  1:1114906‑M1/1‑17 (MQ=255)
                                 acaacaTCAATGAATACtaacaatgccggacaggagtactatcaggcat                  <  1:1647091‑M1/49‑33 (MQ=255)
                                   aacaTCAATGAATACtaacaatgccggacaggagtactatcaggcattcg               <  1:562597‑M1/50‑36 (MQ=255)
                                    acaTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgt              >  1:1346098‑M1/1‑14 (MQ=255)
                                     caTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgtt             <  1:772090‑M1/50‑38 (MQ=255)
                                     caTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgtt             >  1:2137766‑M1/1‑13 (MQ=255)
                                      aTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgttt            <  1:2010278‑M1/50‑39 (MQ=255)
                                         aaTGAATACtaacaatgccggacaggagtactatcaggcattcgtttta          <  1:1228408‑M1/49‑41 (MQ=255)
                                         aaTGAATACtaacaatgccggacaggagtactatcaggcattcgtttta          >  1:21812‑M1/1‑9 (MQ=255)
                                            gAATACtaacaatgccggacaggagtactatcaggcattcgttttatac       >  1:605388‑M1/1‑6 (MQ=255)
                                             aaTACtaacaatgccggacaggagtactatcaggcattcgttttatacga     >  1:2189576‑M1/1‑5 (MQ=255)
                                             aaTACtaacaatgccggacaggagtactatcaggcattcgttttatacg      <  1:1835675‑M1/49‑45 (MQ=255)
                                               tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc   <  1:215186‑M1/50‑48 (MQ=255)
                                                 ctaacaatgccggacaggagtactatcaggcattcgttttatacgatcg  <  1:650764‑M1/49‑49 (MQ=255)
                                                 |                                                
ATGCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACCTTAGCCACAGGGATGGTTCCTGAAACATCAGTATTGCTGGTTGATGA  >  NC_000913/3362112‑3362209

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: