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.

GAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGCAGCCAGCAGA  >  NZ_CP009273/3440134‑3440294
                                                                                      |                                                                          
gAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCt                                                                         <  2:596644/90‑1 (MQ=255)
gAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCt                                                                         <  2:788223/90‑1 (MQ=255)
 agtaaggtagggttactaaatggctaagcaatcaatgaaagcacgcgaagtaaaacgcgtagctttagctgataaatacttcgcgTGCtt                                                                        <  2:791990‑M2/5‑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:608221/1‑90 (MQ=255)
               gagaTGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTa                                                          <  2:137907/90‑1 (MQ=255)
                 gatgatCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGc                                                        <  2:20116/90‑1 (MQ=255)
                    gatCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCAt                                                     >  2:665289/1‑90 (MQ=255)
                         cTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGt                                                <  2:754354/90‑1 (MQ=255)
                         cTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGt                                                >  2:831486/1‑90 (MQ=255)
                          tttCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGc                                                                          >  2:828429/1‑63 (MQ=255)
                          tttCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGc                                                                          <  1:828429/63‑1 (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)
                                               gaagtaaaacgcgtagctttagctgataaatacttcgcgTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAAc                          <  1:871885‑M2/51‑1 (MQ=255)
                                                   taaaacgcgtagctttagctgataaatacttcgcgTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGa                      <  1:493118‑M2/55‑1 (MQ=255)
                                                          cAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCaaa               <  1:806229/90‑1 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTa                                                          >  1:410094/1‑45 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTa                                                          <  2:410094/45‑1 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGt                                                <  1:182262/55‑1 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGt                                                >  2:182262/1‑55 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             >  2:489812/1‑90 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             >  1:884647/1‑90 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             >  1:851242/1‑90 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             <  1:820416/90‑1 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             <  1:791990/90‑1 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             >  1:500419/1‑90 (MQ=255)
                                                            gCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAgg             <  1:373600/90‑1 (MQ=255)
                                                              tAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGca           >  2:335876/1‑90 (MQ=255)
                                                              tAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGca           >  2:678884/1‑90 (MQ=255)
                                                                     aCGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGCAGCcagca    >  2:600508/1‑90 (MQ=255)
                                                                     aCGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGCAGCcagca    >  2:571777/1‑90 (MQ=255)
                                                                       gcgTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGCAGCCAGCaga  >  1:113348/1‑90 (MQ=255)
                                                                                      |                                                                          
GAAGCGTTCACATCAGAGATGATCGCTTTCAGTTCAGCGCGTTTCGCGAAGTATTTATCAGCTAAAGCTACGCGTTTTACTTCGCGTGCTTTCATTGATTGCTTAGCCATTTAGTAACCCTACCTTACTTGCGGAACGGGAAGTCAAAGGCAGCCAGCAGA  >  NZ_CP009273/3440134‑3440294

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