Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,911,491 Δ7,528 bp [relA][gudD] [relA], rlmD, barA, [gudD]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2911491 2919018 7528 49 [0] [1] 50 [relA]–[gudD] [relA], rlmD, barA, [gudD]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 29114900 (0.000)49 (1.030) 43/280 0.2 99.0% coding (2162/2235 nt) relA (p)ppGpp synthetase I/GTP pyrophosphokinase
?NC_000913 2919019 = 1 (0.020)coding (367/1341 nt) gudD D‑glucarate dehydratase 1

CCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATC  >  NC_000913/2911342‑2911629
                                                                                                                                                    |                                                                                                                                           
cccTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAga                                                                                                                                           <  1:307502‑M1/151‑3 (MQ=255)
         cAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcat                                                                                                                                   <  2:405347‑M1/150‑11 (MQ=255)
          aaaCCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattg                                                                                                                                 >  2:296941‑M1/1‑139 (MQ=255)
            aCCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgcc                                                                                                                               <  2:352862‑M1/151‑15 (MQ=255)
            aCCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgcc                                                                                                                               <  2:546842‑M1/151‑15 (MQ=255)
             ccGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccg                                                                                                                              >  1:375961‑M1/1‑136 (MQ=255)
                  aTCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttct                                                                                                                         <  2:405573‑M1/151‑21 (MQ=255)
                   tCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttct                                                                                                                         >  1:68203‑M1/1‑130 (MQ=255)
                    cATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTTCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttcta                                                                                                                        <  2:793766‑M1/150‑22 (MQ=255)
                    cATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctat                                                                                                                       >  1:657637‑M1/1‑129 (MQ=255)
                       atGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatccc                                                                                                                    <  1:236921‑M1/151‑26 (MQ=255)
                                 aGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactac                                                                                                           >  2:579328‑M1/1‑116 (MQ=255)
                                              cTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTAGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtac                                                                                             >  1:556772‑M1/1‑103 (MQ=255)
                                                 cATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgt                                                                                           >  1:230031‑M1/1‑100 (MQ=255)
                                                 cATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgt                                                                                           >  1:38827‑M1/1‑100 (MQ=255)
                                                   ttGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtagg                                                                                        >  2:273534‑M1/1‑98 (MQ=255)
                                                    tGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggt                                                                                       >  1:479999‑M1/1‑97 (MQ=255)
                                                      tAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtc                                                                                      >  2:551393‑M1/1‑95 (MQ=255)
                                                       aGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtca                                                                                     >  2:621925‑M1/1‑94 (MQ=255)
                                                           aCGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgt                                                                                <  1:208384‑M1/151‑62 (MQ=255)
                                                            cGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtc                                                                               <  1:163112‑M1/151‑63 (MQ=255)
                                                             gAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtct                                                                              >  2:353604‑M1/1‑88 (MQ=255)
                                                             gAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtct                                                                              <  1:178601‑M1/151‑64 (MQ=255)
                                                               gCAAATTTCGGCCGAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGTCGCCTGGTTGAGTTTTCCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctac                                                                            <  2:68348‑M1/151‑66 (MQ=255)
                                                                       cGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcg                                                                    >  2:732352‑M1/1‑78 (MQ=255)
                                                                        ggCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgc                                                                   >  2:771732‑M1/1‑77 (MQ=255)
                                                                            tAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccac                                                               >  2:321892‑M1/1‑73 (MQ=255)
                                                                            tAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccac                                                               >  1:1251‑M1/1‑73 (MQ=255)
                                                                             aaCTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccacc                                                              <  2:1674‑M1/151‑80 (MQ=255)
                                                                                  ccGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcat                                                         >  1:149862‑M1/1‑67 (MQ=255)
                                                                                    gTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatca                                                       <  2:153586‑M1/151‑87 (MQ=255)
                                                                                      gCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacg                                                     <  1:530568‑M1/151‑89 (MQ=255)
                                                                                        aaCCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgat                                                   >  1:559813‑M1/1‑61 (MQ=255)
                                                                                        aaCCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacga                                                    >  2:250509‑M1/1‑61 (MQ=255)
                                                                                        aaCCGACGCGCGTCCATAACATCCGTCACCTGGTTGAGTTTGCCCATCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtatctacaggaatagtggtacgtaggtcaaatgtcgggaaaccgcgcccacccgcatcacg                                                     <  1:689968‑M1/149‑89 (MQ=255)
                                                                                           cGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcgg                                                <  1:51085‑M1/151‑94 (MQ=255)
                                                                                               gcgcgTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaa                                             >  2:235845‑M1/1‑54 (MQ=255)
                                                                                                   gTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagta                                        <  2:571544‑M1/151‑102 (MQ=255)
                                                                                                         aaCATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacg                                   <  1:423214‑M1/150‑107 (MQ=255)
                                                                                                              ccGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcacca                              <  2:11162‑M1/150‑112 (MQ=255)
                                                                                                                    ccTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtca                        <  2:500412‑M1/150‑118 (MQ=255)
                                                                                                                    ccTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                         >  1:502099‑M1/1‑33 (MQ=255)
                                                                                                                     cTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                         >  1:751729‑M1/1‑32 (MQ=255)
                                                                                                                      tGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcaga                      >  1:631358‑M1/1‑31 (MQ=255)
                                                                                                                       ggTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaa                     <  1:790616‑M1/150‑121 (MQ=255)
                                                                                                                         ttGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacg                   <  1:333073‑M1/150‑123 (MQ=255)
                                                                                                                                       cACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcaccc    >  2:134667‑M1/1‑14 (MQ=255)
                                                                                                                                        aCGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcaccc    >  1:456408‑M1/1‑13 (MQ=255)
                                                                                                                                         cgcgCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccag  <  2:230066‑M1/151‑140 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATC  >  NC_000913/2911342‑2911629

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