New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 215975 =42 (1.190)3 (0.100) 3/146 2.7 8.6% coding (107/708 nt) tsaA tRNA (N6‑threonylcarbamoyladenosine(37)‑N6)‑methyltransferase TrmO
?NZ_CP009273 216027 = 27 (0.880)coding (55/708 nt) tsaA tRNA (N6‑threonylcarbamoyladenosine(37)‑N6)‑methyltransferase TrmO
Rejected: Frequency below/above cutoff threshold.

AACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATACATAC  >  NZ_CP009273/215938‑216116
                                                                                         |                                                                                         
aaCGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtt                                                                                           <  2:647944/90‑1 (MQ=255)
aaCGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtt                                                                                           <  1:7830/90‑1 (MQ=255)
  cGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtttt                                                                                         <  1:104933/90‑1 (MQ=255)
                 tGTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGgcgagc                                                                          <  2:235547/90‑1 (MQ=255)
                  gTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGgcgagcg                                                                         <  1:411277/90‑1 (MQ=255)
                   tAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGgcgagcga                                                                        <  1:936813/90‑1 (MQ=255)
                        aGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtttt                                                                                         >  1:260466/1‑68 (MQ=255)
                        aGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtttt                                                                                         <  2:260466/68‑1 (MQ=255)
                        aGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAAc                                                                   >  1:225048/1‑90 (MQ=255)
                        aGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAAc                                                                   >  2:179783/1‑90 (MQ=255)
                        aGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAAc                                                                   >  1:892410/1‑90 (MQ=255)
                         gCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACg                                                                  >  1:164984/1‑90 (MQ=255)
                         gCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACg                                                                  >  2:716232/1‑90 (MQ=255)
                               agaTGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTAtt                                                            >  1:923549/1‑90 (MQ=255)
                               agaTGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTAtt                                                            >  1:643220/1‑90 (MQ=255)
                                   gCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCt                                                        <  1:83662/90‑1 (MQ=255)
                                   gCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCt                                                        <  1:626199/90‑1 (MQ=255)
                                     AGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCg                                                      >  1:1022328‑M1/1‑90 (MQ=255)
                                      GTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGa                                                     <  2:988509‑M1/90‑1 (MQ=255)
                                        TCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAAc                                                   <  2:30029‑M1/90‑1 (MQ=255)
                                          GCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTg                                                 >  2:923549‑M1/1‑90 (MQ=255)
                                          GCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTg                                                 >  2:643220‑M1/1‑90 (MQ=255)
                                                ggCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtt                                                                                           >  1:207538/1‑42 (MQ=255)
                                                ggCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACtt                                                                                           <  2:207538/42‑1 (MQ=255)
                                                         tACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTc                                  <  2:171533/90‑1 (MQ=255)
                                                           ccAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGc                                <  1:775606/90‑1 (MQ=255)
                                                           ccAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGc                                <  1:389672/90‑1 (MQ=255)
                                                             aGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCg                              <  2:626417/90‑1 (MQ=255)
                                                               ctggtaaaaagcgccaacggcgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGa                                               <  2:394092‑M2/45‑1 (MQ=255)
                                                               ctggtaaaaagcgccaacggcgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGa                                               >  1:394092‑M2/27‑71 (MQ=255)
                                                                 ggtaaaaagcgccaacggcgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAAt                          >  1:259504‑M2/25‑90 (MQ=255)
                                                                          gcgGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAgcgc                 >  2:1093577/1‑90 (MQ=255)
                                                                          gcgGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAgcgc                 <  1:981328/90‑1 (MQ=255)
                                                                             gAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGa              <  2:133839/90‑1 (MQ=255)
                                                                             gAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAAAGTTAAGCGCAGa              <  2:164984/90‑1 (MQ=255)
                                                                               acggcgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAAc            <  2:1047186‑M2/80‑1 (MQ=255)
                                                                               acggcgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAAc            <  2:184720‑M2/80‑1 (MQ=255)
                                                                                  gcgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGa         <  1:354699‑M2/83‑1 (MQ=255)
                                                                                   cgaactTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGat        <  1:496540‑M2/84‑1 (MQ=255)
                                                                                       ctTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGatacat    >  1:703960‑M2/3‑90 (MQ=255)
                                                                                        tTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGatacata   <  1:360973‑M2/89‑1 (MQ=255)
                                                                                         tttCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGatacatac  <  1:290027/90‑1 (MQ=255)
                                                                                         |                                                                                         
AACGGCGTCGGCCTGGTTGTAGGGAGCAATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATACATAC  >  NZ_CP009273/215938‑216116

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.