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] [0] 49 [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)48 (1.160) 45/280 0.1 100% coding (2162/2235 nt) relA (p)ppGpp synthetase I/GTP pyrophosphokinase
?NC_000913 2919019 = 0 (0.000)coding (367/1341 nt) gudD D‑glucarate dehydratase 1

CCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTACCA  >  NC_000913/2911342‑2911640
                                                                                                                                                    |                                                                                                                                                      
cccTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAga                                                                                                                                                      <  1:387929‑M1/151‑3 (MQ=255)
     tCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatcca                                                                                                                                                  <  2:622356‑M1/150‑7 (MQ=255)
       cTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagc                                                                                                                                                >  2:430608‑M1/1‑142 (MQ=255)
         cAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcatt                                                                                                                                             <  1:628131‑M1/151‑12 (MQ=255)
         cAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcatt                                                                                                                                             <  2:468220‑M1/151‑12 (MQ=255)
         cAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcatt                                                                                                                                             <  1:560791‑M1/151‑12 (MQ=255)
          aaaCCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattg                                                                                                                                            >  1:382130‑M1/1‑139 (MQ=255)
             ccGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccg                                                                                                                                         <  1:273816‑M1/151‑16 (MQ=255)
                     atatGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatc                                                                                                                                 >  2:648236‑M1/1‑128 (MQ=255)
                       atGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatccc                                                                                                                               <  2:291549‑M1/151‑26 (MQ=255)
                          tAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggt                                                                                                                            >  2:591127‑M1/1‑123 (MQ=255)
                           aGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggt                                                                                                                            >  2:198192‑M1/1‑122 (MQ=255)
                            gATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaa                                                                                                                          <  2:92210‑M1/151‑31 (MQ=255)
                                    atatatCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGGGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatg                                                                                                                   <  1:49233‑M1/150‑38 (MQ=255)
                                           aaTCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtgg                                                                                                           <  2:118223‑M1/151‑46 (MQ=255)
                                             tCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggt                                                                                                          >  2:428986‑M1/1‑104 (MQ=255)
                                                aCATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGACACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgt                                                                                                      <  1:519333‑M1/151‑51 (MQ=255)
                                                 cATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgt                                                                                                      >  2:555936‑M1/1‑100 (MQ=255)
                                                    tGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggt                                                                                                  >  1:116822‑M1/1‑97 (MQ=255)
                                                        gATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaa                                                                                               <  1:195030‑M1/150‑58 (MQ=255)
                                                         aTACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaat                                                                                             >  2:217541‑M1/1‑92 (MQ=255)
                                                             gAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtct                                                                                         >  2:678742‑M1/1‑88 (MQ=255)
                                                              aGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctg                                                                                        >  2:302914‑M1/1‑87 (MQ=255)
                                                                     ttCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaacc                                                                                  <  2:324224‑M1/150‑71 (MQ=255)
                                                                      tCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccg                                                                                 <  2:556055‑M1/150‑72 (MQ=255)
                                                                        ggCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgc                                                                              <  2:288554‑M1/151‑75 (MQ=255)
                                                                        ggCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgcagcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgc                                                                              >  2:598642‑M1/1‑77 (MQ=255)
                                                                          ccTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgccc                                                                            >  2:538276‑M1/1‑75 (MQ=255)
                                                                           cTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgccc                                                                            >  1:311969‑M1/1‑74 (MQ=255)
                                                                             aaCTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccacc                                                                         <  1:428949‑M1/151‑80 (MQ=255)
                                                                             aaCTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccac                                                                          <  1:78353‑M1/150‑79 (MQ=255)
                                                                                               gcgcgTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaa                                                        >  2:291048‑M1/1‑54 (MQ=255)
                                                                                                  cgTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagt                                                    <  2:586119‑M1/151‑101 (MQ=255)
                                                                                                   gTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagta                                                   <  2:195699‑M1/151‑102 (MQ=255)
                                                                                                   gTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaag                                                     >  2:372101‑M1/1‑50 (MQ=255)
                                                                                                     cGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtatt                                                 >  2:95177‑M1/1‑48 (MQ=255)
                                                                                                      gATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtatta                                                <  1:158481‑M1/151‑105 (MQ=255)
                                                                                                       aTAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgccccccagcatcacgatcggcaaaagtattac                                               <  1:346689‑M1/151‑106 (MQ=255)
                                                                                                        tAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattac                                               >  2:551465‑M1/1‑45 (MQ=255)
                                                                                                          aCATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgc                                             >  1:495710‑M1/1‑43 (MQ=255)
                                                                                                                 gCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgt                                     <  2:275047‑M1/151‑116 (MQ=255)
                                                                                                                 gCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgt                                     >  1:420776‑M1/1‑36 (MQ=255)
                                                                                                                    ccTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                                    >  1:647912‑M1/1‑33 (MQ=255)
                                                                                                                           gAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtt                            <  1:591086‑M1/150‑125 (MQ=255)
                                                                                                                                tACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgt                       >  2:33251‑M1/1‑21 (MQ=255)
                                                                                                                                tACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgt                       <  2:659525‑M1/150‑130 (MQ=255)
                                                                                                                                         cgcgCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccag             <  2:228431‑M1/151‑140 (MQ=255)
                                                                                                                                         cgcgCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccag             <  1:389741‑M1/151‑140 (MQ=255)
                                                                                                                                                    agatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccagcgttttaccta  <  2:431619‑M1/151‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTACCA  >  NC_000913/2911342‑2911640

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