New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1852976 =39 (1.320)4 (0.140) 4/236 NT 10.0% coding (2020/2241 nt) relA ATP:GTP 3'‑pyrophosphotransferase
?NC_002947 1853007 = 35 (1.240)coding (2051/2241 nt) relA ATP:GTP 3'‑pyrophosphotransferase

AGGACTGCGCCTCGGTACTGCAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCT  >  NC_002947/1852869‑1853150
                                                                                                                                          |                                                                                                                                               
aGGACTGCGCCTCGGTACTGCAGCTGGCGGGCAAAGAGCCAGAGCGCAGTATCCAGGTGAGCTGGGGGCGGATCCCGGTGCAGACCTCCCCGGTCGACATCGACATCCGAGCATACGACTGCCCCGTACGGCCGCGCGATGTGTCGCAgg                                                                                                                                      >  2:333438/1‑150 (MQ=255)
      gcgcCTCGGTATTTCAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGtgctgc                                                                                                                                <  1:721067/150‑1 (MQ=255)
      gcgcCTCGGTACTGCAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGtgctgc                                                                                                                                <  1:702680/150‑1 (MQ=255)
          cTCGGTACTGCAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGaa                                                                                                                            >  2:202496/1‑150 (MQ=255)
          cTCGGTACTGCAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGaa                                                                                                                            >  2:588558/1‑150 (MQ=255)
                   ctgatccgcccggtgttctcttccttgttcgaacgggcgttcaccgccagcacgttgatcttctcgttcagcagcccctgcgccaaatCGCGCCCACGAGCGCCCGGGCCGGCTTCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGAt                                                                                                                   >  2:360020‑M1/89‑150 (MQ=255)
                    cAGCTGGCGGGCAAAGGGCCAGAGCGCATGATCCAGGTGCGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATc                                                                                                                  <  1:333438/150‑1 (MQ=255)
                    cAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATc                                                                                                                  >  2:492309/1‑150 (MQ=255)
                      gCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCaa                                                                                                                >  2:453113/1‑150 (MQ=255)
                      atcagcgcggtgttgtcttccttgttcgaacgggtgttcaccgccagcacgttgatcttctcgttcagcagcacctgcgacacatCGCGCCTACGACGGCCCGTTCCTGCTGCGCGATGTGTGGCAGGTGCGGCTGAACGAGAAGATcag                                                                                                                <  1:360020‑M1/65‑2 (MQ=255)
                       cTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAAc                                                                                                               <  1:686764/150‑1 (MQ=255)
                                     gCCAGAGCGCATGATCCAGGTGAGCGGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGaa                                                                                                 >  1:705668/1‑150 (MQ=255)
                                              ttcgaacgggtgttcaccgccagcacgttgatcttctcgttcagcagcacctgcgacacatCGCGCCTACGACCGCCCGGGCCTGCTGCGCGAtgt                                                                                                                                              <  2:526766‑M1/35‑1 (MQ=255)
                                              ttcgaacgggtgttcaccgccagcacgttgatcttctcgttcagcagcacctgcgacacatCGCGCCTACGACCGCCCGGGCCTGCTGCGCGAtgt                                                                                                                                              >  1:526766‑M1/62‑96 (MQ=255)
                                              cATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTc                                                                                        <  2:356406/150‑1 (MQ=255)
                                              cATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTc                                                                                        <  2:347321/150‑1 (MQ=255)
                                               aTGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCg                                                                                       >  1:591272/1‑150 (MQ=255)
                                                 gATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGaa                                                                                     <  2:591272/150‑1 (MQ=255)
                                                      aGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAgg                                                                                >  2:158616/1‑150 (MQ=255)
                                                          gAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGa                                                                            <  1:158616/150‑1 (MQ=255)
                                                           aGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACCCCCGTTCGAACAAGGTAGac                                                                           <  2:647942/150‑1 (MQ=255)
                                                             cTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCGTCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGacaa                                                                         >  1:141591/1‑150 (MQ=255)
                                                                        tCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGacaac                                                                        <  2:752015/140‑1 (MQ=255)
                                                                        tCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGacaac                                                                        >  1:752015/1‑140 (MQ=255)
                                                                         ctgatcttctcgttcagccgcacctgccacacatcgcgcagcaggaacgggccgtcgtaggcgcgATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGAt                                                             <  1:360020‑M2/85‑1 (MQ=255)
                                                                            atcttctcgttcagcagcacctgcgacacatcgcgaagccggcccgggcgctcgtgggcgcgATTTGGCGCAGGGGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACGCCCGTTCGAACAAGGAAGAGAACACCGGGCggatcag                                                          >  2:360020‑M2/63‑147 (MQ=255)
                                                                                cAGACCTACCCAGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCACCGTGCTGTCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTg                                                      >  2:48495/1‑150 (MQ=255)
                                                                                  gACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAAc                                                                                                               <  1:445914/91‑1 (MQ=255)
                                                                                  gACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAAc                                                                                                               >  2:445914/1‑91 (MQ=255)
                                                                                        cccGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGac                                                                           <  2:50763/121‑1 (MQ=255)
                                                                                        cccGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGac                                                                           >  1:50763/1‑121 (MQ=255)
                                                                                          cGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGt                                                                                          <  1:113855/104‑1 (MQ=255)
                                                                                          cGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGt                                                                                          >  2:113855/1‑104 (MQ=255)
                                                                                                 cATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTTCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAgg                                     >  1:789044/1‑150 (MQ=255)
                                                                                                  aTCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGc                                    >  2:640442/1‑150 (MQ=255)
                                                                                                  aTCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAACACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGc                                    >  2:544657/1‑150 (MQ=255)
                                                                                                     gTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACTAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTg                                 <  2:180885/150‑1 (MQ=255)
                                                                                                       acatcgcgcagcaggcccgggcggtcgtaggcgcgATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGaa                                                                                     >  1:526766‑M2/36‑96 (MQ=255)
                                                                                                       acatcgcgcagcaggcccgggcggtcgtaggcgcgATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGaa                                                                                     <  2:526766‑M2/61‑1 (MQ=255)
                                                                                                                cTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTggg                      >  2:435718/1‑150 (MQ=255)
                                                                                                                    gACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGgcgc                  <  1:202496/150‑1 (MQ=255)
                                                                                                                         cccGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCctgct             >  2:173312/1‑150 (MQ=255)
                                                                                                                           cGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGAGCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTgg           <  2:401291/150‑1 (MQ=255)
                                                                                                                                 tgctgcGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGACCTGGACGCGCTGGGGCGCCTGCTGGGGCGGa     >  2:405832/1‑150 (MQ=255)
                                                                                                                                  gctgcGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGAt    <  1:588558/150‑1 (MQ=255)
                                                                                                                                    tgcGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGAtct  >  2:722015/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                               
AGGACTGCGCCTCGGTACTGCAGCTGGCGGGCAAAGAGCCAGAGCGCATGATCCAGGTGAGCTGGGGGCCGATCCCGGTGCAGACCTACCCGGTCGACATCGTCATCCGCGCCTACGACCGCCCGGGCCTGCTGCGCGATGTGTCGCAGGTGCTGCTGAACGAGAAGATCAACGTGCTGGCGGTGAACACCCGTTCGAACAAGGAAGACAACACCGCGCTGATGTCGCTGACCATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCT  >  NC_002947/1852869‑1853150

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