New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4306903 =20 (0.690)3 (0.110) 3/158 2.1 9.9% coding (133/1137 nt) phnM alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
?NZ_CP009273 4306956 = 36 (1.330)coding (80/1137 nt) phnM alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Rejected: Frequency below/above cutoff threshold.

ATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTC  >  NZ_CP009273/4306871‑4307045
                                                                                     |                                                                                         
ggttcgctggaggtgcagaacggcgaaatccgcgcctttgccgaaagccagagccgcctgccggaggcgatggacggcgaaggcgCGGAt                                                                                       <  2:102126‑M2/5‑1 (MQ=255)
aTGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                                                                                       >  1:260585/1‑90 (MQ=255)
  gCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAttt                                                                                     >  2:777470/1‑90 (MQ=255)
        cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                                                                               >  2:834022/1‑90 (MQ=255)
        cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                                                                               >  1:398988/1‑90 (MQ=255)
             tgcagaacggcgaaatccgcgcctttgccgaaagccagagccgcctgccggaggcgatggacggcgaaggcgCGGATTTCGCCGTTCTGc                                                                          <  2:183504‑M2/18‑1 (MQ=255)
              gcagaacggcgaaatccgCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCa                                                                         <  2:183504‑M1/72‑1 (MQ=255)
              gCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCa                                                                         >  2:674511/1‑90 (MQ=255)
               ccccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAc                                                                        >  2:433745/1‑90 (MQ=255)
                cccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAcc                                                                       <  2:719999/90‑1 (MQ=255)
                cccGGCAGCAGCCAGCCGCCTTCGCCGTCCACCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAcc                                                                       <  1:39/90‑1 (MQ=255)
                  cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTc                                                                     >  1:724853/1‑90 (MQ=255)
                        cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                                                                                       >  2:832042/1‑66 (MQ=255)
                        cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                                                                                       <  1:832042/66‑1 (MQ=255)
                          gccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAc                                                             >  2:345890/1‑90 (MQ=255)
                          gccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAc                                                             <  2:764964/90‑1 (MQ=255)
                           atccgCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAcc                                                            <  2:102126‑M1/85‑1 (MQ=255)
                                      cGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCt                                                 >  1:690190/1‑90 (MQ=255)
                                          gaaagccagagccgcctgccggaggcgatggacggcgaaggcgCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGtt                                             >  2:90873‑M2/44‑90 (MQ=255)
                                                   cTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAc                                    <  2:71216/90‑1 (MQ=255)
                                                   cTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAc                                    <  2:260585/90‑1 (MQ=255)
                                                   cTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAc                                    <  2:377740/90‑1 (MQ=255)
                                                    tCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAcc                                   <  1:20255/90‑1 (MQ=255)
                                                      cGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAcc                                   >  1:778306/1‑88 (MQ=255)
                                                      cGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCAcc                                   <  2:778306/88‑1 (MQ=255)
                                                         cAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCTCCAGCtt                              <  1:163515/90‑1 (MQ=255)
                                                         cAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt                              <  1:345890/90‑1 (MQ=255)
                                                            gcggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt                              >  2:271220/1‑87 (MQ=255)
                                                            gcggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt                              >  2:286152/1‑87 (MQ=255)
                                                            gcggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt                              <  1:286152/87‑1 (MQ=255)
                                                            gcggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt                              <  1:271220/87‑1 (MQ=255)
                                                             cggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg                          <  1:336370/90‑1 (MQ=255)
                                                              ggcTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGt                         <  2:333104/90‑1 (MQ=255)
                                                                ctctGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTa                       <  1:823315/90‑1 (MQ=255)
                                                                ctctGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTa                       >  2:817199/1‑90 (MQ=255)
                                                                   tGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTg                    <  1:196818/90‑1 (MQ=255)
                                                                     gCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGAt                  >  2:220929/1‑90 (MQ=255)
                                                                     gCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGAt                  <  2:690190/90‑1 (MQ=255)
                                                                       tttCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATaa                <  1:728835/90‑1 (MQ=255)
                                                                             cAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGa          <  2:438744/90‑1 (MQ=255)
                                                                              aaaGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGaa         >  2:723084/1‑90 (MQ=255)
                                                                                 ggcgCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAagag      <  2:732016‑M2/86‑1 (MQ=255)
                                                                                 ggcgCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAagag      <  1:90873‑M2/86‑1 (MQ=255)
                                                                                   cgCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGc    <  2:281269‑M2/88‑1 (MQ=255)
                                                                                    gCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCt   >  1:758981‑M2/2‑90 (MQ=255)
                                                                                     cgGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTc  <  1:170999/90‑1 (MQ=255)
                                                                                     |                                                                                         
ATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACGTTATTGATAATCATGAAGAGGCTC  >  NZ_CP009273/4306871‑4307045

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