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 51 [0] [0] 51 [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)51 (0.990) 43/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

GCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACCTTAGCCACAGGGATGGTTCCTGAAACATCAGTATTGCTGGTTGATG  >  NC_000913/3362114‑3362208
                                               |                                               
gCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACt                                                >  1:1678426‑M1/1‑48 (MQ=255)
 cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaa                                              <  1:3264555‑M1/50‑4 (MQ=255)
 cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaa                                              >  1:796400‑M1/1‑47 (MQ=255)
 cTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACta                                               <  1:1888705‑M1/49‑3 (MQ=255)
    aGACACATTACATTCACTGTATTTATAACAACATCAATAAATACtaac                                             <  1:692170‑M1/48‑5 (MQ=255)
     gACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaat                                          <  1:972128‑M1/50‑8 (MQ=255)
     gACACATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaat                                          <  1:2705522‑M1/50‑8 (MQ=255)
      acacATTACATTCACTGTATTTATAACAACATCAATGAATACtaacaatg                                         <  1:1270826‑M1/50‑9 (MQ=255)
        acattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgcc                                       >  1:809527‑M1/1‑40 (MQ=255)
        acattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgcc                                       <  1:3176262‑M1/50‑11 (MQ=255)
        acattacattCACTGTATTTATAACAACATCAATGAATACtaacaatgcc                                       >  1:1813505‑M1/1‑40 (MQ=255)
           ttacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                    >  1:3021199‑M1/1‑37 (MQ=255)
            tacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                    >  1:2968490‑M1/1‑36 (MQ=255)
            tacattCACTGTATTTATAACAACATCAATGAATACtaacaatgccgga                                    >  1:2569470‑M1/1‑36 (MQ=255)
             acattCACTGTATTTATAACAACATCAATGAATACtaacaatgccggac                                   <  1:2287291‑M1/49‑15 (MQ=255)
              cattCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                 <  1:3825724‑M1/50‑17 (MQ=255)
              cattCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                 >  1:2242055‑M1/1‑34 (MQ=255)
              cattCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                 >  1:3502722‑M1/1‑34 (MQ=255)
              cattCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                 <  1:2674866‑M1/50‑17 (MQ=255)
               attCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                 >  1:1607232‑M1/1‑33 (MQ=255)
                ttCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacagga                               <  1:1149533‑M1/50‑19 (MQ=255)
                ttCACTGTATTTATAACAACATCAATGAATACtaacaatgccggacag                                 >  1:3759826‑M1/1‑32 (MQ=255)
                  cACTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggag                              >  1:2352192‑M1/1‑30 (MQ=255)
                   aCTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagta                            <  1:1418154‑M1/50‑22 (MQ=255)
                   aCTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagta                            >  1:1842370‑M1/1‑29 (MQ=255)
                    cTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                           >  1:746636‑M1/1‑28 (MQ=255)
                    cTGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                           <  1:987122‑M1/50‑23 (MQ=255)
                     tGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                           >  1:1563630‑M1/1‑27 (MQ=255)
                     tGTATTTATAACAACATCAATGAATACtaacaatgccggacaggagtac                           >  1:1617453‑M1/1‑27 (MQ=255)
                       tATTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                        <  1:1038439‑M1/50‑26 (MQ=255)
                       tATTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                        >  1:2445489‑M1/1‑25 (MQ=255)
                        aTTTATAACAACATCAATGAATACtaacaatgccggacaggagtactat                        >  1:3629045‑M1/1‑24 (MQ=255)
                          ttATAACAACATCAATGAATACtaacaatgccggacaggagtactatcag                     >  1:1053332‑M1/1‑22 (MQ=255)
                          ttATAACAACATCAATGAATACtaacaatgccggacaggagtactatcag                     >  1:530398‑M1/1‑22 (MQ=255)
                           tataACAACATCAATGAATACtaacaatgccggacaggagtactatcag                     <  1:2983084‑M1/49‑29 (MQ=255)
                            ataACAACATCAATGAATACtaacaatgccggacaggagtactatcagg                    <  1:1778488‑M1/49‑30 (MQ=255)
                            ataACAACATCAATGAATACtaacaatgccggacaggagtactatcagg                    >  1:3267886‑M1/1‑20 (MQ=255)
                            ataACAACATCAATGAATACtaacaatgccggacaggagtactatcagg                    <  1:3064107‑M1/49‑30 (MQ=255)
                             taACAACATCAATGAATACtaacaatgccggacaggagtactatcaggca                  <  1:19818‑M1/50‑32 (MQ=255)
                             taACAACATCAATGAATACtaacaatgccggacaggagtactatcaggc                   <  1:2223887‑M1/49‑31 (MQ=255)
                               acaacaTCAATGAATACtaacaatgccggacaggagtactatcaggcatt                <  1:408935‑M1/50‑34 (MQ=255)
                                  acaTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgt             >  1:3449887‑M1/1‑14 (MQ=255)
                                    aTCAATGAATACtaacaatgccggacaggagtactatcaggcattcgttt           >  1:3732348‑M1/1‑12 (MQ=255)
                                     tCAATGAATACtaacaatgccggacaggagtactatcaggcattcgtttt          <  1:2091585‑M1/50‑40 (MQ=255)
                                      cAATGAATACtaacaatgccggacaggagtactatcaggcattcgtttta         >  1:719204‑M1/1‑10 (MQ=255)
                                       aaTGAATACtaacaatgccggacaggagtactatcaggcattcgttttat        <  1:520336‑M1/50‑42 (MQ=255)
                                          gAATACtaacaatgccggacaggagtactatcaggcattcgttttatac      >  1:568627‑M1/1‑6 (MQ=255)
                                           aaTACtaacaatgccggacaggagtactatcaggcattcgttttatacga    <  1:3564087‑M1/50‑46 (MQ=255)
                                           aaTACtaacaatgccggacaggagtactatcaggcattcgttttatacg     <  1:790935‑M1/49‑45 (MQ=255)
                                            aTACtaacaatgccggacaggagtactatcaggcattcgttttatacgat   >  1:3960502‑M1/1‑4 (MQ=255)
                                             tACtaacaatgccggacaggagtactatcaggcattcgttttatacgatc  <  1:3499530‑M1/50‑48 (MQ=255)
                                               |                                               
GCTCAGACACATTACATTCACTGTATTTATAACAACATCAATGAATACCTTAGCCACAGGGATGGTTCCTGAAACATCAGTATTGCTGGTTGATG  >  NC_000913/3362114‑3362208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: