New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 548215 =23 (1.160)4 (0.250) 4/132 1.1 18.8% coding (342/510 nt) purE 5‑(carboxyamino)imidazole ribonucleotide mutase
?NZ_CP009273 548231 = 16 (1.010)coding (326/510 nt) purE 5‑(carboxyamino)imidazole ribonucleotide mutase
Rejected: Frequency below/above cutoff threshold.

AATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAGCACCGGCACCAGCGTTTT  >  NZ_CP009273/548149‑548310
                                                                                  |                                                                               
aaTTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGAt                                                                          <  2:18283/90‑1 (MQ=255)
  tttGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATgg                                                                        <  2:140444/90‑1 (MQ=255)
        ctgcCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTaga                                                                  <  1:59087/90‑1 (MQ=255)
        ctgcCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGAATTTGTACGATGGAGTaga                                                                  <  1:212728/90‑1 (MQ=255)
           tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaaatgccgcgcggcattTGTACGAt                                                                          <  1:307192‑M2/8‑1 (MQ=255)
           tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaaatgccgcgcggcattTGTACGAt                                                                          >  2:307192‑M2/72‑79 (MQ=255)
           tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaaatgccgcgcggcattTGTACGATGGAGTAGAgg                                                                <  2:376131‑M2/18‑1 (MQ=255)
           tgccagtacagagcgccgcactgagcggtgtcgatagcctctactccatcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAgg                                                                <  1:376131‑M1/34‑1 (MQ=255)
           ccAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGc                                                               <  2:2896/90‑1 (MQ=255)
              gTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                                           <  1:553458/90‑1 (MQ=255)
               tAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                                           <  2:595806/90‑1 (MQ=255)
                    cctcgTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGAt                                                                          <  1:329645/67‑1 (MQ=255)
                    cctcgTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGAt                                                                          >  2:329645/4‑70 (MQ=255)
                                     gCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGctc                                     <  2:514313/90‑1 (MQ=255)
                                               tCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGa                                                                           <  1:318792/42‑1 (MQ=255)
                                               tCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGa                                                                           >  2:318792/1‑42 (MQ=255)
                                               tCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           <  1:370245/90‑1 (MQ=255)
                                                cctctactccatcgtacaaatgccgcgcggcattTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           <  1:376131‑M2/55‑1 (MQ=255)
                                                cctctactccatcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           <  2:376131‑M1/71‑1 (MQ=255)
                                                      cGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAg                    <  2:538117/90‑1 (MQ=255)
                                                          atcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           <  1:307192‑M1/71‑1 (MQ=255)
                                                          atcgtacaAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           >  2:307192‑M1/9‑79 (MQ=255)
                                                            cACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                                           <  1:143638/45‑1 (MQ=255)
                                                            cACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATc                                                           >  2:143638/1‑45 (MQ=255)
                                                                   ATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGctc                                     <  2:580635‑M1/60‑1 (MQ=255)
                                                                   ATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGctc                                     >  1:580635‑M1/1‑60 (MQ=255)
                                                                        gcgcggcattTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           >  2:42998‑M2/11‑65 (MQ=255)
                                                                        gcgcggcattTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCa                           <  1:42998‑M2/55‑1 (MQ=255)
                                                                        GCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAGCACCGGCACCAGCGtttt  <  2:23972‑M1/90‑1 (MQ=255)
                                                                                  |                                                                               
AATTTGTGCTGCCAGTAACGCCGCGTTTGCCGCGCCAGCTTTACCAATCGCCAGCGTACCCACCGGAATGCCGCGCGGCATTTGTACGATGGAGTAGAGGCTATCGACACCGCTCAGTGCGGCGCTCTGTACTGGCACGCCCAGCACCGGCACCAGCGTTTT  >  NZ_CP009273/548149‑548310

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