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 50 [0] [0] 50 [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)45 (0.900) 34/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

TGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGG  >  NC_000913/2336585‑2336682
                                                |                                                 
gatggcaatcctgatggtcattgtcccaacgcttactattggttatattt                                                  >  1:2187750‑M2/49‑50 (MQ=255)
gatggcaatcctgatggtcatcgtcccaacgcttactattggttatattt                                                  >  1:828103‑M2/49‑50 (MQ=255)
 atggcaatcctgatggtcattgtcccaacgcttactattggttatatt                                                   <  1:3752490‑M2/1‑1 (MQ=255)
  tggcaatcctgatggtcattgtcccaacgcttactattggttatatTTca                                                >  1:3489815‑M2/47‑50 (MQ=255)
  tggcaatcctgatggtcattgtcccaacgcttactattggttatatTTca                                                >  1:611364‑M2/47‑50 (MQ=255)
    gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                               <  1:1977404‑M2/5‑1 (MQ=255)
    gcaatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                               <  1:1028858‑M2/5‑1 (MQ=255)
      aatcctgatggtcattgtcccaacgcttactattggttatatTTcaccac                                            >  1:1474662‑M2/43‑50 (MQ=255)
      aatcctgatggtcattgtcccaacgcttactattggttatatTTcaccac                                            >  1:2912697‑M2/43‑50 (MQ=255)
       atcctgatggtcattgtcccaacgcttactattggttatatTTCACCACt                                           <  1:412991‑M2/9‑1 (MQ=255)
          ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAt                                        >  1:2197278‑M2/39‑50 (MQ=255)
          ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAt                                        >  1:2141730‑M2/39‑50 (MQ=255)
           tgatggtcattgtcccaacgcttactattggttatatTTCACCACTCAtt                                       >  1:2534743‑M2/38‑50 (MQ=255)
             atggtcattgtcccaacgcttactattggttatatTTCACCACTCAtt                                       <  1:1319771‑M2/13‑1 (MQ=255)
             atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                     >  1:1813642‑M2/36‑50 (MQ=255)
                 tcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcaca                                 >  1:1586470‑M2/32‑50 (MQ=255)
                    ttgtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTc                              <  1:1200983‑M2/22‑1 (MQ=255)
                      gtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTcc                             <  1:363672‑M2/23‑1 (MQ=255)
                      gtcccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCt                            >  1:2639375‑M2/27‑50 (MQ=255)
                       tcccaacgcttaccattggttatatTTCACCACTCATTGCCACATTCCt                            <  1:2117678‑M2/24‑1 (MQ=255)
                         ccaacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgt                         <  1:2747992‑M2/27‑1 (MQ=255)
                          caacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgtg                        >  1:2066928‑M2/23‑50 (MQ=255)
                          caacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgtg                        >  1:1207553‑M2/23‑50 (MQ=255)
                            acgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                      <  1:1570638‑M2/30‑1 (MQ=255)
                            acgcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                      >  1:1717121‑M2/21‑50 (MQ=255)
                              gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtata                    <  1:2679823‑M2/32‑1 (MQ=255)
                              gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGtata                    <  1:3439132‑M2/32‑1 (MQ=255)
                               cttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATag                   <  1:1774590‑M2/33‑1 (MQ=255)
                               cttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATag                   <  1:296405‑M2/33‑1 (MQ=255)
                               cttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATag                   >  1:1270014‑M2/18‑50 (MQ=255)
                                ttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                  >  1:3003751‑M2/17‑50 (MQ=255)
                                  actattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                  <  1:2861732‑M2/34‑1 (MQ=255)
                                    tattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc              >  1:3541224‑M2/13‑50 (MQ=255)
                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagcc             <  1:3927243‑M2/39‑1 (MQ=255)
                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc              <  1:778488‑M2/38‑1 (MQ=255)
                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAgccagcca            <  1:2162467‑M2/40‑1 (MQ=255)
                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAgccagcca            <  1:2326735‑M2/40‑1 (MQ=255)
                                        ggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtt          >  1:760967‑M2/9‑50 (MQ=255)
                                         gttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttt         >  1:3618934‑M2/8‑50 (MQ=255)
                                         gttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttt         >  1:3181257‑M2/8‑50 (MQ=255)
                                           tatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt       >  1:1856184‑M2/6‑50 (MQ=255)
                                            atatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt       <  1:1085855‑M2/45‑1 (MQ=255)
                                             tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTa     >  1:614713‑M2/4‑50 (MQ=255)
                                              atTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc    >  1:2988457‑M2/3‑50 (MQ=255)
                                              atTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc    >  1:1269824‑M2/3‑50 (MQ=255)
                                               tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc    <  1:343571‑M2/48‑1 (MQ=255)
                                               tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc    <  1:3934612‑M2/48‑1 (MQ=255)
                                               tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc    <  1:651205‑M2/48‑1 (MQ=255)
                                               tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACg   <  1:729787‑M2/49‑1 (MQ=255)
                                                ttCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACgg  >  1:3662171/1‑50 (MQ=255)
                                                |                                                 
TGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGG  >  NC_000913/2336585‑2336682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: