New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP000730 = 16014613 (0.040)6 (0.080) 6/262 12.0 67.8% coding (132/177 nt) USA300HOU_1484 hypothetical bacteriophage protein
?CP000730 = 1601491 NA (NA)coding (102/177 nt) USA300HOU_1484 hypothetical bacteriophage protein
Rejected: Coverage evenness skew score above cutoff.

TTATCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAGTTGTTTATGTCCATT  >  CP000730/1601459‑1601627
                                |                                                                                                                                        
ttATCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTtgttgt                              <  2:495257/141‑1 (MQ=31)
  aTCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtc                            >  1:164413/1‑141 (MQ=32)
  aTCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtc                            >  2:320643/1‑141 (MQ=32)
   tCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCCAGAAGTCAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtcc                           <  2:217420/141‑1 (MQ=17)
   tCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtcc                           <  2:574858/141‑1 (MQ=32)
   tCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtcc                           <  2:168255/141‑1 (MQ=32)
    cACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtcct                          <  2:822860/141‑1 (MQ=32)
     aCTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTtcctc                         >  2:379579/1‑141 (MQ=32)
         gCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCttt                     <  1:223640/141‑1 (MQ=32)
          cAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCtttt                    <  1:371682/141‑1 (MQ=32)
            aGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTaa                  <  2:35844/141‑1 (MQ=32)
             gTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAg                 >  1:531684/1‑141 (MQ=32)
                cTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAgttg              >  1:236836/1‑141 (MQ=35)
                cTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAgttg              >  2:734786/1‑141 (MQ=35)
                    cTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAGTTGTTTa          <  2:340597/141‑1 (MQ=37)
                       gTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAGTTGTTTATGt       <  1:304217/141‑1 (MQ=39)
                            cTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAGTTGTTTATGTCCAtt  >  1:552850/1‑141 (MQ=255)
                                |                                                                                                                                        
TTATCACTGGCAAGTCCTACCTTGTCAGCTTGGTAACTGTATTGATTATTTACTTCATTTACTAATCGATTCCATTCGTCTCTTGGTATATCTTCAAGAAGTAAGAGAATCATCTTTAATTTTTCTTCGTTCATTTGTTGTTCCTCCTTTTAAGTTGTTTATGTCCATT  >  CP000730/1601459‑1601627

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