New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4623274 =26 (0.860)4 (0.150) 4/150 1.7 14.5% coding (289/513 nt) yjjX inosine/xanthosine triphosphatase
?NZ_CP009273 4623410 = 24 (0.890)coding (153/513 nt) yjjX inosine/xanthosine triphosphatase
Rejected: Frequency below/above cutoff threshold.

CATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATC  >  NZ_CP009273/4623322‑4623494
                                                                                        |                                                                                    
cATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGcgcgcg                                                                                     >  1:156045/1‑90 (MQ=255)
     aTGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTcc                                                                                >  1:721396/1‑90 (MQ=255)
          aGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACt                                                                           <  1:98907/90‑1 (MQ=255)
           gCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCAc                                                                            <  1:298375/88‑1 (MQ=255)
           gCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCAc                                                                            >  2:298375/1‑88 (MQ=255)
               cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCaa                                                                      <  1:345753/90‑1 (MQ=255)
                aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCaaa                                                                     <  2:835210/90‑1 (MQ=255)
                    gCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGc                                                                 <  1:194335/90‑1 (MQ=255)
                     ccACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCt                                                                <  1:287429/90‑1 (MQ=255)
                     attgaagctggcatcgatggtgacagcactttcagttgggtggtgattgaaaacgccagccagcgcgCGTTTCCTCACTGCCAAAGGGCt                                                                <  2:611528‑M2/23‑1 (MQ=255)
                       tgaagctggcatcgatggtgacagcactttcagttgggtggtgattgaaaacgccagccagcgcgCGTTTCCTCACTGCCAAAGGGCTGt                                                              <  1:434541‑M2/25‑1 (MQ=255)
                       aCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGt                                                              <  1:162927/90‑1 (MQ=255)
                       aCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGt                                                              <  2:163031/90‑1 (MQ=255)
                       aCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGt                                                              <  1:596630/90‑1 (MQ=255)
                       aCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGt                                                              <  1:276500/90‑1 (MQ=255)
                               aTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTac                                                      >  1:614874/1‑90 (MQ=255)
                                 atcgatggtgacaccatttccagttgggtggtgattgaaaacgccagccagcgcgCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTacac                                                    <  1:92539‑M2/35‑1 (MQ=255)
                                  aGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACAcc                                                   >  2:693916/1‑90 (MQ=255)
                                  aGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACAcc                                                   >  1:448302/1‑90 (MQ=255)
                                                acttgcagttgggtggtgattgaaaacgccagccagcgcgCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCaa                                     <  1:52536‑M2/50‑1 (MQ=255)
                                                    gcggGCATTGGCTACCTGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGa                                 >  1:432379/1‑90 (MQ=255)
                                                     cggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGAt                                <  2:231580/90‑1 (MQ=255)
                                                       ggCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGc                              <  1:432293/90‑1 (MQ=255)
                                                             ggCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCaa                                                                      <  2:691489/44‑1 (MQ=255)
                                                             ggCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCaa                                                                      >  1:691489/1‑44 (MQ=255)
                                                                       tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCt              >  1:58925/1‑90 (MQ=255)
                                                                       tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCt              >  2:543246/1‑90 (MQ=255)
                                                                       tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCt              >  2:743808/1‑90 (MQ=255)
                                                                               gccagcgcgCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGaa      >  1:289703‑M2/10‑90 (MQ=255)
                                                                                 cagcgcgCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGa    <  2:448302‑M2/83‑1 (MQ=255)
                                                                                   gcgcgCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGAtc  <  1:611528‑M2/85‑1 (MQ=255)
                                                                                        |                                                                                    
CATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATC  >  NZ_CP009273/4623322‑4623494

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