New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP000730 = 114451428 (0.810)3 (0.090) 3/248 5.0 9.8% coding (29/108 nt) USA300HOU_1070 hypothetical protein
?CP000730 = 1144528 29 (0.900)coding (43/108 nt) USA300HOU_1070 hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCCCGATTATAAATAACAAAGTAATTACATATGATATAGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGATTTATTTTGAAGAAGAAAA  >  CP000730/1144391‑1144650
                                                                                                                           |                                                                                                                                        
gCCCGATTATAAATAACAAAGTAATTACATATGATATAGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGatgatg                                                                                                                         <  2:420895/141‑1 (MQ=255)
  ccGATTATAAATAACAAAGTAATTACATATGATATAGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGt                                                                                                                       <  1:245281/141‑1 (MQ=255)
     aTTATAAATAACAAAGTAATTACATATGATATAGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGt                                                                                                                    <  1:367978/141‑1 (MQ=255)
                       aTTACATATGATATAGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACaaa                                                                                                  >  1:209468/1‑141 (MQ=255)
                                   ataGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTaacaac                                                                                      >  2:378616/1‑141 (MQ=255)
                                        cTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAAtt                                                                                 <  2:152548/141‑1 (MQ=255)
                                                  aTGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACAt                                                                       <  2:177962/141‑1 (MQ=255)
                                                           ttGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTATCTTCCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGtttgatt                                                              <  1:27095/141‑1 (MQ=255)
                                                            tGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGtttgattt                                                             <  2:258221/141‑1 (MQ=255)
                                                             gAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGtttgatttg                                                            <  2:90847/141‑1 (MQ=255)
                                                             gAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGtttgatttg                                                            >  2:366034/1‑141 (MQ=255)
                                                                 ttAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAAtg                                                        <  1:88359/141‑1 (MQ=255)
                                                                 ttAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAAtg                                                        <  1:88305/141‑1 (MQ=255)
                                                                      gagGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCtt                                                   <  1:317877/141‑1 (MQ=255)
                                                                        ggTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTg                                                 <  1:269765/141‑1 (MQ=255)
                                                                               gTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAaataa                                          <  2:311489/141‑1 (MQ=255)
                                                                                     gTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTATTTAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTTCTTTGAAAATAAtcatca                                    >  2:417998/1‑141 (MQ=255)
                                                                                           aCAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATtata                              <  1:339502/141‑1 (MQ=255)
                                                                                              aTTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCa                           >  1:100254/1‑141 (MQ=255)
                                                                                                 ggTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCAttt                        <  2:327260/141‑1 (MQ=255)
                                                                                                    aTTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGa                     <  1:380298/141‑1 (MQ=255)
                                                                                                        ttAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTg                      >  1:355278/1‑136 (MQ=255)
                                                                                                        ttAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTg                      <  2:355278/136‑1 (MQ=255)
                                                                                                          aGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGatttatt               <  2:198245/141‑1 (MQ=255)
                                                                                                                tGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGATTTATTTTgaag         >  2:333275/1‑141 (MQ=255)
                                                                                                                 gTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGATTTATTTTgaaga        >  2:220519/1‑141 (MQ=255)
                                                                                                                 gTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGATTTATTTTgaaga        >  2:62790/1‑141 (MQ=255)
                                                                                                                       aCCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGATTTATTTTGAAGaagaaaa  <  2:296964/141‑1 (MQ=255)
                                                                                                                           |                                                                                                                                        
GCCCGATTATAAATAACAAAGTAATTACATATGATATAGGCTATAGTTACATGAAATAATTGAATTTAAAGAGGTTGCAGTACTTGTTATGACAATTGGTATTGTTAGCTTTTGTAACGACCATACATGATACCGATGATGGTCGTTTTTTTAATGAACACAAACATGCTAACAACAAATTGCTAAAACATAGTTTGATTTGAATGTGCTTTGAAAATAATCATCAATTATATCATTTTGATTTATTTTGAAGAAGAAAA  >  CP000730/1144391‑1144650

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