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 47 [1] [1] 48 [relA]–[gudD] [relA], rlmD, barA, [gudD]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 29114901 (0.020)45 (0.900) 43/280 0.3 97.9% 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

ACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTACCA  >  NC_000913/2911341‑2911640
                                                                                                                                                     |                                                                                                                                                      
aCCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAg                                                                                                                                                       <  2:221658‑M1/151‑2 (MQ=255)
      tCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatcca                                                                                                                                                  <  1:155192‑M1/150‑7 (MQ=255)
        cTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagc                                                                                                                                                >  1:309007‑M1/1‑142 (MQ=255)
         tCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccggcat                                                                                                                                              >  1:750150‑M1/1‑141 (MQ=255)
          cAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcatt                                                                                                                                             <  2:121275‑M1/151‑12 (MQ=255)
              ccGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccg                                                                                                                                         >  2:40953‑M1/1‑136 (MQ=255)
                gCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgct                                                                                                                                       >  1:202487‑M1/1‑134 (MQ=255)
                gCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgct                                                                                                                                       <  2:483626‑M1/151‑18 (MQ=255)
                   aTCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttct                                                                                                                                    >  1:53523‑M1/1‑131 (MQ=255)
                   aTCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttc                                                                                                                                     <  2:691066‑M1/150‑20 (MQ=255)
                      atatGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatc                                                                                                                                 >  2:563344‑M1/1‑128 (MQ=255)
                        atGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatccc                                                                                                                               <  2:691399‑M1/151‑26 (MQ=255)
                          gTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccg                                                                                                                              <  2:810623‑M1/150‑27 (MQ=255)
                                aCAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaact                                                                                                                        <  2:13629‑M1/150‑33 (MQ=255)
                                  aGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactac                                                                                                                      >  2:51723‑M1/1‑116 (MQ=255)
                                   gTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacat                                                                                                                    <  1:390467‑M1/151‑37 (MQ=255)
                                     atatatCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatg                                                                                                                   <  2:620936‑M1/150‑38 (MQ=255)
                                           cAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtg                                                                                                            >  1:272838‑M1/1‑107 (MQ=255)
                                             aTCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggt                                                                                                          >  2:50958‑M1/1‑105 (MQ=255)
                                                tACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacg                                                                                                       >  2:371701‑M1/1‑102 (MQ=255)
                                                 aCATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggt                                                                                                          <  1:235571‑M1/147‑47 (MQ=255)
                                                  cATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgt                                                                                                      >  1:1790‑M1/1‑100 (MQ=255)
                                                     tGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtagg                                                                                                   <  1:535428‑M1/150‑54 (MQ=255)
                                                        aGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaa                                                                                               <  1:488890‑M1/151‑58 (MQ=255)
                                                                gCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCACCACGCGCCCCAGCAgatccagctttgccccttctttcccggtaactacatgaatagtggtacgtaggtcaattgtctg                                                                                        >  2:230085‑M1/1‑86 (MQ=255)
                                                                  aaaTTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaa                                                                                     <  2:188913‑M1/151‑68 (MQ=255)
                                                                        cGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcg                                                                               <  1:454743‑M1/151‑74 (MQ=255)
                                                                              aaCTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccacc                                                                         <  2:635900‑M1/151‑80 (MQ=255)
                                                                                cTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgccca                                                                           >  2:256727‑M1/1‑70 (MQ=255)
                                                                                 tCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagc                                                                      >  1:136500‑M1/1‑69 (MQ=255)
                                                                                  cccGTGCACCCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgccctccag                                                                       <  1:417047‑M1/149‑82 (MQ=255)
                                                                                       gCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacg                                                                <  2:510220‑M1/151‑89 (MQ=255)
                                                                                       gCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattaccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacg                                                                <  1:133270‑M1/151‑89 (MQ=255)
                                                                                        cAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacg                                                                <  2:228628‑M1/150‑89 (MQ=255)
                                                                                           ccGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcg                                                            <  2:531808‑M1/151‑93 (MQ=255)
                                                                                             gACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggc                                                          >  2:580675‑M1/1‑57 (MQ=255)
                                                                                                  gcgTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaa                                                       >  1:259960‑M1/1‑52 (MQ=255)
                                                                                                    gTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagt                                                    >  1:373950‑M1/1‑50 (MQ=255)
                                                                                                      cGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtatt                                                 <  1:557769‑M1/151‑104 (MQ=255)
                                                                                                        aTAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattac                                               >  2:505564‑M1/1‑46 (MQ=255)
                                                                                                                    aCCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                                    >  2:41891‑M1/1‑34 (MQ=255)
                                                                                                                           tGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtt                            >  1:833936‑M1/1‑27 (MQ=255)
                                                                                                                                ttACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgt                       >  2:45417‑M1/1‑22 (MQ=255)
                                                                                                                                   cccAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtatt                    >  2:502649‑M1/1‑19 (MQ=255)
                                                                                                                                         aCGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcaccc               >  1:699305‑M1/1‑13 (MQ=255)
                                                                                                                                               cccAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAg        <  2:687228/151‑1 (MQ=255)
                                                                                                                                                     agatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccagcgttttaccta  <  1:169465‑M1/151‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTACCA  >  NC_000913/2911341‑2911640

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ 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.