New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 120326121 (0.700)5 (0.210) 5/136 1.5 26.3% coding (305/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 = 1205058 11 (0.450)coding (57/540 nt) BW25113_RS06005 phage tail protein

GCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAT  >  NZ_CP009273/1204969‑1205142
                                                                                         |                                                                                    
gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                      >  2:131594/1‑90 (MQ=255)
         gCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccc                                                                             >  2:128762‑M2/1‑81 (MQ=255)
             cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac                                                                         <  1:703667‑M2/90‑14 (MQ=255)
              aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                                        >  1:470488‑M2/1‑76 (MQ=255)
              aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                                        >  2:505398‑M2/1‑76 (MQ=255)
                   tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                                   >  1:130953‑M2/1‑71 (MQ=255)
                         tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                             <  1:845767‑M2/90‑26 (MQ=255)
                                                     aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                 <  2:106117‑M2/90‑54 (MQ=255)
                                                      ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                <  2:703667/90‑1 (MQ=255)
                                                          aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                              <  2:832513/72‑1 (MQ=255)
                                                          aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                              >  1:832513/1‑72 (MQ=255)
                                                           cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                           <  2:470528/90‑1 (MQ=255)
                                                             gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc                         >  1:106117‑M2/1‑29 (MQ=255)
                                                                ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATa                                                >  1:540684/1‑64 (MQ=255)
                                                                ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATa                                                <  2:540684/64‑1 (MQ=255)
                                                                        ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTTTTACCAgg              <  2:124925/90‑1 (MQ=255)
                                                                                  aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctta    <  2:863955‑M2/90‑83 (MQ=255)
                                                                                  aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                           <  1:704571/67‑1 (MQ=255)
                                                                                  aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                           >  2:704571/1‑67 (MQ=255)
                                                                                    gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt  <  2:158245/90‑1 (MQ=255)
                                                                                    gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt  <  1:131594/90‑1 (MQ=255)
                                                                                         |                                                                                    
GCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAT  >  NZ_CP009273/1204969‑1205142

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.