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

TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACG  >  NC_000913/2336584‑2336681
                                                 |                                                
tgatggcaatcctgatggtcattgtcccaacgcttactattggttatatt                                                  >  1:229115‑M2/50‑50 (MQ=255)
  atggcaatcctgatggtcattgtcccaacgcttactattggttatattt                                                 >  1:1372820‑M2/48‑49 (MQ=255)
     gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcacc                                             <  1:2587514‑M2/6‑1 (MQ=255)
     gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                              <  1:857824‑M2/5‑1 (MQ=255)
     gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                              <  1:263570‑M2/5‑1 (MQ=255)
       aatcctgatggtcattgtcccaacgcttactattggttatatTTcaccac                                           >  1:886582‑M2/43‑50 (MQ=255)
        atcctgatggtcattgtcccaacgcttactattggttatatTTCACCACt                                          >  1:1188184‑M2/42‑50 (MQ=255)
         tcctgatggtcattgtcccaacgcttactattggttatatTTCACCACTc                                         <  1:1991023‑M2/10‑1 (MQ=255)
           ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCa                                        <  1:1636732‑M2/11‑1 (MQ=255)
           ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAt                                       >  1:1805944‑M2/39‑50 (MQ=255)
             gatggtcattgtcccaacgcttactattggttatatTTCACCACTCATTg                                     >  1:2192454‑M2/37‑50 (MQ=255)
               tggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGcc                                   <  1:787367‑M2/16‑1 (MQ=255)
               tggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGcc                                   <  1:2045477‑M2/16‑1 (MQ=255)
               tggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                    >  1:553873‑M2/35‑49 (MQ=255)
                 gtcattgtcccaacgcttactattggttatatTTCACCACTCATTGCca                                  <  1:1046424‑M2/17‑1 (MQ=255)
                  tcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcaca                                <  1:485173‑M2/19‑1 (MQ=255)
                  tcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcaca                                >  1:2563211‑M2/32‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCGTTGCCACAt                               <  1:322172‑M2/20‑1 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              <  1:1426407‑M2/21‑1 (MQ=255)
                      tgtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTcc                            >  1:1900236‑M2/28‑50 (MQ=255)
                      tgtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTcc                            >  1:397114‑M2/28‑50 (MQ=255)
                      tgtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTcc                            >  1:2657093‑M2/28‑50 (MQ=255)
                       gtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTcc                            <  1:800079‑M2/23‑1 (MQ=255)
                       gtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTcc                            <  1:2681716‑M2/23‑1 (MQ=255)
                        tcccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCtt                          <  1:1473488‑M2/25‑1 (MQ=255)
                         cccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtg                         >  1:771440‑M2/25‑50 (MQ=255)
                          ccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgt                        <  1:2586102‑M2/27‑1 (MQ=255)
                          ccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtg                         <  1:464472‑M2/26‑1 (MQ=255)
                            aacgcttactattggttatatTTCACCACTCATTGCCACATTCCTTgtgt                      <  1:1893687‑M2/29‑1 (MQ=255)
                                cttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATag                  >  1:727032‑M2/18‑50 (MQ=255)
                                 ttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                 >  1:946307‑M2/17‑50 (MQ=255)
                                  tactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                 >  1:1982373‑M2/16‑49 (MQ=255)
                                   actattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagcc                <  1:530599‑M2/35‑1 (MQ=255)
                                   actattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                 >  1:1149721‑M2/15‑48 (MQ=255)
                                      attggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc             >  1:1592019‑M2/12‑49 (MQ=255)
                                      attggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc             >  1:2566987‑M2/12‑49 (MQ=255)
                                      attggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc             >  1:900671‑M2/12‑49 (MQ=255)
                                        tggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAt          >  1:1003320‑M2/10‑50 (MQ=255)
                                         ggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtt         >  1:332030‑M2/9‑50 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       <  1:531810‑M2/44‑1 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       >  1:2108142‑M2/7‑50 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       <  1:2060637‑M2/44‑1 (MQ=255)
                                            tatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt      >  1:283294‑M2/6‑50 (MQ=255)
                                            tatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt      >  1:2550709‑M2/6‑50 (MQ=255)
                                             atatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt      <  1:1946362‑M2/45‑1 (MQ=255)
                                              tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttttt     <  1:476434‑M2/46‑1 (MQ=255)
                                                tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc   <  1:1422319‑M2/48‑1 (MQ=255)
                                                tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc   <  1:107214‑M2/48‑1 (MQ=255)
                                                tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACg  >  1:181522‑M2/2‑50 (MQ=255)
                                                 |                                                
TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACG  >  NC_000913/2336584‑2336681

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: