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 49 [0] [0] 49 [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)46 (1.550) 37/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

TGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATCGTAGAAACGGAAGGACGTTCAGCAGTCTTATCTGAAAAAGAGAAAAAA  >  NC_000913/2320095‑2320192
                                                |                                                 
tGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATt                                                  >  1:229115‑M1/1‑49 (MQ=255)
  aTGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATtt                                                 >  1:1372820‑M1/1‑47 (MQ=255)
     gCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcacc                                             <  1:2587514‑M1/50‑7 (MQ=255)
     gCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcac                                              <  1:857824‑M1/49‑6 (MQ=255)
     gCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcac                                              <  1:263570‑M1/49‑6 (MQ=255)
       aaTCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccac                                           >  1:886582‑M1/1‑42 (MQ=255)
        aTCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccact                                          >  1:1188184‑M1/1‑41 (MQ=255)
         tCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactc                                         <  1:1991023‑M1/50‑11 (MQ=255)
           cTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcat                                       >  1:1805944‑M1/1‑38 (MQ=255)
           cTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactca                                        <  1:1636732‑M1/49‑12 (MQ=255)
             gATGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattg                                     >  1:2192454‑M1/1‑36 (MQ=255)
               tGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgcc                                   <  1:787367‑M1/50‑17 (MQ=255)
               tGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgcc                                   <  1:2045477‑M1/50‑17 (MQ=255)
               tGGTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgc                                    >  1:553873‑M1/1‑34 (MQ=255)
                 gTCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgcca                                  <  1:1046424‑M1/49‑18 (MQ=255)
                  tCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccaca                                <  1:485173‑M1/50‑20 (MQ=255)
                  tCATTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccaca                                >  1:2563211‑M1/1‑31 (MQ=255)
                    aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcgttgccacat                               <  1:322172‑M1/49‑21 (MQ=255)
                    aTTGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacatt                              <  1:1426407‑M1/50‑22 (MQ=255)
                      tGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcc                            >  1:1900236‑M1/1‑27 (MQ=255)
                      tGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcc                            >  1:397114‑M1/1‑27 (MQ=255)
                      tGTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcc                            >  1:2657093‑M1/1‑27 (MQ=255)
                       gTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcc                            <  1:800079‑M1/49‑24 (MQ=255)
                       gTCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcc                            <  1:2681716‑M1/49‑24 (MQ=255)
                        tCCCAACGCTTACTATTGGTTATATttcaccactcattgccacattcctt                          <  1:1473488‑M1/50‑26 (MQ=255)
                         cccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttg                         >  1:771440‑M1/1‑24 (MQ=255)
                          ccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttgt                        <  1:2586102‑M1/50‑28 (MQ=255)
                          ccAACGCTTACTATTGGTTATATttcaccactcattgccacattccttg                         <  1:464472‑M1/49‑27 (MQ=255)
                            aaCGCTTACTATTGGTTATATttcaccactcattgccacattccttgtgt                      <  1:1893687‑M1/50‑30 (MQ=255)
                                cTTACTATTGGTTATATttcaccactcattgccacattccttgtgtatag                  >  1:727032‑M1/1‑17 (MQ=255)
                                 ttACTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 >  1:946307‑M1/1‑16 (MQ=255)
                                  tACTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 >  1:1982373‑M1/1‑15 (MQ=255)
                                   aCTATTGGTTATATttcaccactcattgccacattccttgtgtatagcc                <  1:530599‑M1/49‑36 (MQ=255)
                                   aCTATTGGTTATATttcaccactcattgccacattccttgtgtatagc                 >  1:1149721‑M1/1‑14 (MQ=255)
                                      aTTGGTTATATttcaccactcattgccacattccttgtgtatagccagc             >  1:1592019‑M1/1‑11 (MQ=255)
                                      aTTGGTTATATttcaccactcattgccacattccttgtgtatagccagc             >  1:2566987‑M1/1‑11 (MQ=255)
                                      aTTGGTTATATttcaccactcattgccacattccttgtgtatagccagc             >  1:900671‑M1/1‑11 (MQ=255)
                                        tGGTTATATttcaccactcattgccacattccttgtgtatagccagccat          >  1:1003320‑M1/1‑9 (MQ=255)
                                         ggTTATATttcaccactcattgccacattccttgtgtatagccagccatt         >  1:332030‑M1/1‑8 (MQ=255)
                                           ttATATttcaccactcattgccacattccttgtgtatagccagccatttt       <  1:531810‑M1/50‑45 (MQ=255)
                                           ttATATttcaccactcattgccacattccttgtgtatagccagccatttt       >  1:2108142‑M1/1‑6 (MQ=255)
                                           ttATATttcaccactcattgccacattccttgtgtatagccagccatttt       <  1:2060637‑M1/50‑45 (MQ=255)
                                            tataTttcaccactcattgccacattccttgtgtatagccagccattttt      >  1:283294‑M1/1‑5 (MQ=255)
                                            tataTttcaccactcattgccacattccttgtgtatagccagccattttt      >  1:2550709‑M1/1‑5 (MQ=255)
                                             atatttcaccactcattgccacattccttgtgtatagccagccattttt      <  1:1946362‑M1/49‑46 (MQ=255)
                                              tatttcaccactcattgccacattccttgtgtatagccagccatttttt     <  1:476434‑M1/49‑47 (MQ=255)
                                                tttcaccactcattgccacattccttgtgtatagccagccattttttacg  >  1:181522‑M1/1‑1 (MQ=255)
                                                tttcaccactcattgccacattccttgtgtatagccagccattttttac   <  1:1422319‑M1/49‑49 (MQ=255)
                                                tttcaccactcattgccacattccttgtgtatagccagccattttttac   <  1:107214‑M1/49‑49 (MQ=255)
                                                |                                                 
TGATGGCAATCCTGATGGTCATTGTCCCAACGCTTACTATTGGTTATATCGTAGAAACGGAAGGACGTTCAGCAGTCTTATCTGAAAAAGAGAAAAAA  >  NC_000913/2320095‑2320192

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: 
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.