New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4306900 =51 (0.960)3 (0.060) 3/150 4.0 5.4% coding (136/1137 nt) phnM alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
?NZ_CP009273 4306937 = 60 (1.250)coding (99/1137 nt) phnM alpha‑D‑ribose 1‑methylphosphonate 5‑triphosphate diphosphatase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGC  >  NZ_CP009273/4306812‑4306989
                                                                                        |                                                                                         
gAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGcagccag                                                                                          >  1:275347/1‑90 (MQ=255)
     ggCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccgcc                                                                                     >  2:1355481/1‑90 (MQ=255)
        aGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCtta                                                                                  <  1:917151/90‑2 (MQ=255)
          gCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGc                                                                                <  1:142956/90‑1 (MQ=255)
                 aaCTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt                                                                         <  2:404285/90‑1 (MQ=255)
                 aaCTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt                                                                         <  1:1610317/90‑1 (MQ=255)
                   cTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccgc                                                                                      >  2:1043043/1‑75 (MQ=255)
                   cTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAgccgc                                                                                      <  1:1043043/75‑1 (MQ=255)
                              ggggTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAgg                                                            >  1:227127/1‑90 (MQ=255)
                               gggTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAggc                                                           <  1:776617/90‑1 (MQ=255)
                               gggTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGAAggc                                                           <  1:1361789/90‑1 (MQ=255)
                                     agaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGctct                                                     <  1:1168051/90‑1 (MQ=255)
                                         tttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCt                                                 >  2:337984/1‑90 (MQ=255)
                                         tttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCt                                                 >  1:456799/1‑90 (MQ=255)
                                         tttATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGACTCCGGCAGGCGGCTCTGGCt                                                 >  2:968047/1‑90 (MQ=255)
                                           tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCttt                                               >  2:318058/1‑90 (MQ=255)
                                           tATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCttt                                               >  2:378993/1‑90 (MQ=255)
                                             tCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGcc                                                                               <  1:1526862/56‑1 (MQ=255)
                                             tCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGcc                                                                               >  2:1526862/1‑56 (MQ=255)
                                             tCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCg                                             >  1:425250/1‑90 (MQ=255)
                                                aGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCa                                          >  2:703457/1‑90 (MQ=255)
                                                aGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCa                                          >  1:1132049/1‑90 (MQ=255)
                                                aGATTATCGGTATGCAACTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCa                                          >  1:1237346/1‑90 (MQ=255)
                                                  aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCaaa                                        <  2:275347/90‑1 (MQ=255)
                                                  aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCaaa                                        >  2:1144478/1‑90 (MQ=255)
                                                  aTTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAAGCGGCTCTGGCTTTCGGCaaa                                        >  2:71349/1‑90 (MQ=255)
                                                       cGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGgcgc                                   <  2:820632/90‑1 (MQ=255)
                                                           aTGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                               >  2:11687/1‑90 (MQ=255)
                                                           aTGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                               >  1:1378470/1‑90 (MQ=255)
                                                           aTGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                               >  2:1525253/1‑90 (MQ=255)
                                                             gCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAttt                             >  1:1118218/1‑90 (MQ=255)
                                                               aGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCg                           >  2:387688/1‑90 (MQ=255)
                                                               aGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCg                           >  2:1585299/1‑90 (MQ=255)
                                                               aGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCg                           >  2:245792/1‑90 (MQ=255)
                                                                gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt                                                                         <  1:1216836/43‑1 (MQ=255)
                                                                gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCAt                                                                         >  2:1216836/1‑43 (MQ=255)
                                                                gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGc                          >  2:552290/1‑90 (MQ=255)
                                                                gCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGc                          >  1:11801/1‑90 (MQ=255)
                                                                   cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTc                            >  2:1402638/1‑85 (MQ=255)
                                                                   cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTc                            <  1:1402638/85‑1 (MQ=255)
                                                                   cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                       <  1:318058/90‑1 (MQ=255)
                                                                   cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                       >  1:1170239/1‑90 (MQ=255)
                                                                   cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                       >  2:850382/1‑90 (MQ=255)
                                                                   cAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGACTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGt                       <  1:968047/90‑1 (MQ=255)
                                                                          ccccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAc                >  2:1279926/1‑90 (MQ=255)
                                                                           cccGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAcc               <  2:1502431/90‑1 (MQ=255)
                                                                             cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAggc                                                           >  1:1179032/1‑44 (MQ=255)
                                                                             cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAggc                                                           <  2:1179032/44‑1 (MQ=255)
                                                                             cGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTc             >  2:1071202/1‑90 (MQ=255)
                                                                                   cagccagcCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGa       <  1:1279926/90‑1 (MQ=255)
                                                                                      ccagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAAcc    <  2:227127/90‑1 (MQ=255)
                                                                                       gagccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCg   >  2:637559/2‑90 (MQ=255)
                                                                                        agccGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGc  >  2:529735/1‑90 (MQ=255)
                                                                                        |                                                                                         
GAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGC  >  NZ_CP009273/4306812‑4306989

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.