New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3440178 =30 (0.940)3 (0.110) 3/150 2.5 8.6% coding (66/306 nt) rpsN 30S ribosomal protein S14
?NZ_CP009273 3440220 = 37 (1.300)coding (24/306 nt) rpsN 30S ribosomal protein S14
Rejected: Frequency below/above cutoff threshold.

GAGTCTGCAGCTTGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTG  >  NZ_CP009273/3440093‑3440264
                                                                                     |                                                                                      
gagTCTGCAGCTTGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGa                                                                                    <  2:583485/90‑1 (MQ=255)
         gCTTGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTAt                                                                           <  1:38233/90‑1 (MQ=255)
            tGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGa                                                                                    <  2:400882/78‑1 (MQ=255)
            tGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGa                                                                                    >  1:400882/1‑78 (MQ=255)
            tGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAg                                                                        >  2:316902/1‑90 (MQ=255)
                   aGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGa                                                                                    >  1:480808/1‑71 (MQ=255)
                   aGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGa                                                                                    <  2:480808/71‑1 (MQ=255)
                     cGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTa                                                               <  2:204084/90‑1 (MQ=255)
                              tcccgttccgcaagtaaggtagggttactaaatggctaagcaatcaatgaaagcaCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTa                                                      <  1:493118‑M1/35‑1 (MQ=255)
                                  gttccgcaagtaaggtagggttactaaatggctaagcaatcaatgaaagcaCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTc                                                  <  1:871885‑M1/39‑1 (MQ=255)
                                         gAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCt                                           <  2:788223/90‑1 (MQ=255)
                                         gAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCt                                           <  2:596644/90‑1 (MQ=255)
                                                tCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCattg                                    <  2:148583/90‑1 (MQ=255)
                                                  acaTCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCattgat                                  <  1:51589/90‑1 (MQ=255)
                                                   caTCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCattgatt                                 <  2:262227/90‑1 (MQ=255)
                                                       agagATGATCGCTTTCAGTTCAGCGCGTTTCTCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCtt                             <  2:782913/90‑1 (MQ=255)
                                                       agagATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCtt                             <  2:429764/90‑1 (MQ=255)
                                                        gagaTGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTa                            <  2:137907/90‑1 (MQ=255)
                                                        gagaTGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTa                            >  2:608221/1‑90 (MQ=255)
                                                          gatgatCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGc                          <  2:20116/90‑1 (MQ=255)
                                                             gatCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTAt                                                                           >  1:522904/1‑38 (MQ=255)
                                                             gatCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCAt                       >  2:665289/1‑90 (MQ=255)
                                                             gatCGCTTTCAGTTCAGCGCGCTTCGCGAAGTATTTAt                                                                           <  2:522904/38‑1 (MQ=255)
                                                                  cTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGt                  >  2:831486/1‑90 (MQ=255)
                                                                  cTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGt                  <  2:754354/90‑1 (MQ=255)
                                                                   tttCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGc                                            <  1:828429/63‑1 (MQ=255)
                                                                   tttCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGc                                            >  2:828429/1‑63 (MQ=255)
                                                                   tttCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTa                 <  2:874992/90‑1 (MQ=255)
                                                                        gttcagCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCt            >  2:655276/1‑90 (MQ=255)
                                                                            agCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCt        >  1:31728/1‑90 (MQ=255)
                                                                            agCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCt        >  2:342608/1‑90 (MQ=255)
                                                                                aagcaCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACt    <  2:791990‑M1/85‑1 (MQ=255)
                                                                                  tttCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTg  >  2:579191/1‑90 (MQ=255)
                                                                                     |                                                                                      
GAGTCTGCAGCTTGAGAACAGCGTTCCAACGATCTTCGTCGGAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTG  >  NZ_CP009273/3440093‑3440264

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