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 61 [0] [0] 60 [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)58 (1.020) 42/96 0.0 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:2399014‑M1/50‑1 (MQ=255)
 tGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATt                                                  >  1:2849915‑M1/1‑49 (MQ=255)
  gATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATtt                                                 >  1:3209594‑M1/1‑48 (MQ=255)
   aTGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttc                                                <  1:3519215‑M1/50‑4 (MQ=255)
   aTGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATt                                                  <  1:1665589‑M1/48‑2 (MQ=255)
      gCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcac                                              <  1:1897955‑M1/49‑6 (MQ=255)
        aaTCCTGATGGTCATTGTCCCCACGCTTACTATTGGTTATATttcaccac                                           >  1:4017648‑M1/1‑42 (MQ=255)
         aTCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccact                                          <  1:2463621‑M1/50‑10 (MQ=255)
          tCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccact                                          >  1:809120‑M1/1‑40 (MQ=255)
           ccTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactca                                        >  1:2509779‑M1/1‑39 (MQ=255)
            cTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcat                                       <  1:2379080‑M1/50‑13 (MQ=255)
               aTGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgc                                    >  1:3580462‑M1/1‑35 (MQ=255)
               aTGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgc                                    >  1:1783754‑M1/1‑35 (MQ=255)
               aTGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgc                                    <  1:3227906‑M1/50‑16 (MQ=255)
               aTGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcatt                                      <  1:2167495‑M1/48‑14 (MQ=255)
                 ggTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgcca                                  >  1:2644426‑M1/1‑33 (MQ=255)
                  gTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgcca                                  >  1:1815640‑M1/1‑32 (MQ=255)
                  gTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgcca                                  <  1:2633733‑M1/49‑18 (MQ=255)
                   tCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccaca                                <  1:3451773‑M1/50‑20 (MQ=255)
                     aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacatt                              >  1:4303243‑M1/1‑29 (MQ=255)
                     aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacatt                              >  1:103171‑M1/1‑29 (MQ=255)
                     aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacatt                              >  1:2452782‑M1/1‑29 (MQ=255)
                     aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacatt                              >  1:3282321‑M1/1‑29 (MQ=255)
                     aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacat                               <  1:321352‑M1/49‑21 (MQ=255)
                     aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacat                               <  1:1103994‑M1/49‑21 (MQ=255)
                        gTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcct                           >  1:528819‑M1/1‑26 (MQ=255)
                         tCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcctt                          >  1:3539058‑M1/1‑25 (MQ=255)
                          cccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttg                         >  1:3842894‑M1/1‑24 (MQ=255)
                          cccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttg                         >  1:3211678‑M1/1‑24 (MQ=255)
                           ccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttgt                        >  1:3808373‑M1/1‑23 (MQ=255)
                           ccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttg                         <  1:1281923‑M1/49‑27 (MQ=255)
                           ccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttg                         <  1:353900‑M1/49‑27 (MQ=255)
                            cAACGCTTACTATTGGTTATATttcaccactcattgccacattccttgtg                       <  1:1298055‑M1/50‑29 (MQ=255)
                             aaCGCTTTCTATTGGTTATATttcaccactcattgccacattccttgtgt                      <  1:2943122‑M1/50‑30 (MQ=255)
                             aaCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgt                      <  1:543344‑M1/50‑30 (MQ=255)
                             aaCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtg                       >  1:433386‑M1/1‑21 (MQ=255)
                              aCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     <  1:4003724‑M1/50‑31 (MQ=255)
                              aCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     >  1:2284458‑M1/1‑20 (MQ=255)
                              aCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     <  1:1522398‑M1/50‑31 (MQ=255)
                               cGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgtat                    <  1:3706991‑M1/50‑32 (MQ=255)
                               cGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgtat                    <  1:3099914‑M1/50‑32 (MQ=255)
                               cGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgtat                    <  1:1821753‑M1/50‑32 (MQ=255)
                               cGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     >  1:653127‑M1/1‑19 (MQ=255)
                                gCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     >  1:2962714‑M1/1‑18 (MQ=255)
                                gCTTACTATTGGTTATATttcaccactcattgccacattccttgtgta                     >  1:3576808‑M1/1‑18 (MQ=255)
                                 cTTACTATTGGTTATATttcaccactcattgccacattccttgtgtatag                  >  1:1177580‑M1/1‑17 (MQ=255)
                                   tACTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 >  1:2868702‑M1/1‑15 (MQ=255)
                                   tACTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 >  1:1848316‑M1/1‑15 (MQ=255)
                                     cTATTGGTTATATttcaccactcattgccacattccttgtgtatagccag              >  1:1651412‑M1/1‑13 (MQ=255)
                                          ggTTATATttcaccactcattgccacattccttgtgtatagccagccat          <  1:4039452‑M1/49‑42 (MQ=255)
                                           gTTATATttcaccactcattgccacattccttgtgtatagccagccattt        >  1:1272411‑M1/1‑7 (MQ=255)
                                            ttATATttcaccactcattgccacattccttgtgtatagccagccatttt       >  1:3528440‑M1/1‑6 (MQ=255)
                                            ttATATttcaccactcattgccacattccttgtgtatagccagccatttt       >  1:660328‑M1/1‑6 (MQ=255)
                                             tataTttcaccactcattgccacattccttgtgtatagccagccattttt      >  1:749944‑M1/1‑5 (MQ=255)
                                               tatttcaccactcattgccacattccttgtgtatagccagccatttttta    <  1:2898683‑M1/50‑48 (MQ=255)
                                               tatttcaccactcattgccacattccttgtgtatagccagccatttttt     <  1:4226236‑M1/49‑47 (MQ=255)
                                                atttcaccactcattgccacattccttgtgtatagccagccattttttac   >  1:1281301‑M1/1‑2 (MQ=255)
                                                atttcaccactcattgccacattccttgtgtatagccagccattttttac   >  1:3891344‑M1/1‑2 (MQ=255)
                                                atttcaccactcattgccacattccttgtgtatagccagccattttttac   >  1:3780908‑M1/1‑2 (MQ=255)
                                                 tttcaccactcattgccacattccttgtgtatagccagccattttttacg  <  1:3050221‑M1/50‑50 (MQ=255)
                                                 tttcaccactcattgccacattccttgtgtatagccagccattttttac   <  1:973363‑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 < 18 ≤ ATCG/ATCG < 33 ≤ 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.