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 27 [0] [0] 27 [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)27 (1.010) 23/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

TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTT  >  NC_000913/2336584‑2336677
                                                 |                                            
tgatggcaatcctgatggtcattgtcccaacgcttactattggttatatt                                              >  1:1162845‑M2/50‑50 (MQ=255)
     gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                          >  1:1801298‑M2/45‑49 (MQ=255)
      caatcctgatggtcattgtcccaacgcttactattggttatatTTcacca                                        >  1:1238886‑M2/44‑50 (MQ=255)
      caatcctgatggtcattgtcccaacgcttactattggttatatTTcacca                                        >  1:2185805‑M2/44‑50 (MQ=255)
           ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAt                                   <  1:622902‑M2/12‑1 (MQ=255)
            tgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAtt                                  >  1:2187992‑M2/38‑50 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCAtt                                  <  1:1429113‑M2/13‑1 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                >  1:1999062‑M2/36‑50 (MQ=255)
                 gtcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcac                             <  1:640645‑M2/18‑1 (MQ=255)
                  tcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcaca                            >  1:396541‑M2/32‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                          >  1:1760663‑M2/30‑50 (MQ=255)
                        tcccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCtt                      <  1:1829262‑M2/25‑1 (MQ=255)
                            aacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgtg                   >  1:913359‑M2/22‑49 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtat                <  1:1143931‑M2/31‑1 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                 >  1:183157‑M2/20‑49 (MQ=255)
                              cgcttactattggttatatTTCACCACTCATTGCCAAATTCCTTGTGtat                <  1:1909645‑M2/31‑1 (MQ=255)
                                cttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATag              >  1:1969082‑M2/18‑50 (MQ=255)
                                  tactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc             >  1:2036934‑M2/16‑49 (MQ=255)
                                   actattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc             >  1:2186722‑M2/15‑48 (MQ=255)
                                   actattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc             >  1:716375‑M2/15‑48 (MQ=255)
                                    ctattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccag          >  1:2136529‑M2/14‑50 (MQ=255)
                                       ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc         <  1:826000‑M2/38‑1 (MQ=255)
                                       ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc         <  1:940464‑M2/38‑1 (MQ=255)
                                        tggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAt      >  1:354591‑M2/10‑50 (MQ=255)
                                         ggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAt      <  1:1848886‑M2/41‑1 (MQ=255)
                                            tatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  >  1:1194436‑M2/6‑50 (MQ=255)
                                             atatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt  <  1:1547840‑M2/45‑1 (MQ=255)
                                                 |                                            
TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTT  >  NC_000913/2336584‑2336677

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.