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 54 [0] [0] 53 [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)50 (0.830) 46/280 0.5 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

ACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTA  >  NC_000913/2911341‑2911637
                                                                                                                                                     |                                                                                                                                                   
aCCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAg                                                                                                                                                    <  2:21388‑M1/151‑2 (MQ=255)
    tttCCTCAAACCGCTTTCATATGTAGCTCCAGTATATATCAATCTACATTGTAGACACGCGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCCCGTCGATAACATCAGGCACCTGGATGAGTTTACACAGCACGCGCCccaccagaacc                                                                                                                                                >  2:820602/1‑146 (MQ=255)
       ccTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccag                                                                                                                                              <  2:239909‑M1/150‑8 (MQ=255)
         tCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcat                                                                                                                                           >  1:948594‑M1/1‑141 (MQ=255)
             aCCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgcc                                                                                                                                       >  2:219976‑M1/1‑137 (MQ=255)
              ccGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccg                                                                                                                                      >  1:611689‑M1/1‑136 (MQ=255)
               cGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgc                                                                                                                                     >  1:409257‑M1/1‑135 (MQ=255)
                gCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgct                                                                                                                                    >  2:887722‑M1/1‑134 (MQ=255)
                  tATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttc                                                                                                                                  >  1:292409‑M1/1‑132 (MQ=255)
                   aTCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttc                                                                                                                                  <  2:290041‑M1/150‑20 (MQ=255)
                     cATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctat                                                                                                                               >  1:323045‑M1/1‑129 (MQ=255)
                        atttagaTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCCCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatccc                                                                                                                            >  2:700550‑M1/4‑126 (MQ=255)
                         tGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccg                                                                                                                           >  2:432279‑M1/1‑125 (MQ=255)
                           tAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggt                                                                                                                         >  2:110191‑M1/1‑123 (MQ=255)
                            aGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggta                                                                                                                        <  2:795219‑M1/151‑30 (MQ=255)
                                                tACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacg                                                                                                    >  1:272118‑M1/1‑102 (MQ=255)
                                                 aCATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacg                                                                                                    <  1:519967‑M1/150‑50 (MQ=255)
                                                 aCATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGACCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggta                                                                                                      <  2:932957‑M1/148‑48 (MQ=255)
                                                   aTTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgt                                                                                                   >  2:903154‑M1/1‑99 (MQ=255)
                                                     tGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggt                                                                                               >  1:346060‑M1/1‑97 (MQ=255)
                                                        aGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaa                                                                                            <  1:224448‑M1/151‑58 (MQ=255)
                                                        aGATACGAGCAAATTTCGGCCTAACTCCCGTACAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtatctacatgaatagtggtacgtaggtcaa                                                                                            <  2:885301‑M1/151‑58 (MQ=255)
                                                            aCGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgt                                                                                        >  1:1010528‑M1/1‑90 (MQ=255)
                                                                      ttCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaacc                                                                               >  1:409415‑M1/1‑80 (MQ=255)
                                                                       tCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccg                                                                              <  1:587854‑M1/150‑72 (MQ=255)
                                                                         ggCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgc                                                                           >  2:645338‑M1/1‑77 (MQ=255)
                                                                         ggCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgc                                                                           <  2:594646‑M1/151‑75 (MQ=255)
                                                                               aCTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccacca                                                                     <  1:588287‑M1/151‑81 (MQ=255)
                                                                                    cGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatc                                                                >  1:706223‑M1/1‑66 (MQ=255)
                                                                                    cGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatc                                                                <  2:48716‑M1/151‑86 (MQ=255)
                                                                                           ccGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcg                                                         <  2:581866‑M1/151‑93 (MQ=255)
                                                                                                  gcgTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaag                                                  <  2:557649‑M1/151‑100 (MQ=255)
                                                                                                      cGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtat                                               <  2:437801‑M1/150‑103 (MQ=255)
                                                                                                      cGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtat                                               <  2:650166‑M1/150‑103 (MQ=255)
                                                                                                          aaCATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgctgctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgc                                          <  1:370893‑M1/151‑108 (MQ=255)
                                                                                                           aCATCCGGCACCTGGTTGAGTTTACCCAGCGCGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccatcatcacgatcggcaaaagtattacgc                                          <  2:450301‑M1/150‑108 (MQ=255)
                                                                                                           aCATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgca                                         <  1:231851‑M1/151‑109 (MQ=255)
                                                                                                            cATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcac                                        <  1:337718‑M1/151‑110 (MQ=255)
                                                                                                             aTCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcacc                                       >  2:604692‑M1/1‑41 (MQ=255)
                                                                                                                   cACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                                 >  2:211265‑M1/1‑35 (MQ=255)
                                                                                                                   cACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                                 <  1:951357‑M1/151‑117 (MQ=255)
                                                                                                                    aCCTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtca                                <  1:761364‑M1/151‑118 (MQ=255)
                                                                                                                     ccTGGTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtc                                 >  2:878115‑M1/1‑33 (MQ=255)
                                                                                                                        ggTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcaga                              <  1:563926‑M1/149‑120 (MQ=255)
                                                                                                                        ggTTGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcagcaaaagtattacgcaccagcgtcagaa                             <  2:225349‑M1/150‑121 (MQ=255)
                                                                                                                           tGAGTTTACCCAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgt                          <  1:471191‑M1/150‑124 (MQ=255)
                                                                                                                                     cAGCACGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattca               <  2:524815‑M1/151‑135 (MQ=255)
                                                                                                                                         aCGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcaccca           <  2:566092‑M1/151‑139 (MQ=255)
                                                                                                                                         aCGCGCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcaccc            >  1:914187‑M1/1‑13 (MQ=255)
                                                                                                                                          cgcgCCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccag          <  2:95845‑M1/151‑140 (MQ=255)
                                                                                                                                            cgcCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccagcg        <  1:511101‑M1/151‑142 (MQ=255)
                                                                                                                                             gcCCCAGCAgatccagcattgccgcttctatcccggtaactacatgaatagaggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccagcgt       >  2:1019895‑M1/1‑9 (MQ=255)
                                                                                                                                                ccAGCAgatccagcattgccccttctatcccggtaactacatgaatagtggtacgtagatcaaatgtctgcaaaccgcgcccgccagcatcacgatcggcaaaagtattacgcagcagcgtcagaacgttattgtattcccccagcgtttt    <  2:107202‑M1/151‑146 (MQ=255)
                                                                                                                                                  agcagatccagcattgccgcttctatcccggtaactacatgaatagtggtacgtaggtcaaatgtctgcaaaccgcgcccaccagcatcacgatcggcaaaagtattacgcaccagcgtcagaacgtttttgtattcacccagcgttttac  <  1:939070‑M1/151‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTACGGCTGGCAACGCCAAGCACGTTCACCTTCTCGTTGGCGAGAATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTA  >  NC_000913/2911341‑2911637

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