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.

GTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACC  >  NZ_CP009273/4306814‑4306987
                                                                                         |                                                                                    
agcctcttcatgattatcaataacgttaagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgc                                                                                      >  1:758981‑M1/90‑90 (MQ=255)
 gcctcttcatgattatcaataacgttaagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgcg                                                                                     <  2:281269‑M1/2‑1 (MQ=255)
   ggCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccgcc                                                                                   >  1:438744/1‑90 (MQ=255)
   ctcttcatgattatcaataacgttaagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgCGcc                                                                                   <  2:732016‑M1/4‑1 (MQ=255)
   ctcttcatgattatcaataacgttaagctggtgctggaaaacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgCGcc                                                                                   <  1:90873‑M1/4‑1 (MQ=255)
                    tCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCt                                                                  >  2:170999/1‑90 (MQ=255)
                      gggCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTcc                                                                <  2:828484/90‑1 (MQ=255)
                           cgGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAg                                                           <  1:16098/90‑1 (MQ=255)
                                         tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCttt                                             >  2:703203/1‑90 (MQ=255)
                                         tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCttt                                             >  2:152438/1‑90 (MQ=255)
                                          aacgaggtggtaagcggttcgctggaggtgcagaacggcgaaatccgCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTc                                            >  2:90873‑M1/48‑90 (MQ=255)
                                           tCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCACCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCg                                           >  2:534845/1‑90 (MQ=255)
                                                         aTGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                             >  1:260585/1‑90 (MQ=255)
                                                           gCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAttt                           >  2:777470/1‑90 (MQ=255)
                                                                 cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                     >  1:398988/1‑90 (MQ=255)
                                                                 cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                     >  2:834022/1‑90 (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                             <  1:832042/66‑1 (MQ=255)
                                                                                 cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                             >  2:832042/1‑66 (MQ=255)
                                                                                   gccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAc   >  2:345890/1‑90 (MQ=255)
                                                                                   gccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAc   <  2:764964/90‑1 (MQ=255)
                                                                                    atccgCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAcc  <  2:102126‑M1/85‑1 (MQ=255)
                                                                                         |                                                                                    
GTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACC  >  NZ_CP009273/4306814‑4306987

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