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 26 [0] [0] 26 [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)25 (0.850) 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

TTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACG  >  NC_000913/2336584‑2336681
                                                 |                                                
tgatggcaatcctgatggtcattgtcccaacgcttactattggttatatt                                                  >  1:1030964‑M2/50‑50 (MQ=255)
  atggcaatcctgatggtcattgtcccaacgcttactattggttatatTTc                                                <  1:122991‑M2/3‑1 (MQ=255)
   tggcaatcctgatggtcattgtcccaacgcttactattggttatatTTca                                               <  1:502499‑M2/4‑1 (MQ=255)
      caatcctgatggtcattgtcccaacgcttactattggttatatTTcacca                                            <  1:927711‑M2/7‑1 (MQ=255)
        atcctgatggtcattgtcccaacgcttactattggttatatTTCACCACt                                          >  1:1241145‑M2/42‑50 (MQ=255)
        atcctgatggtcattgtcccaacgcttactattggttatatTTCACCACt                                          <  1:898545‑M2/9‑1 (MQ=255)
              atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                    >  1:1410805‑M2/36‑50 (MQ=255)
               tggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                    >  1:1044494‑M2/35‑49 (MQ=255)
                   cattgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAt                               >  1:1238198‑M2/31‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              <  1:350851‑M2/21‑1 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                              >  1:1937373‑M2/30‑50 (MQ=255)
                    attgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAt                               <  1:1832855‑M2/20‑1 (MQ=255)
                       gtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCt                           >  1:1336708‑M2/27‑50 (MQ=255)
                            aacgcttactattggttatatTTCACCACTCATTGCCACATTCCTTgtgt                      <  1:773619‑M2/29‑1 (MQ=255)
                               gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtata                   <  1:1910261‑M2/32‑1 (MQ=255)
                               gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtata                   <  1:873405‑M2/32‑1 (MQ=255)
                                  tactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                 >  1:1487532‑M2/16‑49 (MQ=255)
                                     tattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc             >  1:1511780‑M2/13‑50 (MQ=255)
                                     tattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc             >  1:841343‑M2/13‑50 (MQ=255)
                                       ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAgccagcca           <  1:931153‑M2/40‑1 (MQ=255)
                                        tggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAt          >  1:1750810‑M2/10‑50 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       >  1:1950463‑M2/7‑50 (MQ=255)
                                           ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt       <  1:1341727‑M2/44‑1 (MQ=255)
                                             atatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt      >  1:975877‑M2/5‑49 (MQ=255)
                                              tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttttt     <  1:1119916‑M2/46‑1 (MQ=255)
                                                tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACg  >  1:1789034‑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: