Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,320,144 Δ16,489 bp [atoS][yfaA] [atoS], atoC, atoD, atoA, atoE, atoB, yfaP, yfaQ, yfaS, yfaT, [yfaA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2320144 2336632 16489 50 [0] [0] 50 [atoS]–[yfaA] [atoS],atoC,atoD,atoA,atoE,atoB,yfaP,yfaQ,yfaS,yfaT,[yfaA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 23201430 (0.000)45 (0.900) 34/96 0.1 100% coding (101/1827 nt) atoS sensory histidine kinase AtoS
?NC_000913 2336633 = 0 (0.000)coding (12/1689 nt) yfaA DUF2138 domain‑containing protein YfaA

CTGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATCGTAGAAACGGAAGGACGTTCAGCAGTCTTATCTGAAAAAGAGAAAAAA  >  NC_000913/2320094‑2320192
                                                 |                                                 
cTGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATAT                                                   >  1:3415464‑M1/1‑50 (MQ=255)
  gATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATtt                                                 >  1:2187750‑M1/1‑48 (MQ=255)
  gATGGCAATCCTGATGGTCATCGTCCCAACGCTTACTATTGGTTATATtt                                                 >  1:828103‑M1/1‑48 (MQ=255)
   aTGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATt                                                  <  1:3752490‑M1/48‑2 (MQ=255)
    tGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttca                                               >  1:611364‑M1/1‑46 (MQ=255)
    tGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttca                                               >  1:3489815‑M1/1‑46 (MQ=255)
      gCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcac                                              <  1:1028858‑M1/49‑6 (MQ=255)
      gCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcac                                              <  1:1977404‑M1/49‑6 (MQ=255)
        aaTCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccac                                           >  1:1474662‑M1/1‑42 (MQ=255)
        aaTCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccac                                           >  1:2912697‑M1/1‑42 (MQ=255)
         aTCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccact                                          <  1:412991‑M1/50‑10 (MQ=255)
            cTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcat                                       >  1:2197278‑M1/1‑38 (MQ=255)
            cTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcat                                       >  1:2141730‑M1/1‑38 (MQ=255)
             tGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcatt                                      >  1:2534743‑M1/1‑37 (MQ=255)
               aTGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgc                                    >  1:1813642‑M1/1‑35 (MQ=255)
               aTGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcatt                                      <  1:1319771‑M1/48‑14 (MQ=255)
                   tCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccaca                                >  1:1586470‑M1/1‑31 (MQ=255)
                      ttGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattc                             <  1:1200983‑M1/50‑23 (MQ=255)
                        gTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcct                           >  1:2639375‑M1/1‑26 (MQ=255)
                        gTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcc                            <  1:363672‑M1/49‑24 (MQ=255)
                         tCCCAACGCTTACCATTGGTTATATttcaccactcattgccacattcct                           <  1:2117678‑M1/49‑25 (MQ=255)
                           ccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttgt                        <  1:2747992‑M1/50‑28 (MQ=255)
                            cAACGCTTACTATTGGTTATATttcaccactcattgccacattccttgtg                       >  1:1207553‑M1/1‑22 (MQ=255)
                            cAACGCTTACTATTGGTTATATttcaccactcattgccacattccttgtg                       >  1:2066928‑M1/1‑22 (MQ=255)
                              aCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     <  1:1570638‑M1/50‑31 (MQ=255)
                              aCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     >  1:1717121‑M1/1‑20 (MQ=255)
                                gCTTACTATTGGTTATATttcaccactcattgccacattccttgtgtata                   <  1:2679823‑M1/50‑33 (MQ=255)
                                gCTTACTATTGGTTATATttcaccactcattgccacattccttgtgtata                   <  1:3439132‑M1/50‑33 (MQ=255)
                                 cTTACTATTGGTTATATttcaccactcattgccacattccttgtgtatag                  >  1:1270014‑M1/1‑17 (MQ=255)
                                 cTTACTATTGGTTATATttcaccactcattgccacattccttgtgtatag                  <  1:296405‑M1/50‑34 (MQ=255)
                                 cTTACTATTGGTTATATttcaccactcattgccacattccttgtgtatag                  <  1:1774590‑M1/50‑34 (MQ=255)
                                  ttACTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 >  1:3003751‑M1/1‑16 (MQ=255)
                                    aCTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 <  1:2861732‑M1/48‑35 (MQ=255)
                                      tATTGGTTATATttcaccactcattgccacattccttgtgtatagccagc             >  1:3541224‑M1/1‑12 (MQ=255)
                                        ttGGTTATATttcaccactcattgccacattccttgtgtatagccagcca           <  1:2162467‑M1/50‑41 (MQ=255)
                                        ttGGTTATATttcaccactcattgccacattccttgtgtatagccagcca           <  1:2326735‑M1/50‑41 (MQ=255)
                                        ttGGTTATATttcaccactcattgccacattccttgtgtatagccagcc            <  1:3927243‑M1/49‑40 (MQ=255)
                                        ttGGTTATATttcaccactcattgccacattccttgtgtatagccagc             <  1:778488‑M1/48‑39 (MQ=255)
                                          ggTTATATttcaccactcattgccacattccttgtgtatagccagccatt         >  1:760967‑M1/1‑8 (MQ=255)
                                           gTTATATttcaccactcattgccacattccttgtgtatagccagccattt        >  1:3181257‑M1/1‑7 (MQ=255)
                                           gTTATATttcaccactcattgccacattccttgtgtatagccagccattt        >  1:3618934‑M1/1‑7 (MQ=255)
                                             tataTttcaccactcattgccacattccttgtgtatagccagccattttt      >  1:1856184‑M1/1‑5 (MQ=255)
                                              atatttcaccactcattgccacattccttgtgtatagccagccattttt      <  1:1085855‑M1/49‑46 (MQ=255)
                                               tatttcaccactcattgccacattccttgtgtatagccagccatttttta    >  1:614713‑M1/1‑3 (MQ=255)
                                                atttcaccactcattgccacattccttgtgtatagccagccattttttac   >  1:2988457‑M1/1‑2 (MQ=255)
                                                atttcaccactcattgccacattccttgtgtatagccagccattttttac   >  1:1269824‑M1/1‑2 (MQ=255)
                                                 tttcaccactcattgccacattccttgtgtatagccagccattttttacg  <  1:729787‑M1/50‑50 (MQ=255)
                                                 tttcaccactcattgccacattccttgtgtatagccagccattttttac   <  1:343571‑M1/49‑49 (MQ=255)
                                                 tttcaccactcattgccacattccttgtgtatagccagccattttttac   <  1:3934612‑M1/49‑49 (MQ=255)
                                                 tttcaccactcattgccacattccttgtgtatagccagccattttttac   <  1:651205‑M1/49‑49 (MQ=255)
                                                 |                                                 
CTGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATCGTAGAAACGGAAGGACGTTCAGCAGTCTTATCTGAAAAAGAGAAAAAA  >  NC_000913/2320094‑2320192

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: