New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2622492 =87 (0.550)3 (0.020) 3/126 14.6 3.9% coding (1429/2067 nt) ppk polyphosphate kinase
?NC_000913 2622504 = 75 (0.560)coding (1441/2067 nt) ppk polyphosphate kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAAT  >  NC_000913/2622430‑2622574
                                                                            |                                                                      
ggACGGGGAACTTAAACAAAAAACCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCAc                                                                       <  2:4436048/76‑1 (MQ=255)
 ggccGGGGATTTTAAAAAAAAACCCCGCGCGTTTTTATACTAACTATTCGTTGCTGACCGCCGATGCGCGCATCAc                                                                       <  1:2322439/74‑1 (MQ=255)
 gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGTCCGCCCCACCGCGCATcccc                                                                      >  1:1570775/1‑76 (MQ=255)
 gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCAcc                                                                      >  1:654267/1‑76 (MQ=255)
 gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCCATGCGCGCCTCAcc                                                                      >  2:3218302/1‑76 (MQ=255)
  aCCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCa                                                                     >  1:3849562/1‑76 (MQ=255)
   ccGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCac                                                                    >  2:2653607/1‑75 (MQ=255)
   ccGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCaa                                                                    >  1:3188187/1‑76 (MQ=255)
    cGGGACTTTTAGCGAAAACCCCGTGTGTTTTAATCCAGATTCTGCGGTGCGGACCGCCGATGCGCGCATCACCAAc                                                                   <  1:2266405/76‑1 (MQ=255)
    cGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGCTGCGCGCATCACCCAc                                                                   >  2:507628/1‑76 (MQ=255)
    cGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAAc                                                                   >  1:3297613/1‑76 (MQ=255)
     gggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACg                                                                  >  1:5637124/1‑76 (MQ=255)
     gggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACg                                                                  >  1:1859721/1‑76 (MQ=255)
      ggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGa                                                                 >  2:2208869/1‑76 (MQ=255)
        aagtttAACGAAAAACCCGCGCGTCTTTATACTGACTTTTCGTTGCTGCCCGCCGATGCGCGCATCACCAACGAAg                                                               <  2:1980187/73‑1 (MQ=255)
           tttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTAc                                                            >  1:1162383/1‑76 (MQ=255)
           tttAAAAAAAAACCCGCGCGTTTTTATACGGACTATTCGTTGCTGCCCGCCGATGCGCGCATCACCAACGAAGTAc                                                            <  2:2509046/76‑1 (MQ=255)
            ttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTAcg                                                           >  2:4412264/1‑76 (MQ=255)
              aaCGAAAAAACCGCGCGTCTTTATACTGCCTATTCGTTGCTGACCGCCGCTGCCCCCCTCCCCACCGCAGTCcgcg                                                         >  2:160585/1‑74 (MQ=255)
              aaCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTAcggc                                                         >  1:85585/1‑76 (MQ=255)
              aaCAAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTAcggc                                                         <  2:717099/76‑1 (MQ=255)
                caaaaaacCCGCGCGTTTTTATACTGACTATTCGTTGCGGACCGCCGATGCGCGCATCACCAACGAAGTAcggcgg                                                       <  2:4690025/74‑1 (MQ=255)
                cGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCCTCACCACAGAAGCAcggcgg                                                       >  2:1561309/1‑76 (MQ=255)
                aaaaaaacCCGCGCGTCTTTATACGGACTTTGCGTTGCTGCCCGCCGATGCGCGCATCACCAACGAAGTAcggcgg                                                       <  2:951700/74‑1 (MQ=255)
                 aaaaaacCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCggg                                                      <  1:3644588/75‑1 (MQ=255)
                  aaaaaaCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCggga                                                     >  1:151936/1‑75 (MQ=255)
                  aaaaaaCCGAGCGTCTTTATACTGACTATTCGTGGCTGAACGCCGATGGGCGCATCACCAACGAAACAcggcggca                                                     >  1:203250/1‑74 (MQ=255)
                    aaacCCGCGCGTGTTTTTACGGATTTTTCGTTGCTGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTAt                                                   <  2:3269324/76‑1 (MQ=255)
                    aaaaCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGAt                                                   >  1:2875673/1‑76 (MQ=255)
                        aCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTa                                                >  2:4256284/1‑76 (MQ=255)
                         ccGCGCTTTTTTATACTGATTTTTCGTTGCGGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTaa                                               <  2:3147093/76‑1 (MQ=255)
                          cgcgctTTTTTATCCTGTTTATTCGTTGCTGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAc                                              <  1:4317917/76‑1 (MQ=255)
                          cgcgcgTTTTTATACTGACTATTGTTGTCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAAc                                              <  2:3554582/76‑1 (MQ=255)
                           gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACt                                             >  2:3704474/1‑76 (MQ=255)
                           gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACt                                             >  1:5649306/1‑76 (MQ=255)
                           gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACt                                             >  2:3012401/1‑76 (MQ=255)
                             gcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACttt                                           >  1:5234307/1‑76 (MQ=255)
                              cgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTa                                          >  1:234035/1‑76 (MQ=255)
                              cgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAAttaa                                          >  1:4165974/1‑74 (MQ=255)
                                tttttATGCTGTCTATTCGTTGCTCCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTAtt                                        <  2:1570775/74‑1 (MQ=255)
                                tCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACtttttt                                        >  2:3709175/1‑76 (MQ=255)
                                 cTTTATACTGACTATTCGATGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTAATTAACTTTATTg                                       >  1:1120580/1‑76 (MQ=255)
                                  tttATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTGa                                      >  1:4786490/1‑76 (MQ=255)
                                  tttATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTGa                                      >  1:393694/1‑76 (MQ=255)
                                  tttATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTGa                                      >  1:5201716/1‑76 (MQ=255)
                                   ttATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTAtttta                                     >  1:608953/1‑73 (MQ=255)
                                    tatgCTGCCTTTGCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGaaa                                    <  2:4954566/76‑1 (MQ=255)
                                     ttaCAGATTATTCTTTCCGGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGaaaa                                   <  1:3811333/75‑1 (MQ=255)
                                      taCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTTTAAAc                                  >  1:2710389/1‑76 (MQ=255)
                                       aCTGACTATGCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAAcc                                 <  2:2270110/76‑1 (MQ=255)
                                        cTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTTAAAAccc                                >  2:5580378/1‑76 (MQ=255)
                                          gACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTTAAAAcccct                              >  2:5237155/1‑74 (MQ=255)
                                          accTATTCGTTGCGGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCAt                              <  1:4318222/74‑1 (MQ=255)
                                             tATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATAcc                           >  1:5021858/1‑76 (MQ=255)
                                              aTTCGTTGCTGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCg                          <  2:1140621/76‑1 (MQ=255)
                                              aTTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTTTTTAAAACCCATAccc                          >  2:1635658/1‑75 (MQ=255)
                                              aTTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTTTTGAAAACCCCTAccc                          >  2:4616542/1‑75 (MQ=255)
                                              aTTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTGAAAACCCATAccc                          >  2:1678991/1‑75 (MQ=255)
                                               ttCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCgc                         >  2:3000238/1‑75 (MQ=255)
                                               ttCGTGGCTGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGt                         <  2:3668626/76‑1 (MQ=255)
                                                tCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTTTTTAAAACCCCTACCGTc                        >  2:5258581/1‑76 (MQ=255)
                                                tCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAAATCTCTTGAAAACCCATACCGTc                        >  2:3980509/1‑76 (MQ=255)
                                                  gtggCTGCCCGCCTATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAATCCCATCCCGTCCg                      <  2:1797394/73‑1 (MQ=255)
                                                  gTTGCTCCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCg                      <  1:4575847/76‑1 (MQ=255)
                                                    ggCTGCCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGt                    <  2:3749358/75‑1 (MQ=255)
                                                     gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTg                   >  1:3026342/1‑76 (MQ=255)
                                                          aCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTGAAAACCCATACCGTCCGGTGACAt               >  1:3662962/1‑76 (MQ=255)
                                                          aCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTGAAAACCCACACCCTCCGGTGCCac               >  2:468705/1‑75 (MQ=255)
                                                          aCCGACGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTGTTGTGAAAACACATAACGTCCTGTGTCAt               >  1:5191263/1‑76 (MQ=255)
                                                           ccgccgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACAtt              <  2:2833620/76‑1 (MQ=255)
                                                            cgccgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACAttt             <  1:2334829/76‑1 (MQ=255)
                                                             gccgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTCACTTTATTGAAAACCCCTACCGTCCGGGGACATTTg            >  2:4873498/1‑76 (MQ=255)
                                                             gccgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCCTACCGTCCGGTGACCTTTg            >  2:4840467/1‑76 (MQ=255)
                                                             gccgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTg            >  1:5193030/1‑76 (MQ=255)
                                                              ccgATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGa           <  2:3662962/76‑1 (MQ=255)
                                                                       gcATCACCAAGGTCGTGCGGCGGGTATTTAACTTTATGGAACACCCATACCGTCCGGTGACATTTGATTCTTTact  <  1:5527846/76‑3 (MQ=255)
                                                                            |                                                                      
GGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAAT  >  NC_000913/2622430‑2622574

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