New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 382440426 (0.890)3 (0.120) 3/148 2.0 9.7% coding (588/1710 nt) yicH AsmA family protein
?NZ_CP009273 = 3824452 33 (1.300)coding (636/1710 nt) yicH AsmA family protein
Rejected: Frequency below/above cutoff threshold.

GCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3824308‑3824404
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgttacgctgCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAATCAGTA  <  NZ_CP009273/3824452‑3824364
                                                                                                                                                                                          
GCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGT                                                                                                  <  2:94501/90‑1
GCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGT                                                                                                  >  1:767168/1‑90
GCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGT                                                                                                  >  1:556766/1‑90
    ATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACG                                                                                              <  1:327005/90‑1
    ATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACG                                                                                              <  2:51822/90‑1
      TCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCT                                                                                            <  1:621684/90‑1
             AGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGCGCGTT                                                                                     >  1:564497/1‑90
                                   ATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGC                                                                                             <  2:395413/60‑1
                                   ATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGC                                                                                             >  1:395413/1‑60
                                             TGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCAC                                                     >  1:292515/1‑90
                                              GCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACC                                                    <  2:564497/90‑1
                                                                                       CGTTACGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTT           <  2:216097/90‑1
                                                                                       CGTTACGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTT           <  2:611961/90‑1
                                                                                          TACGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAA        >  1:807033/1‑90
                                                                                          TACGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAA        <  2:753704/90‑1
                                                                                            CGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAATC      <  1:469620/90‑1
                                                                                            CGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAATC      <  1:804499/90‑1
                                                                                            CGCTGCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAATC      <  2:422923/90‑1
                                                                                                GCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAATCAGTA  >  1:542145/1‑90
                                                                                                GCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGG                                                        <  1:833988/36‑1
                                                                                                GCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGG                                                        >  2:833988/1‑36
                                                                                                                                                                                          
GCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3824308‑3824404
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgttacgctgCGCGTTTCCGGTTAATGTCCCGCGGGCGATGTCGGCACCCAGGTTAGTCAGCGTAACGCGATCGTTATCAATACTGCCTTCAATCAGTA  <  NZ_CP009273/3824452‑3824364

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