New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 120326138 (0.610)6 (0.120) 6/130 2.5 19.5% coding (305/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 = 1205058 19 (0.380)coding (57/540 nt) BW25113_RS06005 phage tail protein
Rejected: Frequency below/above cutoff threshold.

GCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAG  >  NZ_CP009273/1203172‑1203347
                                                                                         |                                                                                      
gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                                        <  2:1195116‑M1/90‑1 (MQ=255)
  tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                                      <  2:39751‑M1/90‑3 (MQ=255)
     tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaa                                                                                   <  1:950076/90‑1 (MQ=255)
     catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                                                   >  2:120697‑M1/70‑90 (MQ=255)
     catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                                                   >  2:1503243‑M1/70‑90 (MQ=255)
     catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                                                   >  1:1557914‑M1/70‑90 (MQ=255)
      atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                                                  <  1:1643736‑M1/22‑1 (MQ=255)
      atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                                                  <  2:1467659‑M1/22‑1 (MQ=255)
           gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                             >  1:913419/1‑90 (MQ=255)
           gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                             >  1:1113033/1‑90 (MQ=255)
           gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                             >  2:1312485/1‑90 (MQ=255)
             cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                                           <  1:1069665‑M1/90‑14 (MQ=255)
                 ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                       <  1:898983‑M1/33‑1 (MQ=255)
                  cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                      >  2:934625‑M1/57‑90 (MQ=255)
                   aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                                     <  2:1879370/90‑1 (MQ=255)
                    aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                                                    >  2:1734611/1‑90 (MQ=255)
                        gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                <  2:879756/90‑1 (MQ=255)
                        gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                >  1:879755/1‑90 (MQ=255)
                        aaaagcaggaatgatgccattttcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                >  2:180116‑M1/51‑90 (MQ=255)
                         tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                                               <  2:715785/90‑1 (MQ=255)
                            ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtt                                                            >  2:986466‑M1/1‑62 (MQ=255)
                               ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                         <  2:1773348‑M1/47‑1 (MQ=255)
                               ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                         <  1:461302‑M1/47‑1 (MQ=255)
                                     aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                                                   >  2:1515550/1‑90 (MQ=255)
                                        tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                                <  1:854227/90‑1 (MQ=255)
                                        ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                                >  2:762473‑M1/35‑90 (MQ=255)
                                         cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                                               <  1:1734610/90‑1 (MQ=255)
                                          ggATGGCACTGTGCAAACGGCTCTCGAAAATCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                              >  1:1546221/1‑90 (MQ=255)
                                          atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                                     >  1:1426655‑M1/33‑83 (MQ=255)
                                          atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                                     <  2:1426656‑M1/51‑1 (MQ=255)
                                            aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt                                            >  1:292152/1‑90 (MQ=255)
                                              ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact                                          <  1:1548647‑M1/90‑47 (MQ=255)
                                              ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                          <  2:913420/90‑1 (MQ=255)
                                              ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                          <  2:1113034/90‑1 (MQ=255)
                                               gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                         <  1:697964/90‑1 (MQ=255)
                                               cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                                         >  1:164315‑M1/28‑90 (MQ=255)
                                                 aCTGTGCAAACGGCTCTCGAAAATCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                       <  2:1546222/90‑1 (MQ=255)
                                                    gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                    >  1:1229623/1‑90 (MQ=255)
                                                      gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat                                  <  2:1146085/90‑1 (MQ=255)
                                                        aaaCGGCTCTCGAAAACCTTGGTTTGTGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                <  2:1628565/90‑1 (MQ=255)
                                                        aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                <  2:1799635/90‑1 (MQ=255)
                                                        aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                <  1:1362880/90‑1 (MQ=255)
                                                        aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                <  2:1489225/90‑1 (MQ=255)
                                                         aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGATATAc                                <  2:452041/89‑1 (MQ=255)
                                                         aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGATATAc                                >  1:452041/1‑89 (MQ=255)
                                                              ttttcataaacgTTGGTTTGTGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                          >  1:102464‑M1/13‑90 (MQ=255)
                                                              ttttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                          >  2:1643737‑M1/13‑90 (MQ=255)
                                                              ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                          >  1:414652/1‑90 (MQ=255)
                                                              ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                          >  1:1810982/1‑90 (MQ=255)
                                                              ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                          >  2:1251659/1‑90 (MQ=255)
                                                               tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                         <  2:491403/90‑1 (MQ=255)
                                                                 tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                       >  1:1773347‑M1/10‑90 (MQ=255)
                                                                 tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga                       >  2:1413762‑M1/1‑25 (MQ=255)
                                                                 tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                       >  2:190774/1‑90 (MQ=255)
                                                                   gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                     <  1:1515549/90‑1 (MQ=255)
                                                                     aaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatg                   >  2:919271‑M1/1‑21 (MQ=255)
                                                                      aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                  <  1:649671/90‑1 (MQ=255)
                                                                      aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                  <  1:1312484/90‑1 (MQ=255)
                                                                       aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                <  2:855867/75‑1 (MQ=255)
                                                                       aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                >  1:855866/1‑75 (MQ=255)
                                                                        ccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg                <  2:1448362/90‑1 (MQ=255)
                                                                        ccTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg                <  2:112056/90‑1 (MQ=255)
                                                                          ttGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                                               >  2:1914776‑M1/1‑16 (MQ=255)
                                                                          ttGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                                               <  1:1914775‑M1/57‑42 (MQ=255)
                                                                          TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa              >  2:898984‑M1/1‑90 (MQ=255)
                                                                           TGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt             >  2:178791‑M1/1‑90 (MQ=255)
                                                                            GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                >  2:426287‑M1/1‑38 (MQ=255)
                                                                            GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                                <  1:426287‑M1/38‑1 (MQ=255)
                                                                            GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                          >  2:1416017‑M1/1‑60 (MQ=255)
                                                                            GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                                          <  1:1416016‑M1/60‑1 (MQ=255)
                                                                              tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc          <  2:1263300‑M1/90‑79 (MQ=255)
                                                                                 gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgc       >  2:188022‑M1/1‑9 (MQ=255)
                                                                                 GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa       >  1:1653077‑M1/1‑90 (MQ=255)
                                                                                 GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa       >  1:388339‑M1/1‑90 (MQ=255)
                                                                                  GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa      >  2:1229624‑M1/1‑90 (MQ=255)
                                                                                     GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa   >  2:887018‑M1/1‑90 (MQ=255)
                                                                                      aaGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttttt  >  1:757149‑M1/1‑4 (MQ=255)
                                                                                      aaGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc                                        >  1:740200‑M1/1‑4 (MQ=255)
                                                                                      aaGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcc                                        <  2:740200‑M1/52‑49 (MQ=255)
                                                                                      AAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATCGAGCAATCAg  <  2:1777962‑M1/90‑1 (MQ=255)
                                                                                         |                                                                                      
GCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAG  >  NZ_CP009273/1203172‑1203347

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.