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

TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACG  >  NC_000913/2336584‑2336681
                                                 |                                                
tgatggcaatcctgatggtcattgtcccaacgcttactattggttatatt                                                  >  1:2849915‑M2/50‑50 (MQ=255)
 gatggcaatcctgatggtcattgtcccaacgcttactattggttatattt                                                 >  1:3209594‑M2/49‑50 (MQ=255)
  atggcaatcctgatggtcattgtcccaacgcttactattggttatatt                                                  <  1:1665589‑M2/1‑1 (MQ=255)
  atggcaatcctgatggtcattgtcccaacgcttactattggttatatTTc                                                <  1:3519215‑M2/3‑1 (MQ=255)
     gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                              <  1:1897955‑M2/5‑1 (MQ=255)
       aatcctgatggtcattgtccccacgcttactattggttatatTTcaccac                                           >  1:4017648‑M2/43‑50 (MQ=255)
        atcctgatggtcattgtcccaacgcttactattggttatatTTCACCACt                                          <  1:2463621‑M2/9‑1 (MQ=255)
         tcctgatggtcattgtcccaacgcttactattggttatatTTCACCACt                                          >  1:809120‑M2/41‑49 (MQ=255)
          cctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCa                                        >  1:2509779‑M2/40‑50 (MQ=255)
           ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAt                                       <  1:2379080‑M2/12‑1 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCAtt                                      <  1:2167495‑M2/13‑1 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                    <  1:3227906‑M2/15‑1 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                    >  1:1783754‑M2/36‑50 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                    >  1:3580462‑M2/36‑50 (MQ=255)
                ggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGCca                                  >  1:2644426‑M2/34‑50 (MQ=255)
                 gtcattgtcccaacgcttactattggttatatTTCACCACTCATTGCca                                  >  1:1815640‑M2/33‑49 (MQ=255)
                 gtcattgtcccaacgcttactattggttatatTTCACCACTCATTGCca                                  <  1:2633733‑M2/17‑1 (MQ=255)
                  tcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcaca                                <  1:3451773‑M2/19‑1 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              >  1:3282321‑M2/30‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              >  1:2452782‑M2/30‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              >  1:4303243‑M2/30‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              >  1:103171‑M2/30‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAt                               <  1:1103994‑M2/20‑1 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAt                               <  1:321352‑M2/20‑1 (MQ=255)
                       gtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCt                           >  1:528819‑M2/27‑50 (MQ=255)
                        tcccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCtt                          >  1:3539058‑M2/26‑50 (MQ=255)
                         cccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtg                         >  1:3842894‑M2/25‑50 (MQ=255)
                         cccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtg                         >  1:3211678‑M2/25‑50 (MQ=255)
                          ccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgt                        >  1:3808373‑M2/24‑50 (MQ=255)
                          ccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtg                         <  1:1281923‑M2/26‑1 (MQ=255)
                          ccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtg                         <  1:353900‑M2/26‑1 (MQ=255)
                           caacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgtg                       <  1:1298055‑M2/28‑1 (MQ=255)
                            aacgctttctattggttatatTTCACCACTCATTGCCACATTCCTTgtgt                      <  1:2943122‑M2/29‑1 (MQ=255)
                            aacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgtg                       >  1:433386‑M2/22‑49 (MQ=255)
                            aacgcttactattggttatatTTCACCACTCATTGCCACATTCCTTgtgt                      <  1:543344‑M2/29‑1 (MQ=255)
                             acgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                     <  1:1522398‑M2/30‑1 (MQ=255)
                             acgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                     >  1:2284458‑M2/21‑50 (MQ=255)
                             acgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                     <  1:4003724‑M2/30‑1 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtat                    <  1:1821753‑M2/31‑1 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtat                    <  1:3706991‑M2/31‑1 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtat                    <  1:3099914‑M2/31‑1 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                     >  1:653127‑M2/20‑49 (MQ=255)
                               gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                     >  1:3576808‑M2/19‑48 (MQ=255)
                               gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                     >  1:2962714‑M2/19‑48 (MQ=255)
                                cttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATag                  >  1:1177580‑M2/18‑50 (MQ=255)
                                  tactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                 >  1:2868702‑M2/16‑49 (MQ=255)
                                  tactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                 >  1:1848316‑M2/16‑49 (MQ=255)
                                    ctattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccag              >  1:1651412‑M2/14‑50 (MQ=255)
                                         ggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAt          <  1:4039452‑M2/41‑1 (MQ=255)
                                          gttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttt        >  1:1272411‑M2/8‑50 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       >  1:660328‑M2/7‑50 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       >  1:3528440‑M2/7‑50 (MQ=255)
                                            tatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt      >  1:749944‑M2/6‑50 (MQ=255)
                                              tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttttt     <  1:4226236‑M2/46‑1 (MQ=255)
                                              tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTa    <  1:2898683‑M2/47‑1 (MQ=255)
                                               atTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc   >  1:1281301‑M2/3‑50 (MQ=255)
                                               atTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc   >  1:3891344‑M2/3‑50 (MQ=255)
                                               atTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc   >  1:3780908‑M2/3‑50 (MQ=255)
                                                tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc   <  1:973363‑M2/48‑1 (MQ=255)
                                                tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACg  <  1:3050221‑M2/49‑1 (MQ=255)
                                                 |                                                
TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACG  >  NC_000913/2336584‑2336681

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ 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.